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Predgovor 


Getrenntes zvieder — vereinigen. 
R. Courant 


Ova je knjiga prva u nizu udžbenika iz primijenjene matematike, koje na- 
mjenjujemo studentima kako matematike, tako i tehničkih disciplina. 


Razvoj matematike te fizike i tehnike u našem stoljeću doveo je do daleko- 
sežnog razdvajanja tih dvaju dotada usko povezanih područja. Osim općenite 
i posvuda prisutne težnje za što većom specijalizacijom takvom razvoju uvelike 
je pogodovao tzv. »apstraktni pristup« koji omogućuje da se matematika organizira 
na strogo deduktivnom principu, bez pozivanja na ikakav geometrijski, fizikalni 
ili neki drugi zor. Konsekventna primjena apstraktnog pristupa, a usto i novog 
»matematičkog jezika« u sveučilišnoj nastavi činila je gotovo nemogućim spora- 
zumijevanje prosječnog matematičara s inženjerom pa i fizičarom. Još je veći ras- 
korak na području istraživačkog rada. Opasnosti, koje u sebi krije ovakav razvoj 
dobro je uočio R. Courant u predgovoru glasovitom djelu Courant-Hilbert, Metho- 
den der Mathematischen Physik (Springer, Berlin, 1924), već tada propagira- 
jući potrebu ponovne sinteze. 

Burni samosvojni i u sebi pozitivni razvoj apstraktne matematike nimalo ne 
umanjuje potrebu komunikacije između matematičara i stručnjaka u primjenama. 
Dapače, ta potreba raste iz dana u dan. Namjera autora bila je napisati tekst, koji 
će moći čitati studenti kako matematike tako i fizike i tehnike i koji će pridonijeti 
njihovom međusobnom razumijevanju. 


Izbor materijala slijedi tradicionalni okvir matematičke fizike. Nastojali 
smo da budu zastupljena tri osnovna aspekta primijenjene matematike; formulacija, 
matematičko istraživanje i numerička analiza modela. Od drugih knjiga sličnog 
sadržaja ova knjiga se znatno razlikuje po zastupljenosti prvog aspekta, pa u tom 
smislu predstavlja neku vrstu uvoda u mehaniku kontinuuma. Drugi aspekt je u 
velikom dijelu podvrgnut osnovnoj namjeni, ali se idejno približava i suvremenoj 
teoriji. Posebno označeni dijelovi pisani su strogo; njih čitalac ne-matematičar može 
ispustiti u prvom čitanju, a neke od njih i posve. Čitajući barem neke od tih dije- 
lova on će uvidjeti nužnost strožeg pristupa u iole delikatnijim pitanjima teorije. 
S druge strane, matematički orijentiran čitalac će katkada zapeti na ovom ili onom 
»inženjerski« pisanom mjestu; produži li čitanje naići će na posebno označen tekst, 
u kojem stoji formulacija s uobičajenom strogosti. 


Treći, numerički aspekt ograničen je na metodu konačnih elemenata, osnovno 
sredstvo strukturalne mehanike pa i mehanike. fluida. 


Osnovu prezentiranog materijala čine Zinearni problemi. Nelinearni problemi 
tretiraju se uglavnom kao male smetnje linearnih pa se na primjerima ilustrira 
metoda sukcesivnih aproksimacija. 


Posebnu ulogu u našem izlaganju imaju zadaci; oni čine integralni dio teksta. 
Naučiti studenta da riješi zadatak, izražen žargonom modela, osnovna je svrha ove 
knjige. 

Matematičko predznanje, koje se traži od čitatelja uključuje osnove diferen- 
cijalnog i integralnog računa! te linearne algebre (matričnog računa)?. 

U 1. svesku, koji sada predlažemo čitatelju, obrađuju se jednodimenzionalni 
ravnotežni modeli; oni se svode na rubne probleme za obične diferencijalne jednadšbe*. 


Činjenice o običnim diferencijalnim jednadžbama, koje koristimo, izložene su u 
Dodatku. 


Kod izrade zadataka u 1. svesku pomogli su nam Ž. Hanjš i M. Rogina, 
asistenti Matematičkog odjela Prirodoslovno-matematičkog fakulteta u Zagrebu ; 
Ž. Hanjš je osim toga pomno pregledao cijeli rukopis i dao korisne primjedbe. 
Tekst je s velikim strpljenjem tipkala B. Grdović, administrator Matematičkog 
odjela. Svima njima dugujemo veliku zahvalnost. 


Zahvaljujemo također kolektivu Školske knjige, posebno uredniku Zlatku 
poreru. 


U Zagrebu, siječnja 1984. 


I. Aganović 
K. Veselić 


1 Na primjer Ž. Marković, Uvod u višu analizu I dio, Školska knjiga, Zagreb, 1956; S. Kurepa, 
Matematička analiza, Funkcije jedne varijable, Tehnička knjiga, Zagreb, 1971. 

2 Na primjer S. Kurepa, Uvod u linearnu algebru, Školska knjiga, Zagreb, 1975. 

5% Od domaćih knjiga o običnim diferencijalnim jednadžbama v. na primjer F. Križanić, 
Navadne diferencijalne enačbe in varijacijski račun, Državna založba Slovenije, Ljubljana, 1974; 
V. Marić, M. Skendžić, Obične diferencijalne jednadžbe, Naučna knjiga, Beograd, 1980. 
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OBIČNA DIFERENCIJALNA JEDNADŽBA 


S 1. Ravnoteža žice 


U ovoj glavi proučavat ćemo ravnotežni položaj tanke žice na koju djeluje 
zadana vanjska sila (sl. 1). Neka odsječak [0, 7] na osi-x predstavlja nedeformirani 
položaj žice. Pod utjecajem vanjske sile žica se deformira; 
pretpostavit ćemo da je vanjska sila slaba tako da se žica 
malo deformira, to jest da se njezin deformirani položaj malo 
razlikuje od segmenta [0, /]. Radi jednostavnosti pretpostavit | 
ćemo također da je vanjska sila paralelna ravnini xy; tadai | 
deformirana žica leži u ravnini xy. 

Pri deformaciji točka na mjestu x nedeformirane žice Slika 1. 
prijeđe u neki položaj P(x). Pomak xPG) ima uzdušnu* i 
poprečnu komponentu; pretpostavit ćemo da se uzdužni pomak može zanemariti. Po- 


prečni pomak (progib) točke x označit ćemo sa u (x). Deformirana žica tada ima 
jednadžbu (sl. 2) 


YA 
y=u (2). (1) 
Derivacija u' (x) je mjera deformacije (kraće de- Pix) 
formacija) na mjestu x. Pretpostavka da je de- 
formacija mala izražava se uvjetom rutx) 
0 Ž x 
lu' (x)| < 1 za svako x. (2) Slika 2. 
Iz jednakosti 
u(x)=u(0 + [u (AE (3) 
0 
dobivamo 
luQ&) —uQi=|fv&al<[w&la<[da=x<:, (4) 
0 0 0 
ili 
jeu (o). = MSIE 209) 


* Smjer osi x nazivamo uzdužnim ili longitudinalnim, a smjer osi y poprečnim ili transverzalnim. 
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Prema tome, pri maloj deformaciji relativni progib u(x) — u(0) je (po apsolutnoj 
vrijednosti) malen prema duljini žice. 

Ako se žica pod utjecajem vanjskih sila nalazi u ravnoteži, funkcija u (ravnotež- 
ni progib) zadovoljava diferencijalnu jednadžbu koju ćemo izvesti prihvaćajući ovaj 

Princip ravnoteže sile. Ako je tijelo u ravnoteži, zbroj svih sila koje djeluju 
na bilo koji komad (dio) tijela jednak je nuli. 

Da bismo primijenili taj princip, moramo opisati sile koje djeluju na proizvolj- 
ni dio deformirane žice. Te sile su kontaktna i linijska. 

Unutrašnja kontaktna sila opisuje djelovanje jednog komada žice na drugi 
koji mu je susjedan ; ta sila je potpuno određena položajem kontakta (to jest, ne ovisi o 
veličini komada). Označimo sa q (x) silu ko- 
. . ERA . . . PR 
jom dio P(x) P(1) djeluje na dio P(0) P(x) 
(sl. 1.3); tada dio P (0) P (2) djeluje na dio 
P Ta). P (2) silom — q q(%). Funkcija q (x) de- 
finirana je i na krajevima x=0ix=1; 
q()i—q (0) su vanjske kontaktne si- 


le. Ako je x, < x2, na dio PGP (22) 
Slika 3. djeluje ukupna kontaktna sila 


q(x) —q(x0. (6) 


Neka je i (x) jedinični tangencijalni vektor žice u točki P (x): 


i+u(a)j 


(x) =>. (7) 
VI + (0)? 
Prema (2) vrijedi 
1+(q(x)2al, (8) 
pa ćemo pisati 
i(k)=i+y" (0. (9) 
Ako je deformacija mala, kontaktna sila je tangencijalna na žicu: 
a(a)=a(x)1(); (10) 


a (x) je napetost žice u točki P (x). U daljnjem ćemo pretpostavljati da je žica napeta 
u svakoj točki, tj. da vrijedi 


a(x) > 0 za svako x € [0, 2]. (11) 
Iz (9) i (10) slijedi 


do) =a(x)i+a(a)u' (7 (12) 
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Stavljajući g = q.i + q,], iz (12) dobivamo 
dx (X) = a (2), (13) 
1,(x) =a (x) u (2). (14) 


Prema (13) uzdužna kontaktna sila je napetost žice. Pretpostavka (11) posebno znači 
daje q,(1) >0i g,(0) > 0. 

Vanjska linijska sila je raspoređena po žici. Neka je f = ra PJ # limjska 
gustoća te sile, to jest vanjska sila po jedinici duljine žice. Zbog (8) za element luka 
žice imamo 

ds =)/T + (49)? dx as da; (15) 


zato pod veličinom f možemo razumijevati silu po jedinici duljine nedeformirane 
žice. Ukupna linijska sila koja djeluje na komad P(x,) P(x2) (x, > x2) jednaka je 


[oa=č[10a+7]1,60d5 16 


prvi član na desnoj strani je uzdužna, a drugi poprečna limjska sila. 
Sad možemo primijeniti princip ravnoteže sile. Bit će dovoljno razmotriti 


dio P(x)P(1) za proizvoljno x. Na taj dio prema (6) i (16) djeluju kontaktna sila 
X 


q (x) — mi (0) i linijska sila | 1(#) dć. Ako je žica u ravnoteži, onda je 
0 


qa) —4(0)+[fE)d =0, 17) 
0 
ili u komponentama 
ae) — 4:0) + [f,(Bdć =0, (18) 
0 
dy(2)— (0) + [1, (8) df =0. (19) 
0 


Jednadžba (18) izražava napetost q, (x) = a (x) pomoću vanjskih uzdužnih sila: 
X 

a(2) =a(0) — [f,(E)dč. (20) 
0 


ll 


Iz toga slijedi a(1) =a(0) — | f.(&) dg ili a0) =a(D + f.(£) d&; uvršta- 


vajući to u (20) dobivamo drugi odišk za napetost: 
1 
a(%)=a(D + [f.(8) de. (21) 


Prema tome napetost u točki P(x) jednaka je 
ukupnoj uzdužnoj vanjskoj sili koja djeluje na dio 


PGP (1). Budući da vanjske sile smatramo poz- 
natim, u daljnjem ćemo i napetost a (x) smatrati 
zadanom funkcijom. Posebno za f, = 0 dobivamo 
a(x) = a(l); dakle, ako je žica napeta samo vanj- 
skim kontaktnim silama, napetost je konstanta. Naj- 
jednostavniji način da se žica napne je onaj na 
glazbenim instrumentima; napetost u tom slučaju 
nije dana neposredno, već se mora mjeriti. Zgodan 
način na koji možemo zadati napetost jest da se za 
kraj žice objesi uteg (sl. 4). Ako je masa utega 
M > 0 a g ubrzanje teže (g = 9,81 ms-?), onda 
za napetost imamo a(x) =a(0) =g M. Jedna druga mogućnost opisana je u 
slijedećem primjeru. 


Slika 4, 


1.1. Primjer 


Teška žica (to jest žica na koju djeluje sila teža) slobodno visi obješena za jedan 
kraj i zauzima odsječak [0, 7/7] na osi x. Treba odrediti napetost žice. 

Neka je o linijska gustoća mase žice (to jest masa jedinice duljine); tada je gustoća 
vanjske uzdužne linijske sile (to jest težine) jednako f, = og. Kraj x = I slobodno 
visi pa je g,() =a(/) = 0. Prema (21) dobivamo 


a() =gM(%l), (22) 


gdje je M (x,1) masa komada (x, 7). Primijetimo da u ovom primjeru napetost ne 
zadovoljava uvjet (11). Žica će biti »dobro« napeta ako se kraj x = / optereti utegom 
My > 0; tada je 


: a(2) =gM(s) -+gMo. (23) 


Jednadžba (19) opisuje poprečnu ravnotežu 

žice. Ako se žica nalazi u nekom elastičnom 

0 x sredstvu (sl. 5), onda na nju pored linijske sile 
zadane gustoće f (x) djeluje i linijska sila s gus- 
toćom koja je suprotna progibu, to jest oblika 


Slika 5. —b(x)u (2); (24) 
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ovdje je b (x) > 0 zadani koeficijent elastičnosti sredstva. Umjesto (19) tada imamo 


1()—a0+[G1B—BĐ)u(&)dć=0, (25) 
0 
gdje smo radi kratkoće stavili 
g»=4 = (26) 


Vezu (poprečne) kontaktne sile i (poprečne) deformacije daje jednadžba (14), 
koja se zove zakon ponašanja; u novoj oznaci on glasi 


g(x) =a (xu (2). (27) 


Uvrštavajući (27) u (25) dobivamo 


a(x)u(x)—a(0)uw (0+ ij F6 —b(BĐu(E) dg =0. (28) 
0 


To je integralni oblik zakona ravnoteže (integralna jednadžba ravnoteže). 
Derivirajući (28) po x dobivamo 


(a(x)u (x +F() —b(x)u(x) =0, (29) 
ili 


—(a()u (GY +b()u() = f(2). (30) 


To je diferencijalni oblik zakona ravnoteže (diferencijalna jednadžba ra- 
vnoteže) ili kraće jednadžbe ravnoteže. Iz (30) integriranjem dobivamo (28). 
Drugim riječima, jednadžbe (30) i (28) su ekvivalentne. 

Jednadžba ravnoteže (30) je obična linearna diferencijalna jednadžba drugog 
reda za funkciju u; a i b su koeficijenti, a f slobodni član ili desna strana jednadžbe. 
Jednadžba je homogena ako je f = 0 a nehomogena ako je f s 0. Funkcija u koja za- 
dovoljava jednadžbu ravnoteže zove se ravnotežno (stacionarno) stanje ili 
ravnotežni položaj (progib). 

1.2. Zadaci 


1.2.1. Napišite jednadžbu ravnoteže uzdužno deformiranoga cilindričnog štapa 


(sl. 6). 
ne E x 


Slika 6. 


Rješenje. Pri maloj uzdužnoj deformaciji poprečni presjek štapa pomakne se 
od nedeformiranog položaja x u deformirani položaj P (x) = x + u (x); u (x)je po- 


13 


mak, a u' (x) deformacija na mjestu x. Neka je q (x) kontaktna sila na mjestu x, to 
jest (uzdužna) sila kojom komad (P (x), P(1)) djeluje na komad (P (0), P(x)). 
Zakon ponašanja (Hookeov zakon) ima oblik (27), gdje je a = ES; ovdje je E > 0 
Youngov modul materijala, a S površina poprečnog presjeka. Za homogeni materijal 
E je konstanta. Neka je f gustoća vanjske longitudinalne linijske sile koja djeluje na 
štap. Iz principa ravnoteže sile slijedi jednadžba (30) (sa 2 = 0). 


1.2.2. Napišite jednadžbu ravnoteže štapa kružnog presjeka Koji je podvrgnut 
torziji oko svoje osi (sl. 7). 


Rješenje. Pri maloj zorgionoj defovmaciji poprečni presjek štapa na mjestu 
x zakrene se oko osi x za kut u (x); u (x) je deformacija na mjestu x. Uzrok zakreta 
je moment (kao što je uzrok pomaka sila). Neka je q (x) kontaktni moment na mjestu 
x, to jest (uzdužni) moment kojim dio (x, 7) djeluje na dio (0, x). Zakon ponašanja 
ima oblik (27), gdje je a = uR*7]2; ovdje je u > 0 modul smicanja materijala, a R 
radijus štapa. Za homogeni štap u je konstanta. Neka je f gustoća vanjskog (uzduž- 
nog) linijskog momenta koji djeluje na štap. Princip ravnoteže momenta glasi 
ovako: Ako je tijelo u ravnoteži, zbroj svih momenata koji djeluju na bilo koji komad 
tijela jednak je nuli. Iz toga slijedi jednadžba (30) (sa b = 0). 


1.2.3. Napišite jednadžbu stacionarnog provođenja topline kroz štap. 


Rješenje. Označimo sa u (x) temperaturu poprečnog presjeka štapa na mjestu 
x e (0, 7). Neka je q (x) kontaktni tok (fluks) topline na mjestu x, to jest količina to- 
pline koja u jedinici vremena prijeđe iz dijela (x, 7) u dio (0, x). Zakon ponašanja 
(Fourierov zakon) ima oblik (27), gdje jea = S; ovdje je x > 0 koeficijent provođenja 
materijala, a:S površina poprečnog presjeka. (Fourierov zakon odgovara našemu 
svakidašnjem iskustvu po kojem toplina uvijek prelazi s toplijeg na hladnije, to 
brže što temperatura jače varira.) Neka je f gustoća limjskog toka topline, to jest. 
količina topline koja se (izvana) prenese na jedinicu duljine štapa u jedinici vremena. 
Princip stacionarnog provođenja topline glasi ovako: Ako je provođenje top- 
line kroz tijelo stacionarno, ukupni toplinski tok (toplina po jedinici vremena) koji se 
prenosi na bilo koji komad tijela jednak je nuli. Tz toga slijedi jednadžba (30) (sa b = 0). 


S$ 2. Rubni uvjeti 


Uz zadane funkcije a, b i f jednadžba ravnoteže (1.30) ima opće rješenje u koje 
ulaze dvije proizvoljne konstante; prema tome, ona ima beskonačno mnogo rješe- 
nja. Nas zanima rješenje koje predstavlja ravnotežni progib žice u zadanim vanjskim 
uvjetima. Osim podataka a, b i f tim uvjetima pripadaju i poprečne sile na krajevima; 
one se (posredno ili neposredno) opisuju tzv. rubnim uvjetima. 


Ako je na kraju x = 0 zadan progib ua, onda je 
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u(0) = u. (1) 


[0] x 0 X 
Slika 1. 


(U tom slučaju poprečna sila na kraju nije zadana, nego se očituje kao reakcija fiksi- 
ranog položaja.) Ako je ug = 0, kažemo da je kraj učvršćen (sl. 1). Najjednostavniji 
način da se kraj učvrsti jest onaj na glazbenim instrumentima. Ako je za kraj obješen 
uteg mase M > 0 (kojim se realizira napetost) kao na sl. 2, onda je taj kraj nužno 
učvršćen. Naime, pretpostavka daje vanjska popre- 
čna sila slaba znači da je ona mnogo manja od nape- 
tosti (v. $ 5), u našem slučaju od težine utega; 
zbog zakona ravnoteže ta sila ne može podići uteg. 

Ako jena kraju x = 0 zadana poprečna kontakt- 
na sila q,, onda je 


1(0) = qo5 (2)“ 
uzimajući u obzir (1.11) i (1.27), umjesto (2) pišemo Slika 2. 
do 
f = =>, 3 
“O == (3) 


Ako je do = 0, kraj je slobodan; u tom slučaju imamo 

u' (0) =0. (3) 
Prema tome, tangenta u slobodnom kraju žice paralelna je s osi x (primijetimo da je 
to posljedica pretpostavke a(0) > 0). Uvjet (3) (uz go > 0) može se realizirati 


tako da se kraj veže za kotačić koji se može slobodno kotrljati u poprečnom žlijebu 
i da se za kotačić veže uteg težine g, (sl. 3). 


q,:0 


Slika 3. 
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Ako je kraj x = 0 elastično vezan _sl. 4), yi 
onda je 


d(0) — xu (0) = 0, (4) 


gdje je x > 0 zadani koeficijent; uzimajući opet L 
u obzir (1.11) i (1.27), umjesto (4) pišemo 


(0) —ku(0) =0, (5) Slika 4. 


gdje je & = x/a (0). 
Isti tipovi rubnih uvjeta definiraju se za kraj x = 1. 
Poopćenje uvjeta (1), (3) i (5) su rubni uvjeti: 


at" (0) — Bu) =c (6) 
Put (D+8u()=d, (7) 
gdje su a, B, y, 8, ci d zadani brojevi koji zadovoljavaju nejednakosti 
a>0,8>0a+/B>0, (8) 
9>0,8>0y+8>0. (9) 


Ako je a =0, rubni uvjet (6) je geometrijski ili Dirichletov ; ako je a # 0, uvjet 
(6) je prirodni ili Neumannov. Učvršćenom kraju odgovara Dirichletov, a slobod- 
nom i elastično vezanom Neumannov rubni uvjet. Uvjet (6) je komogen ako je c = 
= 0, a nehomogen ako je c #£ 0. 


Problem nalaženja rješenja jednadžbe (1.30) koje zadovoljava uvjete (6) i (7) 
zove se ravnotežni rubni problem. Očekujemo da će uvjeti (6) i (7) biti dovoljni 
za jednoznačno određivanje proizvoljnih konstanta u općem rješenju jednadžbe 
(1.30). Drugim riječima, očekujemo da rubni problem ima jedinstveno rješenje. 


2.1. Zadaci 


2.1.1. Interpretirajte rubne uvjete (1), (3) i (5) u slučaju (7) longitudinalne de- 
formacije štapa, (72) torzione deformacije štapa i (27) stacionarnog provođenja top- 
line. 


Rješenje. (7) Uvjet (1) označuje da je presjek x = 0 pričvršćen za krutu sti- 
jenku; uvjet (3) odnosno (5) označuje da je presjek slobodan odnosno elastično 
vezan (sL 5). (22) Uvjet (1) označuje da je presjek x = 0 pričvršćen za krutu stijenku ; 
uvjet (3) odnosno (5) označuje da je presjek x = 0 slobodan odnosno elastično ve- 
zan (sl. 6). (771) Uvjet (1) označuje da se presjek x = 0 održava na temperaturi 0“ 
(na primjer, pomoću toplinskog spremnika); uvjet (3) označuje da je presjek top- 
linski izoliran; uvjet (5) označuje da se na presjeku nalazi termički regulator. 
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Slika 5. 


Slika 6. 


Slika 7. 


2.1.2. Homogena teška žica napeta je horizontalno utegom mase M > 0* 
na kraju x = 0. Odredite ravnotežni progib ako j šk drugi kraj (2) učvršćen (sl. 7), 
(2) slobodan (sl. 8). 


Slika 8. 


Rješenje. Zbog homogenosti linijska gustoća mase o je konstanta. Vanjska 
sila je težina pa je f = o g. Iz jednadžbe 


—g Mu =2o0g (10) 

dobivamo 
ug) =—4* +0. 1 
u(%) = — +Cix +C,, (12) 


* Ovdje kao i u daljnjem pretpostavljamo da je težina žice mnogo manja od težine utega 
kojim je ona napeta. 
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gdjesu C, i C, konstante. (i) Iz uvjeta u (0) = u (2) = 01(12) nalazimo C, = ol/2M, 
C, = 0, 


__0 
u (x) =5MĆl-» (13) 
(ii) Iz uvjeta u (0) = v' (2) = 0, (11) i (12) nalazimo C, = o1/[M, C, = 0, 
cE 1 
u (x) =sMmd — x). (14) 


2.1.3. Teška homogena žica napeta je koso utegom mase M > 0 na kraju x = 
= 0. Odredite ravnotežni progib ako je drugi kraj slobodan (sl. 9). 
Rješenje. Imamo 


fx (0) = eg cosp,f () = f, ()= og sin 9. (15) 


Slika 9. 
Prema (1.20) napetost je 
a(x) = Mg > (1 —x)og cos g. (16) 
Jednadžba ravnoteže glasi 
(Mg + —2x)ogcosp)u (x) = —ogsiny. (17) 


Iz toga i uvjeta y' (7) = 0, u (0) = 0 dobivamo 


ća) Lo egsinp *(1— x) 
drogi Mg +(1—x)0g00sq? no 


“ ogsino *(Z — x) dx 


0 
M sino x 0C0Sp 
7 o cosžp m ( "M ee 
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2.1.4. Teška homogena žica napeta je horizontalno utegom mase M na kraju 
x = 0, a nalazi se u homogenom sredstvu s koeficijentom elastičnosti b. 


Odredite ravnotežni progib ako je drugi kraj slobodan. 
Rješenje. Jednadžba ravnoteže je 


"_ku= SQ 
u u M (20) 


gdjejek = VoJMg. Partikularno rješenje je konstanta o g/b. Opće rješenje homogene 
jednadžbe je A ch kx + B sh kx, pa opće rješenje jednadžbe (20) glasi 


u(*) = A ch kx + B sh x + 2£. (21) 


Iz uvjeta u (0) = v' (7) = 0 dobivamo 


A=—t. Boina (22) 


Ravnotežni progib je 


u(x) = 


(1 Se) (23) 


KE 
& 3. Postavka rubnog problema 


Budući da u našim razmatranjima upotrebljavamo postupke diferencijalnog i 
integralnog računa, moramo pretpostavljati da funkcije a, b, f i u posjeduju »dovolj- 
nu glatkost«. U tom pogledu integralni i diferencijalni oblik zakona ravnoteže nisu 
potpuno ekvivalentni. Poći ćemo od integralnog oblika (1.28), jer on stavlja slabije 
zahtjeve na pomenute funkcije. Najjednostavnije je pretpostaviti da funkcija a 
ima neprekidnu prvu derivaciju, da su funkcije b i f neprekidne i da funkcija u ima 
neprekidnu prvu derivaciju (v. zadatke 2.2—2.4). Iz jednadžbe (1.28) odmah 
slijedi (v. 6 1) da funkcija a (x) v' (x) ima neprekidnu prvu derivaciju, to jest da 
funkcija u ima neprekidnu drugu derivaciju te da zadovoljava diferencijalnu jed- 
nadžbu (1.30). Te pretpostavke su međutim preuske za slučajeve koji se pojavljuju 
u praksi. Često funkcije a', b i f imaju u konačnom broju točaka skokove (prekide 
prve vrste, sl. 1). Integralna jednadžba ravnoteže iu tom slučaju ima smisla, čak 
ako pretpostavimo da i derivacija funkcije u ima u konačnom broju točaka skokove. 


Slika 1. 
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Iz jednadžbe odmah slijedi da je funkcija a (x) u' (x) neprekidna (jer je integral fun- 
kcije koja ima konačan broj skokova neprekidna funkcija gornje granice), a to znači 
da je i funkcija u' (x) nužno neprekidna. Međutim, u točki x u kojoj neka od fun- 
kcija a', b i f ima skok, funkcija u (x) nema drugu derivaciju (u"“ (x) ima skok), 
pa u takvoj točki nije zadovoljena diferencijalna jednadžba ravnoteže. Umjesto toga 
imamo samo uvjete neprekidnosti progiba i njegove derivacije: 


u(ko—0)=u(xo +0), (1) 
u (ko — 0) = u (x9 +0). 2) 
3.1. Primjer 


Teška žica sastavljena je od dvaju homogenih komada (0,x9) i (9, 1) (s linijsk-. 
im gustoćama mase 0, i 02) i napeta horizontalno utegom mase M > 0 na kraju 
x = 0. Treba odrediti ravnotežni progib ako je drugi kraj učvršćen. 


Gustoća vanjske linijske sile ima u točki x, skok i dana je formulom 


p Pi 8 0 < X < Xa 
HOJE ra de dejt (3) 
Prema tome, jednadžba ravnoteže glasi ovako: 
—gMu" (x) = 0g,% € (0, x) (4) 
—gMu" () =08x € (1). (5) 
Uzimajući u obzir uvjete u (0) = u (7) = 0 iz (4) i (5) dobivamo 
v&=-%+06 (6) 
01%? 7 
u (x) = — ZM + Cx (D) 
za x € (0, x0) i 
. x 
"&)=-5+D (8) 
u E zibeseje set 
uQ&)=—Me /?2+DE-) (9) 


za x € (x, 1); ovdje su € i D konstante. Uzimajući u obzir uvjete (1) i (2), iz (8) 
—(9) dobivamo za € i D sustav 
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C-p=izay, (10) 


e kT (1) 


Iz toga nalazimo 


ZM MI u e 
l = 
pk ef os 


2 1 om ik: 
- M (Er MI + OM 10) € (0, 5) 


u (x) = 
Ena a [01 mem gle (14) 
SME 1) + (657 MI 4) 1), x € (X, 1). 
Uzmemo li (u (14) limes za o» >0, a zatim / —> oo, dobivamo 
2 
_ SK Oo # €(0%0) 
u(x) = (15) 
01% 
IM 2 X € (0, X0). 
X 
x 
Mg 
y 
Slika 2. 


Usporedivši to s rezultatom zadatka 2.1.2, vidimo da se slobodni kraj može realizirati 
i tako da se zadana žica veže za dodatnu dugačku žicu zanemarive mase (sl. 2,3). 


kO 


Slika 3. 
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* Da bismo preciznije postavili rubni problem, uvest ćemo pojam po dijelovima 
glatke funkcije. Funkcija definirana na segmentu [0, 7] je po dijelovima neprekid- 
na, ako je neprekidna svuda osim u konačno mnogo točaka intervala (0, /) u kojima 
ima prekid prve vrste i u kojima je neprekidna zdesna (sl. 4). Skup svih takvih 
funkcija označavat ćemo sa 3% (0, 2) (kraće sa 5%). Funkcija definirana na segmentu 
[0, 7) je po dijelovima neprekidno diferencijabilna, ako je neprekidna svuda i 
ako ima neprekidnu derivaciju svuda osim u konačno mnogo točaka intervala (0, 7) 
u kojima derivacija ima prekid prve vrste (sl. 5). Skup svih takvih funkcija označa- 
vat ćemo sa 271 (0, 7) (kraće sa 271). Slično se definira skup 9" (0, 2) svih funkcija 
kojima je &-ta derivacija po dijelovima neprekidna. 


* 
i 


mbe 


Slika 4. Slika 5. 


Ako je v€ sf! (0,7), onda je v' očito restrikcija neke funkcije iz 9 (0, 2), 
koju također označujemo sa v' i zovemo derivacija funkcije v na segmentu [0, 7]. 
Slično postupamo i s &-tom derivacijom funkcije iz 3#* (0, 7). 

Navest ćemo pregled osnovnih svojstava skupa 2% (0, 1). 

(1) Funkcija v e 36" ima sve derivacije reda nižeg od k neprekidne na [0, /]**. 

(i) Ako je v e", onda je Ve "(m < k). 


(i22) Akojev, we", ondajev+weH",v we"; iz toga posebno slijedi 
da je 2" linearni prostor. Vrijede formule 


sa . m 
(v + 20) m) — v m) + zu m), (V m 20) m) — > (1 POSL k) WU). (16) 


r"=0 
(tv) Ako je v € 21, onda je 


[7 (x) dx = v(x2) — V (21). (17) 
(v) AkojeveH" i 


x 


vw (x) = il v (£) dć, (18) 


0 


ondajeweHliwW =v. 
) 
* "Tekst koji slijedi (do točke 3.2) može se ispustiti kod prvog čitanja. 
** Pod nultom derivacijom razumijevamo samu funkciju. 
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(vi) Za v, w € #7" vrijedi formula parcijalne integracije 


[voeod=(x)we)—ofk)o)— [vo)wG)dw (1) 


(vii) Za v EH! vrijedi nejednakost 
VX) — v(x)| € |x, — x2| sup |v' (91. (20) 
xe[0,/] 
U daljnjem ćemo pretpostavljati da je 
ae (0, 1),b e (0, 1),f e" (0, 1). (21) 


Polazeći od integralne jednadžbe ravnoteže (1.28), pretpostavit ćemo da je u € 
€ 41 (0, 7). Kao što smo zaključili na početku ovog paragrafa, iz toga odmah slijedi 
da jeu €? (0,17). Diferencijalna jednadžba ravnoteže narušena je u točkama pre- 
kida funkcija a', b i f; u tim točkama ona nema drugu derivaciju. Taj »nedostatak« 
u konačnom broju točaka ćemo zanemariti i reći da funkcija u zadovoljava diferen- 
cijalnu jednadžbu ravnoteže na segmentu [0, 1]. Uz takav dogovor jednadžbe (1.28) i 
(1.30) postaju ekvivalentne. 


Rubni problem možemo sada postaviti ovako: 


Odrediti funkciju ue? (0,1) koja zadovoljava jednadžbu ravnoteže (1.30) 
* na segmentu [0,1] i rubne uvjete na krajevima, 


3.2. Zadaci 


3.2.1. Teška homogena žica napeta je horizontalno utegom mase M >0 
na kraju x = 0. Na dio (49,1) (0 < x9 < 7) djeluje elastična sila s. koeficijentom 
b = const. > 0. Odredite progib, ako je drugi kraj žice učvršćen. 


Rješenje. Jednadžba ravnoteže glasi ovako: 


—u"' (x) = < >x€(0,x0), (22) 
—u' () +Bug)=>,1e (201) (23) 


gdje je & = (b/gM)"'*. Uzimajući u obzir uvjete u (0) = 0 i u (7) = 0 iz toga dobi- 
vamo 


(8) 
= u + Cas €(0, 0) 


«a-| ME s 
1 (ch kx — cth &/ sh kx) + a (1 — 5%) ze E(Ko 1) 


gdje su C i D konstante. Iz (1), (2) i (24) dobivamo za C i D sustav 
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h &x 
m DE = 109: X MELEE 2B ") 
xoC — (ch k xo —cth klsh ko) D =5% + pz (1 a (25) 


0% o chkx, 


— k(shkx,— = ng 
C — k(sh kto — cth kl ch kxo) D => — gg ga EI (26) 
Iz toga nalazimo 
o ++ (1- 2) o. HE 1 Fat) X 

E. a: EJ) OM; haj A 10 Q1 
Donk ĐEKI ) 

D=2 *0 _ Zočhkx 1 (, _shkxo sh kl 
mla k shkl E S) duG= = Fani 


(28) 


3.2.2. Teška žica sastavljena je od dvaju homogenih komada (0, 9) i (Xo> 2), 
s linijskim gustoćama mase 0; i 02, i napeta je vertikalno utegom mase M > 0. 
Odredite ravnotežni progib ako je u (0) =0,u(71) =d > 0. 


Rješenje. Prema (1.22) napetost je 


o og (o — X) +028(1—x) X < % 
S (alo a 20) 


Jednadžba ravnoteže glasi ovako: 
(Mg + 018 (X0 — 2) + 028 (1 — X0)) u (2) =0, x €(0, x), (30) 
((Mg + 028 (1 — 2))u' (XY =0, xe€ (Xx, 1). (31) 


Uzimajući u obzir uvjete u (0) = 0i u (2) = d, iz toga dobivamo 


[o a 
ln (1 x),s<(0,s9 : 
M + +odl—: 
u(x) = | _ Q1Xo + 02 ( X) (32) 
— Zn (148 9) +4 xe(9D), 
1 028 
gdje su Ci D konstante. Iz (1), (2) i (32) dobivamo 
Ratni 1 Oi Z ia o : 
— — hn [1 — Lo BUG) = 
01 ( M +o1X0 +02 (1 si. ln (1 Neris 2) 
(33) 
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S 4. Princip superpozicije 


Ako funkcije 41, 42, ...> 4, imaju isto područje definicije i ako su A1, A2, .--> An 
brojevi, onda se funkcija 


PE = Au + A2u2 slike + Anua (1) 


zove iznearna kombinacija ili superpogicija funkcija 2,, 42, ..-> 4,3 brojevi 4,, 22, ...> 
An Su koeficijenti te linearne kombinacije. 


4.1. Lema 
Proizvoljna superpozicija rješenja homogene jednadžbe 
(auY —bu=0 (2) 


također je rješenje te jednadšbe. 


4.2. Lema 
Proizvoljna superpozicija funkcija koje zadovoljavaju homogeni rubni uvjet 


au'(0)—Bu(O)=0 (3) 
odnosno 
yu (d)+čud)=0 (4) 


također zadovoljava taj uvjet. 

Gornje leme se nazivaju zajedničkim imenom princip superpozicije. "Taj 
princip kazuje da su jednadžba ravnoteže (1.30) i rubni uvjeti (2.6) i (2.7) linearni; 
kraće, da je rubni problem (1.30), (2.6), (2.7) linearan. Općenitiji oblik principa 
superpozicije daje sljedeća lema. 


4.3. Lema 


Ako je funkcija uy (t=1,2,..., 1) rješenje rubnog problema 


(a) — bu + fi = 0 (5) 
a44(0) — Pu (0) = (6) 
Pu, (2) +8u(l)=d, (7) 


onda je superpozicija 
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u=A u +A uzb + At (8) 


s proizvoljnim koeficijentima Ay, Aa> ..> An» rješenje problema 


(au) — bu + * Afi=0 (9) 
izi 
aw (—Pu(O=XAc (10) 
izi 
vu(D+8u() =ž4 nm (11) 


Dokaz ove leme, koji je posve elementaran, prepuštamo čitatelju. 
Princip superpozicije (linearnost problema) omogućuje jednostavan postupak 
homogenizacije rubnih uvjeta, koji se sastoji u sljedećem. Neka je u rješenje problema 


(auY —bu+f=0 (12) 
au'(0) —Bu(0) = c (13) 
vu (d)+8u(]) =d. (14) 


Neka funkcija v zadovoljava uvjete (13) i (14) i neka je u, = u — v. Uvrstivši u = 
= u, +vu (12), dobivamo 


(aug) — buy +fi=0 (15) 

au, (0) — Bu, (0) =0 (16) 

yu (2) +8 (1) =0, (17) 
gdje je 

li =1+(av') — bv. (18) 


“Tako dobivamo rubni problem (15)—(17), u kome su rubni uvjeti homogeni i koji 
je ekvivalentan problemu (12)—(14) u ovom smislu: funkcija u je rješenje problema 
(12)—(14) onda 1 samo onda, ako je funkcija uy = u — v rješenje problema (15)—(17). 


\ 


4.4. Zadatak 


Odredite linearnu funkciju v koja zadovoljava uvjete (13) i (14). 
Rješenje. Neka je 


va)=ilx+8. (19) 
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Uvrstivši (19) u (13) i (14), dobivamo za Č i 9 linearni sustav 
ačt—B8=c (20) 
(py +10) 0+89 =d. (21) 
Zbog (2.8) i (2.9) determinanta tog sustava 
l: A=a8+By+1B8 (22) 
je pozitivna. Dobivamo 


2 cdo+dgf 


Ex . de=c9F19) 
pe g aš Te 


U daljnjem ćemo pretpostavljati da su rubni uvjeti homogeni. 


$ 5. Slučaj b=0 


Ovdje ćemo se zadržati na problemu 


(ay +f=0 (1) 
aw (0) —BPu(0O) =0 (2) 
vu (d)+du(d) =0. (3) 


Iz (1) slijedi 


«= zo [10 4+ (4) 
.0 
— E 
«o-- jeg [rosa ja 56 +e (3) 


0 
gdje su C; i C, konstante. 


Funkcija (5) jest opće rješenje jednadžbe (1)*. Uvrstivši (4) i (5) u (2) i (3) do- 
bivamo 


aC, = 
a (0) — B C2 == 0, (6) 


* Uz pretpostavke (3.21) funkcija (5) pripada skupu 3#72 i zadovoljava jednadžbu (1) u točkama 
neprekidnosti funkcija a'i f. 
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H 
ša * 
ito+o] 5501+80, dn |/04++ [56 G) a(q)dq (D 


"To je linearni sustav za nepoznanice C, i C»; determinanta tog sustava je 


1 
drosarouasi JEL , 


Razlikujemo dva slučaja: 


(i) 8+ 6 > 0 Tada iz (2.8), (2.9) i (8) slijedi A > 0, pa sustav (6), (7) ima 
jedinstveno rješenje. Prema tome, ako je B -- & > 0, problem (1)—(3) ima jedno 1 
samo jedno rješenje; rješenje je dano formulom (5), gdje je (C,, C2) rješenje sustava 
(6), (7). 

(ii) 8 = 8 = 0. Tada je A = 0, pa sustav (6), (7) ili nema rješenja ili rješenje 
nije jedinstveno. Rubni uvjeti su 


u (0) =0, (9) 
u () =0, (10) 
tj. radi se o slučaju kad su oba kraja slobodna. Zbog a, » > 0 iz (6) i (7) slijedi 


ći ti, (11) 
l 
[EZO dx =0. (12) 
0 


Obratno, ako je zadovoljen uvjet (12), onda je za svako C, funkcija 


x € 
"(9 = - [5[/04+c, (13) 
0 


0 


rješenje problema (1)—(3). Imamo ovaj zaključak: U slučaju B= 0 = 0 problem 
(1)—(3) zma rješenje ako i samo ako je zadovoljen uvjet (12); tada problem ima bes- 
konačno mnogo rješenja i ona su dana formulom (13) (to jest, rješenje je određeno do 
na aditivnu konstantu). 

Kako su u slučaju (ii) poprečne kontaktrte sile na krajevima jednake nuli, uvjet 
(12) znači da je ukupna (poprečna) sila koja djeluje na cijelu žicu jednaka nuli; 
kao što znamo od ranije, bez tog uvjeta nema ravnoteže. Ako je taj uvjet ispunjen, 
ravnotežni položaj je određen do na »kruti« pomak cijele žice (sl. 1); zato proizvolj- 
nost konstante C, nema fizikalno značenje. 
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X 


Slika 1. 


Pri izvođenju jednadžbe ravnoteže (8 1) pretpostavili smo malu deformaciju 
Žice. Pokazat ćemo da su formule za rješenje, dobivene u ovom paragrafu, u skladu 
s tom pretpostavkom. Zbog jednostavnosti ograničit ćemo se na slučaj rubnih uvjeta 


u) =y(D=0. (1) 
Iz (4) i (14) dobivamo 
l 
“= 709] #048 (15) 
Neka je š 
aa=min a(x). (16) 
xe[0,1] 
Imamo 
1 I 
115 (8) d&|= | LDL 46 < 1 max FG) (m) 
1z (15)—(17) dobivamo 
1 
lu (| < Ria If (l. (18) 


Ako je max |f| dovoljno malo, desna, a time i lijeva strana u (18) je po volji mala. 
Imamo ovaj zaključak: Akoje vanjska sila dovoljnomala (preciznije: 1 max |f (2)|<a,), 
ravnotežni progib zadovoljava uvjet (1.2) (to jest, deformacija je mala). To svojstvo 
rubnog problema zove se korektnost. 


5.1. Zadaci 


5.1.1. Napišite rješenje problema (1)—(3) u slučaju (i)a =y =0, 
(i) a=6=0. 
Rješenje. (7) Iz (6) i (7) dobivamo C, = 0, 


Lj 


ene ed lEblra dn, (19) 
0 0 


0 


2 LA 


u(*) = — f26|10s+ (jf8y e [roof gy (20 
0 


(i) Iz (6) i (7) dobivamo C, = 0, 


l 
G= J f (1) dn, (21) 


u) = - [se ZO) | f(n) dn + | f(n) %, 55 (22) 


5.1.2. Postavite i riješite rubni problem (u slučaju b = 0) ako je kraj x = 0 
učvršćen, a kraj x = I oslonjen (sl. 2). 


y 


Slika 2. 


Rješenje. Oslonjen kraj se ne može opisati uvjetom oblika (3); za njega postoje 
dvije mogućnosti: ili je u kontaktu s osloncem (sl. 3), ili nije (sl. 4). U prvom slučaju 
jeu(l) =0,gq(1) =a(D)v(1) >0,audrugomu (1) > 0, g(D)=a (1) vu (1) =0. 
Prema tome, u oba slučaja je 


uQ)>04(0)>0u(Du(1)=0. (23) 


y 


Slika 3. Slika 4. 


Taj uvjet je jednostran (unilateralan) i bitno se razlikuje od uvjeta (3) koji je dvostran 
(bilateralan). Uvjet (23) je nelinearan, jer za njega ne vrijedi princip superpozicije. 
Moramo razmotriti dvije mogućnosti. 

(1) u(l) =0. Rješenje je dano formulom 20,4 a nužno je v' (2)>0, to jest 


Fade GU) Ja ajina=o (24) 
0 0 
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(2) u' (1) = 0. Rješenje je dano formulom (22), a nužno je u(2)>0, to jest 


t č U H 


- [29 [104n+ [roman [7520 (25) 
Neka je i l : 
. l 
L= [so s [rnan + ro dy [6 [6 [ra dn = (26) 
- [o SEO 


Iz (24) i (25) zaključujemo: Ako je L < 0, rješenje je dano formulom (20); ako je 
L >0, rješenje je dano formulom (22). 


Ako je a = const., imamo 
1 


aL=([nf(dq: (27) 
0 


to je ukupni moment vanjske sile u odnosu prema polu x = 0. 


8 6. Slučaj b # 0 


Ako je b 5 0, rješenje rubnog problema 


—(auY +bu=f (1) 
aw()—fBu(0) =0 (2) 
yu ()D+8u(]) =0 (3) 


općenito se ne može napisati u zatvorenom obliku (poput onoga u slučaju b = 0). 
Najprije ćemo dokazati jedinstvenost rješenja tog problema. 


6.1. Teorem 
Ako je b #£ 0%, rubni problem (1)—(3) ima najviše jedno rješenje. 


Dokaz. Neka su u, i u» rješenja problema i neka jeu = #, — u,. Zbog linear- 
nosti problema funkcija u zadovoljava uvjete 


* Pretpostavljamo, naravno, da je b (x) > 0 za svako x. 


Sl 


(auY — bu =0, (4) 
aw (0) —Bu(0) =0, (5) 
yu () +8u(1) =0. (6) 


Pomnožimo li jednadžbu (4) sa u i integriramo od 0 do /, dobivamo 
1 ! 
— | (au'Y u dx + | bu? dx = 0. (7) 
0 0 


Prema formuli parcijalne integracije vrijedi 


1 1 1 
| (au'Y u dx = | (au' uY dx — d avudx =a(Duv(Du(dl)— 
0 0 0 (8) 
gd 
_ a(0)w(0)2(0)— [ a (u)? dx. 
0 
Iz (7) i (8) slijedi 


1 
| (a(4)2 + bi2)dx +a(0)r (O+) —aDu(QD)u(l) =0. (9) 
0 


Razlikujemo četiri slučaja : 
(1) B > 0, y > 0. Izračunamo li u (0) iz (5) i u' (2) iz (6) i uvrstimo u (9), do- 
bivamo 
1 


1 (a (u)? + bu?) dx + 


0 


a 


52 (0) (1 (0+ a(Du(0 — 0, (10) 


Svi sumandi u (10) su nenegativni pa je svaki od njih jednak nuli: 


l 


[a (u)? dx =0, (11) 
0] 
H m 
[o dx =0, (12) 
0 
2 (0) (0% — 0, (13) 
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Za (Duž (1) = 0. (14) 


Zbog a (x) > 0 iz (11) slijedi “' = 0, to jest 
u = const. (15) 


Iz (5)1(15) slijedi u (0) = 0,azatimiz(15)u =0tj. u =u,. 
Slučajevi (#) B > 0, 8 >01i (222) a > 0, d > 0 analiziraju se analogno. 


(2) B=8 =0. Tada je u' (0) = u (2) = 0, Iz (9) slijedi (11) i (12); iz (11) 
slijedi (15), a iz (12) slijedi 


u? | bdx =0. (16) 
0 


Kako je po pretpostavci b (x) > 0 i b # 0, b je pozitivno bar na nekom intervalu 
pa je 
l 
[ bdx > 0. (17) 
0 


Iz (16) i (17) slijedi u =0, tj. 4 = u. 


U gornjem dokazu u slučajevima (7)—(122) nismo se koristili pretpostavkom 
b #0. Ti dokazi, dakle, vrijede i za slučaj b = 0 (85). Pretpostavka b # 0 igra, 
međutim, odlučnu ulogu u slučaju (29); prema $ 5. znamo da za b = 0 jedinstvenost 
ne vrijedi. 


$ 7. Svojstvene vrijednosti i svojstvene funkcije 


Kao što smo vidjeli u & 6, za jedinstvenost rješenja problema (6.1)—(6.3) važna 
je pretpostavka b (x) > 0. S druge strane može se očekivati da će se (uz izuzetak 
rubnih uvjeta u' (0) = u' (/) = 0) jedinstvenost sačuvati i za negativno, ali po ap- 
solutnoj vrijednosti ne preveliko b. "To očekivanje podupire činjenica da o jedinstve- 
nosti odlučuje nenegativnost izraza na lijevoj strani u jednakosti (6.9); izraz će 
ostati nenegativan ako b postane malo negativno. 

Najjednostavniji slučaj u kojem se može pojaviti negativno b jest problem rav- 
noteže žice u sustavu koji jednoliko rotira oko osi x kutnom brzinom £2 > 0. Tada 
na žicu djeluje dodatna linijska sila gustoće 


2? o(x)u (x) (1) 
(centrifugalna sila). Tako ukupna vanjska linijska sila ima gustoću 


J(&)—b(2)u(2), (2) 
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gdje je 
b (2) = bo (4) — £2*0 (2), (3) 


a b, (x) > 0 dolazi od vanjskoga elastičnog sredstva. Tako za veliko dobivamo 
negativno 2 (x). 


7.4. Primjer 


Homogena žica napeta je kontaktnim silama (tj. a = const.) i učvršćena na 
krajevima u sustavu koji jednoliko rotira oko osi x kutnom brzinom 2. Treba od- 
rediti kritičvu vrijednost te brzine, to jest onu vrijednost pri kojoj osim nedeformira- 
nog ravnotežnog položaja u = 0 postoji i deformirani ravnotežni položaj u # 0. 
(Tada kažemo da je nedeformirani položaj nestabilan.) 

Neka je 

A = pla. (4) 


Jednadžba ravnoteže u rotirajućem sustavu glasi 


u'+21u=0; (5) 
rubni uvjeti su 


u()=u(l)=0. (6) 


Vrijednost parametra A za koje rubni problem (5), (6) ima nerrivijalno rješenje 
u # 0 zove se svojstvena vrijednost; odgovarajuće rješenje zove se svojstvena 
funkcija. Nas zanima najmanja pozitivna svojstvena vrijednost A1; ona određuje 
kritičnu kutnu brzinu %2, = (24,a/0)'!2. Neka je 2 svojstvena vrijednost i u odgova- 
rajuća svojstvena funkcija; tada je 


u()=Ac0s/Aox+BsinVlox, (7) 


gdje su A i B konstante. Iz uvjeta u (0) = 0 slijedi A = 0; iz uvjeta u (2) = 0 do- 
bivamo 


BsnVlel=0, (8) 
ili (zbog u # 0) 
sin Vlo I =0. (9) 
Iz toga slijedi 
Vžel="nmn=12,.. (10) 
ili 
1-2. (1) 
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Vrijednost A koja odgovara broju 2 označimo sa A,: 


h=h=—=. (12) 
Odgovarajuća funkcija (8) je 


(13) 


u, (X) = B, sin 


gdje je B,, = const. Lako se provjerava da je (12) odnosno (13) zaista svojstvena 
vrijednost odnosno svojstvena funkcija problema (6). Najmanja svojstvena vrijed- 
nost je 


4 = (14) 


odgovarajuća kutna brzina je 


9, «EVE, (15) 
It e 


a odgovarajući deformirani ravnotežni položaj (sl. 1) 


u) =B, sin E. (16) 
y 
: 1/2 1 X. 
Slika I. 


Proizvoljnost konstante B, značila bi da se pri kritičnoj kutnoj brzini žica kida. 
Opažanje, međutim, pokazuje da se pri toj brzini postiže samo relativno velik pro- 
gib. Taj progib se može proračunati iz realističnije nelinearne teorije koja vrijedi 
za veliku deformaciju. 

Činjenica da poslije prekoračenja kritične kutne brzine rješenje postaje opet 
jedinstveno nije važna, jer nelinearna teorija promatrane pojave pokazuje da je od- 
govarajući ravnotežni položaj u = 0 labilan; to znači da je ravnotežu nemoguće 
zadržati. 


Općenito će problem ravnoteže u rotirajućem sustavu biti opisan jednadžbom 
—(a(x)uw (x) +bo(x) u) — Hocu) =f() (17) 
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;rubnim uvjetima 
at (0) —Bu(0) = 0, (18) 
yu(l) +84(D=0 (19) 
O jedinstvenosti rješenja problema (17)—(19) odlučuje homogena jednadžba 
—(au)Y' +bau—2ou=0, (20) 


uz rubne uvjete (18)—(19). Naime, ako homogeni problem (18)—(20) zma samo tri- 
vijalno rješenje u = 0, onda nehomogeni problem (17)—(19) ma najviše jedno vješenje; 
ako homogeni problem ima netrivijalno rješenje ug 5 0 1 ako nehomogeni problem ima 
rješenje u,, onda je 1 funkcija ug -- u, rješenje nehomogenog problema, 


U daljnjem konstantu A smatramo parametrom. Vrijednost A za koju homogeni 
problem ima netivijalno rješenje zove se svojstvena vrijednost problema; 
odgovarajuće rješenje zove se svojstvena funkcija. Zbog linearnosti problema 
svojstvena funkcija je određena do na konstantan faktor. Za vrijednosti /. koje nisu 
svojstvene kažemo da su regularne. Ako je A regularno, onda rubni problem (17)— 
(19) zma najviše jedno rješenje. 


Prema našim zaključcima u $ 6. svako A < 0 je regularno. Ako jeb=01B = 
= 0 =0, onda je A = 0 svojstvena vrijednost, a odgovarajuće svojstvene funkcije 
su konstante. Ako je bar jedan od rubnih uvjeta Dirichletov, vrijednost 2. = 0 je re- 
gularna. 


Svojstvene vrijednosti imaju važnu interpretaciju u vezi s titranjima, o čemu 
će biti govora kasnije. U vezi s jednadžbom ravnoteže svojstvene vrijednosti se po- 
javljuju u problemu stabilnosti. U svakom slučaju posebno je važna najmanja pozi- 
tivna svojstvena vrijednost 4,5 ako takva postoji, sve vrijednosti 0 < 2 < 2, su 
očigledno regularne. 


7.2. Zadaci 
7.2.1. Riješite problem svojstvenih vrijednosti 
u"+2u=0u(0="(7) =0. (21) 
Rješenje. Opće rješenje jednadžbe je 
u(x) = Acos/lx+Bsin/1x. (22) 
Rubni uvjeti daju A =0i 


cos A1 = 0; (23) 


iz toga dobivamo (sl. 2) 


Slika 2. 
(2n +1)? 22 
a (24) 
un) = 8B, snima : (25) 


7.2.2. Dokažite formulu koja povezuje svojstvenu vrijednost A i svojstvenu 
funkciju u problema (18)—(20), a = 6 = 0: 


[ (a (u')2 + buž) dx 
habe, (26) 
[ou dx 
0 


Rješenje. Pomnožimo jednadžbu (4) sa u i integriramo od 0 do /, a zatim prvi 
član na lijevoj strani parcijalno integriramo uzimajući u obzir rubne uvjete. 


7.2.3. Neka je a=86 =0, a, =mina, by = minb, o, = max o. 
Dokažite da su sve vrijednosti 


2 
ia a (7 H bh) (27) 


regularne. 
Rješenje. Iz formule (26) slijedi 


H 1 
af (4)? dx +, fu dx 
Nj. a. o ovna. aja (28) 


— 
1 [uidx 
0 
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Iz jednakosti 


X 


u) = ji uy (£) dć 


0 


dobivamo (koristeći se Cauchyjevom nejednakosti): 


Ga 699% = ( fu 948) < [at fd (BP ae = = [4 Eat < 
0 0 0 0. 


l 
< x | (4 (dg, 
0 


U l 
| (u, (x))* dx < z | (u; (x))2 dx, 
0 Č 


ili 

: i t 
[ (u! (#))2 da > a (eje 
0 


0 


Iz (28) i (32) slijedi 


$ 8. Koncentrirano djelovanje 


(29) 


(30) 


(31) 


(32) 


(33) 


Do sada smo vanjsku silu na žicu opisivali gustoćom f, tako da je ukupna vanj- 


_—_ 
ska sila na luk P(x,)P (x2) bila jednaka 


*2 
ll f(2) dx. 
37 
Stavimo li 
s. 
F (x) = ji Ida ti FG =F() 
0 
= . 
ukupnu vanjsku linijsku silu na luk P(x,) P (x,) možemo izraziti kao 
F(&2) — F(x,). 
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(1) 


(2) 


(3) 


Pokazat ćemo da je formula (3) općenitija nego (1), i to na važnom primjeru kon- 
centrirane sile, koja nema gustoću. 


Kažemo da funkcija F opisuje silu (djelovanje) intenziteta F, #£ 0 koncentriranu 
u točki x, € (0,7) ako je 


0 zax<x, 


F(x) = 4 
& MZax>xy. €) 

Prema (3) ukupna sila na žicu jednaka je 
FQq)—FO=Fi (5) 


S druge strane, uzmemo li bilo koji interval (x,, 2) takav da je udaljenost točke 
Xg Od (X,, x2) veća od nule, dobivamo 


F(&)— F(e)=0. (6) 


"Time se opravdava naziv »koncentrirana« sila. :Pokažimo da za sila nema gustoću. 
Zaista, f (x) postoji svuda osim za x = xy te vrijedi 


I(9) =0, x # x95 (7) 


pretpostavimo li sad f (x) = F' (x), moralo bi biti 


I 
Fo = [f(2) dx, (8) 
0 


što je protuslovlje, jer integral funkcije koja iščezava svuda osim u jednoj točki jed- 
nak je nuli. Dakle, pretpostavljena funkcija f (gustoća) ne postoji. 


= 

Primijenimo princip ravnoteže sile na komad P (x,) P (x2) (x, < x2) uzima- 

jući u obzir da na taj komad osim kontaktne i vanjske elastične sile djeluje i linijska 
sila (3). Dobivamo: 


(2) —a(x)— | budx+F() — F(x) — 0. (9) 
Zbog (4) iz toga slijedi 


x2 
q(x2) — q(x1) — | buda =07zax,,X2 < X 1 ZA Xi, X2 > Xo, (10) 


X1 


d(x) —a(k)— [| budx+ Fo —0 za x, Kay ati: (11) 


E2) 
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Derivirajući (10) po x, i stavljajući x» = x, dobivamo 
q0G%) —dejaci=0mrea, 124 x > ka (12) 
Iz toga i zakona ponašanja q = ar' slijedi 
(a (x) (x) —b(x)u(x) =0 za x < x i Za x > X. (13) 


Drferencijalna jednadžba (13) vrijedi, dakle, ga x # Xx. U točki x, progib je nepreki- 
dan: 


u(xo — 0) = u(x9 + 0). (14) 


Uzmemo li u (+1) limes kad x,, X2, —> o, dobivamo 


100—0)—q(o+0)=F (15) 
ili 
: ; Fo 
Po e EE (16) 


Vidimo da u točki koncentracije vanjskoga linijskog djelovanja derivacija ravnotežnog 
stanja (progiba) ima skok (prekid 1. vrste). Prema tome, u slučaju koncentriranog 
djelovanja za ravnotežno stanje vrijede jednadžbe (13), (14) i (16). Odgovarajući rub- 
ni problem uključuje i rubne uvjete (2.6), (2.7). 


8.1. Primjer 


Neka je a =1,b =0. Odredimo ravnotežno stanje ako je vanjska sila F4 
koncentrirana u točki x € (0,7) i ako je lijevi kraj žice učvršćen, a desni slobodan. 
Imamo ove uvjete: 


uQ&)=0x<x (17) 

u(0) =0 (18) 

u Q&)=0x> x (19) 

w(D=0 (20) 

u(xo—0)—u(xo +0) =0 (21) 

(o—0)—u (x +t0)=1. (22) 

Iz (17) i (18) dobivamo 

uQ)=Axx<xo (23) 

uQ&)=B,x > x, (24) 
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gdje su A i B konstante. Iz (21) i (22) dobivamo A = 1,B = xy. 
Prema tome je (sl. 1) 


Skhika I. 


(25) 


Za razliku od slučaja vanjskog djelovanja s gustoćom (klase 27), 
ravnotežno stanje pod utjecajem koncentriranog djelovanja ne može 
biti klase 272 (0, 7); to se vidi iz uvjeta (16). Prirodan je zahtjev da to 
stanje pripada klasama 2! (0, 7) (time se zadovoljava uvjet (14)), 27% 
(0,9) 1362(x,9,1). Prema tome, precizna postavka rubnog problema u 


ovom slučaju glasi ovako: 


Odrediti funkciju u e 31 (0,1) koja je rješenje jednadžbe (13) na 
intervalima (0, X9) i (X, 1) € koja zadovoljava uvjete (16), (2.6), (2.7). 


8.2. Zadaci 


8.2.1. Za slučaj koncentriranih djelovanja izrecite princip superpozicije i 


ispitajte jedinstvenost rješenja rubnog problema. 


Rješenje. Neka su F, sile intenziteta F;, koncentrirane u točkama x, 2 =1, 


2,...>n. Neka je u, rješenje problema 
(a(k) u; (XY — bu (X) =0, x # x, 
u; (Xi — 0) = u; (2% + 0) = 0 


Fio 
a (x) 


au;(0)—Bu(0) =, 


u —0)—u(w+0) = 


Pu) +8u4(D)=d,. 


Neka su A;> A2, .--> An proizvoljni brojevi. Tada je 


n 
u=x4u 
i=1 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 
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rješenje problema 


(a(X)u' (x) — b(Xu(Q) =0, x # X, 22 ...9 Xu (32) 
u(qq—0)—u(a+0=0i=12,...,n (33) 
“(> 0 — “(+0 — E i bo (34) 
au(0)—Bu(0) =ž4 a; (35) 

yu D+šuD=žAd (36) 


Zbog gornjeg principa superpozicije ispitivanje jedinstvenosti svodi se na pro- 
matranje homogenog problema, tj. slučaja F; = 0 uz homogene rubne uvjete. Pre- 
ma tome, odgovor na pitanje jedinstvenosti je isti kao u slučaju djelovanja s gusto- 
ćom. Posebno, ako je b = B = & = 0 (uz homogene rubne uvjete), rješenje postoji 
onda i samo onda kad je 


žFio =0. (37) 


8.2.2. Teška žica sastavljena je od dvaju homogenih komada jednake duljine 
s linijskim gustoćama mase g, i 02 i napeta horizontalno utegom mase M > 0. 
Odredite ravnotežni progib ako je žica u sredini opterećena utegom mase M, > 0 
i ako je drugi kraj učvršćen (sl. 2). 


"(= _ 5% «< s (38) 
X» 1 
u &=-%.=>5 (39) 
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Iz toga i rubnih uvjeta u (0) = u (/) = 0 dobivamo 


o 1 x2 Cx x 1 
2M A “ 2 
u() = 


gdje su C i D konstante. Iz (14)i(16)za Ci D dobivamo sustav 


C+D=zgla + 3 02) 
a 1 
C—D= T ZM — 02), 
odakle nalazimo 
M 
C = ARO Be 01 + 02), 
M l 
n= =SuraEm o +502). 


S 9. Greenova funkcija 
Ovdje ćemo pretpostavljati da homogeni problem 
(au) —bu=0 
aw(0)—pBu(O) =0 


vt ()D+8u(l) =0 


BR —M)+DGe—D mu 
| 2M“ u LE 


(40) 


(41) 


(42) 


(43) 


(44) 


(1) 
(2) 
(3) 


ima samo trivijalno rješenje u = 0 (pretpostavlja se da je b(x) > 0). Prema $5 


i 86 to je ekvivalentno pretpostavci da je 
iib£0ii B+86>0. 


t 
Slika I. 


(4) 


43 


Razmatrat ćemo slučaj kad na žicu djeluje vanjska jedinična sila koncentrira- 
na u točki x" € (0, 7). Ravnotežni progib koji zadovoljava rubne uvjete (2) i (3) oz- 
načit ćemo sa G (x, x"). Ako ovdje i x' smatramo promjenljivim, dobivamo funkciju 
dviju varijabli definiranu za sve x, x' € (0,7) (sl. 1). Tako definirana funkcija G 
zove se Greenova funkcija rubnog problema (1.30), (2.6), (2.7); prema & 8 ona 
je rješenje sljedećeg problema. 


GOG) —bO GG) =00<x<x<I (5) 
45-000) —8G(0x)=0, x'€(0,0) (6) 
zGOZGEx) —bG)GEx)=00<x<x<1 7 
950 (4x) +86 (1x) —0, x e(01) (8) 
G(x'—0,1)—G(x' +0,2)=0, 1 €(0,1) (9) 

E Gia 2 ZB (5:9 (a; 9) =. 9 sel0D. (0 


9.1. Primjer 


Odredimo Greenovu funkciju u jednostavnom slučaju a = 1, b =0, kad je 
lijevi kraj žice učvršćen, a desni slobodan. Prema 8.1 imamo 


x XxZ£x 


G (x,x) = (11 
em) E x>x'. ) 
Kraće možemo pisati 
G (x, 2") = min (x, x), (12) 
6dakle slijedi da je funkcija G simetrična: 
G(x,x/)=G(x,2x). (13) 


Prema tome zrogib u točki x izazvan jediničnom silom u točki x' jednak je progibu 
u točki x' izazvanom jediničnom silom u točki x. Iz (12) slijedi da je funkcija G po- 
zitivna: 

G (xx) > 0 za sve x, x' €(0,/). (14) 


To svojstvo ima jednostavnu interpretaciju: progib prouzročen koncentriranom si- 
lom ima smjer te sile. Svojstva (13) i (14) Greenova funkcija ima i u općem slučaju, 
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9.2. Zadatak 


Neka je a = const., b = const. Odredite Greenovu funkciju ako su oba kraja 
učvršćena. Provjerite simetričnost i pozitivnost, 


Rješenje. Stavljajući & = (bla)"/?, iz (5) i (6) odnosno (7) i (8) dobivamo za 
b>0 


G(xx)=Bshkxx<x, ' (15) 
G(x,x') = C(chkx — cth klshkx), x < x, (16) 


gdje su B i C konstante koje mogu ovisiti o x". Sad iz (9) i (10) slijedi 


Bshkx' — C(chkx' —cthkish kx'") = 0, (17) 
Bchkx' — C(shkx —cthklchkx') S (18) 
Iz toga dobivamo 
B= Z (eh kx' — cth klshkx), (19) 
l , 
C =—shkx, (20) 
ak 
pa je 
| (eh kx' — cthklshkx'")sh kx, x < x" 
G(ms)=1" (21) 
Se (ch kx — cth kl sh kx) sh kx', x < x. 
a 
Simetrija je očita iz (21). Dalje imamo 
! ik l ' 7 _shk(l—x) 
ch kx' —cth kl sh kx' = rz Sh kl ch kx' — ch kl sh kx') = Pimji 0, 
(22) 


pa je stoga G(x,x) >0za0<x<x<1Za0<x<x<1I pozitivnost se 
provjerava analogno. Za b = 0 dobivamo 


(23) 


Simetričnost i pozitivnost je očita, 
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Prelazimo na konstrukciju Greenove funkcije u općem slučaju. Prema 
$ 8. precizna postavka problema (5)—(10) glasi ovako: 

Odrediti funkciju G (x, x") varijabli x, x' € (0,1) s ovim svojstvima: 

(1) za svako fiksirano x' funkcija G (x, x") pripada klasi 6" (0,1); 

(11) za svako fiksirano x' funkcija G (x, x") je rješenje jednadžbe 
(1) na segmentima [0, x'] i [2', 1]; 

(2:4) funkcija G (x, x") zadovoljava uvjete (6), (8) £ (10). 

Pokazat ćemo da tako postavljen problem ima rješenje (koje je 
prema zadatku 8.2.1 jedinstveno). 


Neka su #, i u, linearno nezavisna rješenja jednadžbe (1) na seg- 
mentu [0, 7]; tada svako rješenje ima oblik 


u=A, U +A,1m, (24) 


gdje su A, i A, konstante (v. Dodatak, $ 2). Odredimo sva rješenja koja 
zadovoljavaju uvjet (2). Uvrstimo li (24) u (2), dobivamo 


aiA, + a,A, = 0, (25) 
gdje je 
a = au, (0) — Bu, (0), az = au, (0) — 8 u2 (0). (26) 


Bar jedan od brojeva (26) je različit od nule, jer bi iz uvjeta a, =0, 
a» =0 i pretpostavke a + f > 0 slijedilo 


uy (0) — u, (0) 


s O pb To) = W(u,,u20) =0, (27) 


a to zbog linearne nezavisnosti funkcija u, i u» nije moguće (v. Doda- 
tak, $ 2). 


Ako je ay # 0, onda je 


čl = = Mia (28) 


a, 
u=A, (- 24 +1), (29) 
đi 
pri čemu je (zbog linearne nezavisnosti u, i #2) 


in ha 2th (30) 
4 


Ako je a, # 0, onda je 


a 
A.=— rme (31) 
u=A, (m - a2) (32) 
a2 
pri čemu je 
U — zu #£0. (33) 


Prema (29) odnosno (32) svako rješenje koje zadovoljava uvjet (2) ima 
oblik 


u=Cv, (34) 


gdje je v jedno fiksirano netrivijalno rješenje koje zadovoljava uvjet (2), 
a C konstanta. Analogno zaključujemo da svako rješenje koje zadovo- 
ljava uvjet (3) ima oblik 


u=Duw, (35) 


gdje je zw jedno fiksirano netrivijalno rješenje koje zadovoljava uvjet (3), 
a D konstanta. 

Uzimajući u obzir gornje zaključke vidimo da Greenova funkcija, 
nužno ima oblik 


Cv(x,O0<x<x<1I 


(36) 
Du(x), 0<x<x<I, 


Ga&r)= 


pri čemu koeficijenti C i D ovise o x. Uvrštavajući (36) u (9) i (10), 
za.te koeficijente dobivamo sustav 


Cv(x) — Dv(x) =0 (37) 


Cv(x) —Dw(x) = a : (38) 


x)" 


Determinanta tog sustava W (x) = W (v, w, x') je različita od nule za 
svako x“. Zaista, kad bi W negdje iščezavalo, funkcije v i w bile bi linear- 
no zavisne pa bi obje zadovoljavale i uvjet (2) i uvjet (3); zbog (4) iz 
toga bi slijedilo v = w = 0, što je protivno pretpostavci. Rješenje sus- 
tava (37), (38) je 


w(x) zo (x) 


ĆE ZG)VGJ TO 


(39) 
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BA DAJ oke via) 
2-— E) VG) GOW e 


ovdje smo iskoristili Liouvilleovu formulu (v. Dodatak, $ 2). Svaka 
od funkcija v i w je određena do na konstantan faktor ; te faktore uvijek 
možemo odabrati tako da bude 


a (0) W (0) = —1. (41) 
Uz takav izbor dobivamo 
w B LE ! 
Go) [da 0<x<x<1 (42) 
v(x)w(x,O0<x/<x<I. 


Sad se lako provjerava da funkcija (42) zadovoljava sve uvjete koji se 
zahtijevaju od Greenove funkcije. Imamo, dakle, ovaj teorem: 


9.3. Teorem 


Greenova funkcija rubnog problema (1.30), (2.6), (2.7) dana je 
formulom (42), gdje su v € w netrivijalna rješenja problema (1), (2) od- 
nosno (1), (3) koja zadovoljavaju uvjet (41). 


9.4. Posljedice 


(2) Greenova funkcija je simetrična, tj. vrijedi G(x,x") = G (x', x). 
(i) Greenova funkcija je neprekidna na zatvorenom kvadratu [0,7] x 
x [0,7]. (ii) Greenova funkcija se ne poništava na kvadratu (0,1) X 
x (0,2). 


Dokaz. Tvrdnje (1) i (i) su očigledne. Dokažimo (12). Pretposta- 
vimo da je G (x,, x) = 0 za neko x,, x, € (0, 0); neka je na primjer, 
Xi Z x. Iz (42) slijedi v (x,) w (x1) = 0, pajeiliv (x,) = O ili w(x4) = 
= 0. Neka je, na primjer, 


v(xi) =0. (43) 
Tada je v rješenje homogenog problema (1), (2), (43), pa je v(x) =0 
za sve x € [0, x,]. Zbog toga jei W (x) =0zasvex € [0,x,],a to pro- 
tuslovi svojstvu Wronskijana linearno nezavisnih rješenja (v. Doda- 
tak, $2). 
S 10. Rješenje rubnog problema pomoću Greenove funkcije 


I ovdje pret postavljamo uvjet (9.4). 
Greenovu funkciju pridružili smo tubnom problemu 
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—(auY +bu=f KL) 
au (0) —Bu(0) =0 (2) 
yu D+8u()=0, (3) 


iako znamo ($ 8) da se ona ne može dobiti kao rješenje tog problema ni za jedno /. 
Razlog zašto smo to ipak učinili jest činjenica da Greenova funkcija omogućuje 
jednostavnu formulu za rješenje rubnog problema za bilo: koje J. 


Podijelivši interval (0, 7) na male intervale oko točaka x; širine /l, očekujemo 
da silu gustoće f možemo aproksimirati zbrojem sila intenziteta 


Xxi1+ E 


[FG)dxmf(x) 4, (4) 
A 


x 5 


koncentriranim u točkama xy, X2, ...> X,« Svaka od tih sila daje rješenje 
Fla): A 6 (8%) (5) 


a sve zajedno (prema principu superpozicije) 
EJ 
: >1(%) "4 +:G (x, x). (6) 
iz 


Ako A teži nuli, gornji izraz teži integralu 


1 
u(x)= [G(mx)f(x)dx. (7) 
0] 


10.1. Teorem 


Rješenje rubnog problema (1)—(3) dano je formulom (7). 
Dokaz. Iz (7) i (9.42) slijedi 


x ; l 
u(x) =v() | v(x)f(x)dx' +o0(2) i v&)f(x)dx, (8) 
0 x 


x l 
u (x) = (| v(x)f (x) dx +v (| w(x)f(x)dx. (9) 
0 x 
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Neka su u točki x € [0, 7] funkcije a', b i f neprekidne. Tada iz (9) do- 
bivamo 


(a(2)u' (= (a(0)w [| vaYf()d' + (aw) 
0 
ti 
JvGa)fG)dr+a(o gw) —vG)owG)f(2) (10) 


Funkcije » i w zadovoljavaju jednadžbu (9.1) pa je 
(a20')' (2) = b (2) (2), (av') (2) = b (2) v (2). (11 
Iz uvjeta (9.41) i Liouvilleove formule (v. Dodatak, $ 2) slijedi 


a((V(x)w(Q) —v(x)w(x) =a(A) W(x) =a (0) W(0) = —1. 


(12) 
Iz (10)—(12) dobivamo 
(aG)w()Y =blo)vo)[|v()f(w)dr +bG)0(x) 
0 
l tu 
[eGafG)d'—f&=V[GGx)fG)d'—f(x)= 
x 0 
— bia)u() —f(). (13) 


Prema tome, funkcija (7) je rješenje jednadžbe (1) na segmentu [0, 7). 
Dalje imamo 


i 
au (O) kato [(«F —BGQ0, #))109 d/'=0,. (4 


0 


1 
“D+2u0=[(7 SE +060) 160) 4 -4 (i 
0 


jer Greenova funkcija zadovoljava uvjete (2) i (3). 
10.2. Posljedica 


Vrijedi 
G(xx') >0 (16) 


za sve x,x' e (0,2). 


Dokaz. Prema posljedici 9.4. funkcija G ima isti znak za sve x, x' e 
€ (0, 7). Neka je funkcija u rješenje rubnog problema u slučaju f(x) = 
= 1. Tada je 


l 
u(x) = | G(x,x')dx', - (17) 
8 


To 


[«tae=[ [06029 dx da“. (18) 
0 0 0 


S druge strane, iz jednakosti 
(au —bu+1=0 (19) 
dobivamo 


f udx = — [ (au'Y u dx + f bu? dx, (20) 
to 0 0 


ili (pretpostavljajući, na primjer, BP > 0, y > 0 i uzimajući u obzir (2) 
i (3)) 


l l 
IE de = [ad Te bu doe Za (Du2(D + rad (0) > 0. (21) 
0 


0 


Iz (18) i (21) slijedi 


2 
f G (x, x") da“ dx > 0. (22) 
0 . 


nora 


Funkcija G ima isti znak kao integral na lijevoj strani u (22); prema 
tome, G je pozitivno na kvadratu (0, 7) x (0, 2). 


Formula (7) jednostavna je i omogućuje neke opće zaključke o rješenju rub- 
nog problema. 


Ako je f >0i/J/ #0, onda je 


l 
u()= [Ga x)fG)dr >0, (23) 
Z | 
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jer je Greenova funkcija pozitivna. Drugim riječima, ako linijska sila djeluje u 
jednom smjeru, taj isti smjer ima i progib ; preciznije, ako je f nenegativno 1 ragličito 
od nule, rješenje rubnog problema (1)—(3) je pozitivno.* 


Drugo važno svojstvo koje slijedi iz formule (7) jest korektnost. Imamo 


“= | ELJEDIA > dx, (24) 
0 
ve oo < [jeGe FG) dar gren jeee (02 
lu' (x)| S h max |f (2)|, (26) 
xe[0,1] 


gdje smo stavili 


1 
h = max Poe) dx. I (27) 
xe[0,1] 2x 
0 


Dakle, ako je f malo (tj. max |f| malo), onda je 1 rješenje rubnog problema malo (mala 
deformacija). 


10.3. Zadaci 


10.3.1. Dokažite da rješenje rubnog problema ovisi monotono o gustoći vanjs- 
ke sile. 


Rješenje. Neka je 
—(auY +bu=f. (28) 


—(au) +bu = (29) 
uz iste rubne uvjete. Oduzimajući (29) od (28) dobivamo 

—(a(u—u)Y +0 —uw)=f1—f. (30) 
Prema (7) iz toga slijedi 


1 
u,(x) —u(x) = / G (x) (fi (X) — f2 (x9) dx. (31) 
0 


* Taj zaključak vrijedi, naravno, uz pretpostavku (9.4) i homogene rubne uvjete. 
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Iz (31) zaključujemo: ako je f, > fx, onda je uy > u2. 


10.3.2. Neka za gustoću vanjske sile f vrijedi ili / > 0 ili f < 0. Dokažite da 
rješenje rubnog problema ovisi monotono o koeficijentu elastičnosti sredstva. 


Rješenje. Neka je 
—(au +bm=f, (32) 
—(au) +bu=f (33) 
uz iste rubne uvjete. Oduzimajući (33) od (32) dobivamo 
—(4(u — u2)Y + — bo, = 0, (34) 
ili 
—(a(u —u)Y +bil—u) =(6—b)u. (35) 
Označimo sa G, odnosno G, Greenovu funkciju jednadžbe (32) odnosno (33). 
Prema (7) iz (35) slijedi 
1 
m (€) — wE) = G,(m,2) (b, (1) — by (u (x) dr = 
0 


11 
[ [ G1(8 2) G2 (2/0) (0, (2) — ba (#0) f (2) da da", (36) 
00 
Iz (36) zaključujemo: ako je bj >b,, onda jeu, Ku zaf>0iuw >uzaf<0. 


S 11. Varijaciona jednadžba 


Ovdje ćemo opisati rubni problem (10.1)—(10.3) na novi način. Neka je u 
ravnotežno stanje, to jest 


q—-bu+f=0q=av. (1) 
Uzmemo li kakvu god funkciju v», pomnožimo (1) sa 7 i integriramo, izlazi 


l l 
[dv dx + | (—bu+ f) v dx = 0. (2) 
0 


0 
Prvi član na lijevoj strani transformirajmo parcijalnom integracijom: 


1 1 1 H 
[go dx = | (qv)' dx — | qvdx=q(D)v(2) —q(0)v (0) =] qvdx. (3) 
0 0 0 0 
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Sad iz (2) dobivamo 


1 1 
1DoQ—q(0v(0)— [q'dx+[(—bu+f)vdx—0. (4) 
0 0 


Interpretiramo li funkciju v kao smetnju ili perturbaciju ravnotežnog položaja. 
članove na lijevoj strani u (4) možemo tumačiti kao rad pojedinih sila na »putu« 
od ravnotežnog položaja u do perturbiranog položaja u -H- v: prva dva člana označuju 
rad vanjskih kontaktnih sila, a zadnja dva člana rad vanjskih linijskih sila. Nazovemo 
li član : 


2 
— [av dx (5) 


0 


radom unutrašnjih sila, jednakost (4) možemo tumačiti kao ovaj Bernoullijev 
princip: 

Ukupan rad vanjskih i unutrašnjih sila na proizvoljnoj perturbaciji ravnotežnog 
stanja jednak je nuli. 

Pretpostavimo sad da funkcija u osim jednadžbe ravnoteže zadovoljava i rub- 
ne uvjete. Budući da će se pokazati da geometrijski i prirodni uvjeti imaju bitno 
različite uloge, mi ćemo, da bismo osvijetlili oba slučaja, pretpostaviti da je, na pri- 
mjer, na lijevom kraju dan geometrijski, a na desnom prirodni uvjet: 


u(0) = 0, (6) 
yu (D+8u(D) =0, y >0. (7) 

Uvjet (7) znači da je (zbog g = au) 
a=—-5a(nuD (8) 


Kazat ćemo da je funkcija v dozvoljena ako zadovoljava iste geometrijske twjete 
kao i funkcija u; kod nas će to značiti 


v (0) = 0. (9) 
Uzimajući u obzir (6), (8) i (9), iz (4) zaključujemo da je zadovoljen sljedeći uvjet: 


1 1 
[lav + buv) dx = [redu — Za ()u(l)v(I), za svako dozvoljeno v. (10) 
0 


0 


Imamo ovaj zaključak: 


Ako je funkcija u rješenje rubnog problema (1), (6), (7), onda ona zadovoljava uvjet 
(10). 
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Pokazat ćemo da (zbog proizvoljnosti funkcije v) vrijedi i obrat tog zaključka. 
Pretpostavimo da dozvoljena funkcija u zadovoljava uvjet (10). Prvi član na lijevoj 
strani transformirajmo parcijalnom integracijom: 


H l 1 
[ arvdx = [ (ao) dx zd (avYvdx=a(Duw(dDov(02)—a(0) (0) v (0) — 
0 0 


0 
l l 
= | (awYodx=a()uw(Dv0(0)— f (au'Y' vdx. (11) 
0 0 


Iz (10) i (11) dobivamo 


€ l ! 
fo) 
a (2) (u' (1) + rk dvi) — [cao — bu +f) vdx = 0, za svako dozvoljeno v. 
% (12) 
Posebno to vrijedi za svako dozvoljeno v koje se poništava i na desnom kraju x = 
=1. Za takvo v imamo 


1 
JI (au) — bu + f) vdx = 0. (13) 
0 


Namjera nam je pokazati da se izraz u zagradi pod znakom integrala poništava. 
Pretpostavimo suprotno, tj. da je na nekom intervalu (X,, x2) taj izraz različit od 
nule, npr. pozitivan. Među dozvoljenim funkcijama koje se poništavaju i na kraju 
x = [ odaberimo funkciju zv koja je jednaka nuli svuda osim na intervalu (x,, X»), 
gdje je pozitivna. Stavljajući u (13) v = w, dobivamo 


| Can) — bu + fu da = 0. | (14) 


%1 


To je kontradiktorno jer je funkcija pod znakom integrala prema pretpostavci po- 
zitivna na cijelom intervalu integracije. Prema tome, ne vrijedi pretpostavka da 
je izraz (au) — bu + f negdje pozitivan; analogno se obara pretpostavka da je taj 
izraz negdje negativan. Funkcija u zadovoljava, dakle, jednadžbu ravnoteže. Uzima- 
jući to u obzir, iz (12) dobivamo 


a(1) (u'(1) + Su (1) v(1) = 0, za svako dozvoljeno v. (15) 
Zbog proizvoljnosti v (i pretpostavke a (1) > 0) zaključujemo 
: Ž 
u Q)+=—u0)=0. (16) 
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Prema tome, funkcija u zadovoljava i rubni uvjet (7). Imamo ovaj zaključak: 
Ako dozvoljena funkcija u zadovoljava uvjet (10), ona je rješenje rubnog problema 
(1), (6), (7). 
Uvjet (10) zove se varijaciona jednadžba* (ili varijacioni rubni problem) 
za funkciju u. Spajajući gornja dva zaključka dobivamo ovaj teorem: 


11.1. Teorem 


Rubni problem (1), (6), (7) ekvivalentan je varijacionoj jednadžbi (10), zo jest, 
vrijedi sljedeće: 

(1) Ako je funkcija u rješenje rubnog problema (1), (6), (7), onda je ona rješenje 
varijacitone jednadžbe (10). 

(11) Ako je dozvoljena funkcija u rješenje varijacione jednadžbe (10), onda je ona 
rješenje rubnog problema (1), (6), (7). 


Zbog tog teorema varijaciona jednadžba (10) često se zove varijacioni oblik 
rubnog problema (1), (6), (7). 

Između rubnog problema (1), (6), (7) i varijacione jednadžbe (10) 

postoji asimetrija u pogledu glatkoće koja se zahtijeva od rješenja u. 
Funkcija u ulazi u rubni problem sa 2. derivacijom, a u varijacionu 
jednadžbu sa 1. derivacijom. Zato je prirodno tražiti rješenje varija- 

cione jednadžbe u klasi 3£! (0, 2). Pri izvođenju zaključka (22) teorema 

11.1 mi smo, međutim, šutke pretpostavili da to rješenje pripada klasi 

KH? (0,1), što nam je omogućilo parcijalnu integraciju (11). Ta pret- 
postavka se može dokazati. Najprije precizirajmo da ćemo dozvolje- 

nom nazivati svaku funkciju v € 3£1 (0, 7) koja zadovoljava uvjet (9). 
Inkluzija u € 562 (0,7) za rješenje varijacione jednadžbe slijedi nepos- 

redno iz ovog teorema regularnosti i naših ranijih zaključaka (v. $ 3). 


11.2. Teorem 


Ako je dozvoljena funkcija u € 3€' (0,1) rješenje varijacione jednad- 
žbe (10), onda ona zadovoljava integralnu jednadžbu ravnoteže (1.28). 
Dokaz. Stavimo 


Je) =a(u(—[GG)u()—f()d'+C, (ID) 
0 
gdje je C konstanta takva da vrijedi 


l 
[36 dx = 0. (18) 
0 


* Istaknimo da je za varijacionu jednadžbu bitan zahtjev »za svako dozvoljeno ve, 
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Neka je 


w (x) = hk (x) da. (19) 
0 


Očito je 7 e 27 (0,2), F — av e (0,1), w eJE! (0, 7), a osim togaje 
zbog (18) 
(0) = v (1) = 0. (20) 
Vrijedi formula 
(F—au)o =(F—auYw+(J—au)w, (21) 
pri čemu je (zbog (20)) integral lijeve strane po intervalu (0, /) jednak 


nuli. Tako imamo 


l 1 
0= [(F—auYw +3 —au)w)dx= -| (bu — f) wdx -+ 
0 0 


1 1 
+172 dx — | aw dx, (22) 
0 0 
odnosno 
1 1 
ka dr = [ (au 20 + bu — fw) dx. (23) 
0 0 : 


Budući da je zv dozvoljena funkcija, u (10) smijemo staviti v = vw, 
što zajedno sa (20) i (23) daje 


1 


ll 32 dx = 0. (24) 


0 


To povlači f = 0, tj. 


a(x)u (x) = ji (b(x)u(x) — f(x) da + C. (25) 
0 


4 


Veoma je važno primijetiti da se u (10) više ne zahtijeva zadovoljenje prirodnog 
uvjeta (7); on je samo ostavio »traga« u zadnjem članu na lijevoj strani. Prema tome 
varijacioni oblik rubnog problema upućuje na važnu razliku među geometrijskim 
i prirodnim rubnim uvjetima. 'Ta će razlika imati veliko značenje pri konstrukciji 
približnog rješenja rubnog problema. 
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11.3. Zadatak 


Opišite varijacioni oblik problema (10.1)—(10.3) ako su 
(2) oba rubna uvjeta geometrijska, 
(22) oba rubna uvjeta prirodna. 


Rješenje. (£) Varijacioni oblik glasi: Odrediti dozvoljenu funkciju u tako da 
vrijedi 
1 l 
| (au'v' + buv) dx — | Jfvdx = 0 za svako dozvoljeno V. (26) 
0 0 


Pri tome je skup dozvoljenih funkcija dan uvjetom v (0) = v (0) = 0. 
(14) Varijacioni oblik glasi: Odrediti funkciju u tako da vrijedi 
1 u ž 
[ad Piti | ses ki fa (0) u (0) 0) +54 (u(Do(D = 0 
0 0 


za svako v. (27) 


$ 12. Funkcional energije 


Ovdje ćemo dati još jedan oblik rubnog problema* 


— (auY +bu=f (1) 
u(0) =0 ' (2) 
vu (l)+6u(F) =0, y > 0. (3) 


Pretpostavimo da funkcija zv; zadovoljava rubne uvjete (2) i (3), to jest 


w(0) =0, (4) 
yw()+65v(]) =0 (5) 
i promotrimo izraz 
1 1 
db (w) — 5 | dane Pajmada — li de (6) 
0 9) 


analogno. 
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Poslije parcijalne integracije u prvom sumandu na desnoj strani dobivamo (uzima- 
jući u obzir (4)) 


l i 
P (o) 23 (aw? + bw?) dx — saa (2) w (1) — [ro dx. (7) 
0 0 


Izračunamo li e (7) iz (5) i uvrstimo u (7), izlazi 


l 1 
b (w) = = | (aw? + bu?) dx + 25200 u (1) — [e de (8) 
0 0 


Izraz (8) ima smisla ne samo za funkcije zv koje zadovoljavaju uvjete (4) i (5), nego i 
za funkcije iz mnogo većeg skupa. Mi ćemo u daljnjem taj izraz promatrati za doz- 
voljene funkcije vw, tj. one koje zadovoljavaju uvjet (4). Svakoj dozvoljenoj funkciji 
zv pridružen je po formuli (8) broj D (w). Preslikavanje koje svakoj funkciji iz nekog 
skupa pridružuje broj zove se funkcional. Mi, dakle, promatramo funkcional D 
na skupu dozvoljenih funkcija. Članovi na desnoj strani u (8) imaju redom ove na- 
zive: energija deformacije (u »stanju« 2v), energija elastične linijske sile, energija vanjske 
kontakine elastične sile i energija vanjske linijske sile. Stoga D zovemo funkcional 
energije. 

Neka je u rješenje problema (1)—(3). Funkcional energije odlikuje se time što 
svoju najmanju vrijednost (na skupu dozvoljenih funkcija) postiže upravo za w = u. 
Preciznije o tome govori ovaj teorem, koji je zapravo princip minimuma poten- 
cijalne energije. 


12.1. Teorem 


Problem (1)—(3) ekvivalentan je minimizaciji funkcionala (8) na skupu dozvolje- 
nih funkcija, tj. vrijedi sljedeće: 

(1) Ako je neka funkcija rješenje problema (1)—(3), onda ona minimizira funkcio- 
nal P na skupu dozvoljenih funkcija (funkcional ma njoj postiže svoju najmanju vri- 
jednost). 


(1) Ako neka funkcija minimizira funkcional PD na skupu dozvoljenih funkcija, 
onda je ona rješenje problema (1H)—(3). 

Dokaz. (1) Neka je u rješenje problema (1)—(3), a w bilo koja dozvoljena fun- 
kcija. Tada je v = w — u također dozvoljena funkcija. Uvrstimo li »v=v-+uu 
(8), dobivamo 


l S“ 
BW)=DPu+v) = 9(0) + [(avv + duo —jo)da+ Ža(Du(Do()+ 
0 


i 1 
+7 [Gao + bo) as +e O=20+7 [Ga0?+oear+ 
0 0 
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Š 2 
rate (2), (9 


gdje smo uzeli u obzir da prema teoremu 11.1 funkcija u zadovoljava varijacionu 
jednadžbu (11.10) Kako su dva posljednja člana na desnoj strani nenegativna, za- 
ključujemo 


B (u) < O (w). (10) 


(££) Neka funkcija w minimizira funkcional PD na skupu dozvoljenih funkcija. 
Uzimajući u obzir zaključak (2), imamo 


D (w) = 0 (u) (11) 


(gdje je u rješenje problema (1)—(3)). Stavljajući v = w» — u i uvrštavajući v» = 
=v+ruu (11), dobivamo 


1 
| ao" + bv?) dx + Za (Q)v2(D) =0 (12) 
0 


ili, zbog nenegativnosti pojedinih sumanada 


1 Ž 
[ e" dx = [o dx =v(l) =0. (13) 
0 0 


Iz toga slijedi v = 0, tj. w = u. 


Ovdje ćemo rezimirati dosadašnje rezultate o raznim oblicima 11: formulacijama 
ravnotežnoga rubnog problema. Te formulacije su sljedeće: 


a) Integralna formulacija. Odredrti funkciju u* koja zadovoljava jednakost 
a()u'()—a(Qw(0—fo&uE)d+[f&da=0 (14) 
0 0 


za svako x € [0, 1] # rubne uvjete na krajevima, 


b) Diferencijalna formulacija. Odrediti funkciju u** koja zadovoljava jed- 
nadžbu ravnoteže i rubne uvjete na krajevima. 


c) Varijaciona formulacija. Odrediti dozvoljenu funkciju koja zadovoljava 
varijacionu jednadžbu. 


* Ovdje se pretpostavlja u € #7! (0,1). 
** Ovdje se pretpostavlja #€2#2(0,7). 
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d) Energetska formulacija. Odrediti funkciju koja minimizira funkcional 
energije na skupu dozvoljenih funkcija. 

Kao što smo vidjeli, formulacije a)—d) su ekvivalentne u ovom smislu: Ako je 
neka funkcija rješenje rubnog problema po jednoj od njih, onda je ona rješenje po svakoj. 


12.2. Zadatak 
Napišite energetsku formulaciju ravnotežnoga rubnog problema u slučaju (1) 
kad su oba rubna uvjeta geometrijska, (17) kad su oba rubna uvjeta prirodna. 


Rješenje. (7) Energetska formulacija glasi: Odrediti dozvoljenu funkciju u 
koja minimizira funkcional 


1 1 
D(w) = | (aw? + bw?) dx — [se dx (15) 
0 0 


na skupu dozvoljenih funkcija zv. Pri tome je skup dozvoljenih funkcija dan uvjetom 
wvO=v(l)=0. 
(22) Energetska formulacija glasi: Odrediti funkciju u koja minimizira funkcional 
1 H 
5 (u) = | (zo: + boa + sa ()w2 (0) _fa (0) 202 (0) — [e de. (16) 
0 o 


$ 13. Varijacioni račun 


Minimizacija funkcionala energije, koju smo razmatrali u $ 12, poseban je 
slučaj važnog problema o ekstremnim vrijednostima funkcionala. "Tim problemom bavi 
se varijacioni račun. Ovdje ćemo na jednom primjeru pružiti uvid u to područje. 


Neka je p (x, y, z) dovoljno glatka funkcija varijabli x, y, ze R. Promotrimo 
funkcional 


l 
Pio= vau) wx) da 0) 
0 


na skupu dozvoljenih funkcija w, to jest funkcija koje zadovoljavaju uvjet zv (0) = 
= 0. Pretpostavimo da na tom skupu prima ekstrem (minimalnu odnosno maksi- 
malnu vrijednost) na funkciji u. Tada za svako dozvoljeno v funkcija 


Pu+2v) (2) 


varijable A e R ima ekstrem u točki A'= 0. Kako znamo iz diferencijalnog računa, 
to povlači 


[2 VWqu+2 2) -o. (3) 


4=0 
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Izraz na lijevoj strani u (3) zove se varijacija funkcionala na funkciji u i obi- 


lježava sa 
iP (u, V), 


pa (3) možemo pisati ovako: 
ŠW(u,v) = 0 za svako dezvoljeno v. 


To je varijaciona jednadžba funkcionala . Izračunajmo & Y/ (u, v). 
Iz 


1 
VW (u-+ 4) = [v (x, u (x) + Av (x), u' (x) + 20 (2)) dx 
o 


slijedi 


1 


8" (u,o) = | ( (52) Gent zv) ) je (52) u(1), 1 (29) o" (0) dz. 


0 
Kraće pišemo 


Padaj [ze v+bo ') a 


0 


Varijaciona jednadžba, dakle, glasi 


t 
[GZ v+5 s o) a« = 0 za svako dozvoljeno v. 
o 


Provedemo li u drugom članu parcijalnu integraciju, dobivamo 


[62-629) ue :+(5% 2) GuOu (Do) =0 


Zbog proizvoljnosti v slijedi (kao i u 511) 


(52) Lu.Dea=0 
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(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Dakle, funkcija u na kojoj funkcional (1) postiže ekstrem zadovoljava diferencijalnu 
jednadžbu (11) £ prirodni rubni uvjet (12), a po pretpostavci i geometrijski uvjet 


u (0) = 0. (13) 


Tako se problem ekstrema funkcinala Y svodi na rješavanje (općenito nelinearnog) 
rubnog problema (11—(13). Jednadžba (11) zove se Eulerova jednadžba funkcio- 
nala Y. 


Najčešće funkcija w ne ovisi o x (2 = ): u tom slučaju Eulerova jednadžba 
se pojednostavnjuje i prima oblik 
Iw 
UP o o = mist. (14) 
du 
Zaista, zbog (11) imamo 
d o fa dy RKT) dp\' 2 
naj r == , r i, Pee dea 1 —_ 15 
al" žr) KOE Pe Ty “ (5%) u 


iz čega slijedi (14). 

Uz neke pretpostavke o funkciji w vrijedi i obrat gornjeg zaključka, tj. funkcio- 
nal w postiže ekstrem na rješenju problema (11)—(13). Da bismo to zaključili, 
često se zadovoljavamo fizikalnim ili geometrijskim argumentima. U nekim sluča- 
jevima rješenje je dano formulom. Mi ćemo navesti jedan takav slučaj, tzv. problem 
brahistohrone, važan za povijest matematike*. 


13.1. Primjer 


U vertikalnoj ravnini xy (sl. 1) zadana je točka (X9, 9). Među krivuljama koje 
prolaze kroz ishodište i točku (X9, Yo) treba odrediti onu po kojoj teška materijalna 
točka, ispuštena iz ishodišta, kližući se bez trenja stigne u točku (%9, Y9) za najkra- 
će vrijeme. "Takva krivulja 4 zove se bralustohrona**. Neka je jednadžba te krivulje 


Xa 


X 


Xo, Ye). 


Skka 1. 


* Problem je postavio J. Bernouli 1696, a riješio ga je sam 1. Newton iste godine. 
** Grčki brahistos znači najkraći, a kronos vrijeme. 


63 


y=y(2). (16) 


Prema zakonu o sačuvanju energije brzina materijalne točke na mjestu (x, y (x)) jed- 
naka je 


Zay (4): (17) 
"To znači 


SQ)", (8) 


gdje je s duljina luka krivulje %#, računata od ishodišta, a £ vrijeme. Iz (15) i (16) 
slijedi 


dr 1 
do ore s 
ili 
dt dx 1 
dx 6 Oo u 
Uzimajući u obzir da je ds = (1 + y'2)!/2 dx, dobivamo 
dr = 1 +7y? 1/2 
Roli) (21) 
ili 
i +9)? a 
i= X. 22 
J 2 gy ) < o 
Prema tome, problem brahistohrone svodi se na minimizaciju funkcionala 
0 
POE | p (99) de, (23) 
0 
gdje je 
A 1+9y'2\1/2 
a (24) 


Prema (14) minimizirajuća funkcija zadovoljava jednadžbu 


1+ !2\ 1/2 y? 
USS oca ie i 


ili 
MERI &. (26) 
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gdje je C konstanta. Stavimo li 


Pe Pa AA 
y =tg0, <V< 7 (27) 
dobivamo 
C C 
Iz (27) i (28) slijedi 
tgo=—Csin2o + o (29) 
ili 
SE O Eb oso) (30) 
do | ' 
Iz (30) dobivamo 
x=(C —Co— S sin2o, (31) 


gdje je C, konstanta. Jednadžbe (28) i (31) parametarski označuju krivulju #*. 
Da bismo lakše prepoznali o kojoj se krivulji radi, uvedimo novi parametar 


#=z—20 (32) 
Dobivamo 
x=D+R(&8—sn8),y=R(—cos9), (33) 


gdjesu Ri D konstante. Uzmemo li da je u ishodištu # = 0, izlazi D= 0. Konstan- 
ta R dobiva se iz uvjeta y(X449) = y9. Krivulja 


x=R(A—sin2)y=R( —cos) (34) 


zove se cikloida. Nju opisuje točka kružnice radiusa R, koja se bez klizanja kotrlja 
po x-osi (sl. 2)**. 


Slika 2. 


* Ovdje nas fizikalni argumenti uvjeravaju da je dobivena krivulja brahistohrona. 
** Način kojim smo riješili problem brahistohrone ima niz nedostataka; opravdanje svih za- 
ključaka zahtijevalo bi mnogo prostora (v. L. Švarc, Analiz, Tom I, Izdat. »Mir«, Moskva, 1972, 
str. 384—389). 
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Neka je na skupu dozvoljenih funkcija zv pored funkcionala (1) zadan i funkcio- 
nal 


l 
PW = [p(w(2), w (a) ds (35) 


0 


gdje je p (x, y, 2) dovoljno glatka funkcija varijabli x, y, z e R. Neka je V skup doz- 
voljenih funkcija zv koje zadovoljavaju uvjet 


D(w) =1. (36) 

Izoperimetrički problem sastoji se u određivanju funkcije na kojoj fun- 
Kkcional / postiže ekstrem na skupu V. 

Neka je ze V rješenje izoperimetričkog problema. Promatrajmo funkciju 

u+aV + av, (37) 


gdje su », i v, dozvoljene funkcije, a a, i a» brojevi. Neka je 


f(a>4) =P(u+ av + a0). (38) 


Imamo 
I(00=00(u) =1. (39) 


Izračunajmo derivacije funkcije f u točki (0, 0). Imamo: 


l 
(54) = (zm E 3") dx =80(4, 0). = (40) 


&i=42=0 


Pretpostavimo da funkcija u ne zadovoljava varijacionu jednadžbu funkcionala D. 
Tada postoji bar jedno v» takvo da je 


8D(u,v,) # 0. (41) 


Uz takav izbor funkcije v»» imamo 
I 40. (42) 


Prema Teoremu o implicitnoj funkciji iz (39) i (42) slijedi da je u nekoj okolini točke 
ay = 0 definirana funkcija a» = a, (a,) takva da vrijedi 


a2 (0) = 0, (43) 
f(apa,(a) =1. (44) 
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Za dovoljno malo |a,| iz (44) slijedi 
D(u+av +a(a)v)=1. 
Iz (43) i (45) zaključujemo da vrijedi 
(se; v (u + av + a2 (a1) 0) =0 
4, ei=0 


za svako dozvoljeno v,. Iz toga dobivamo 


1 1 

da, ov EN) 
[ato oil + (g), ou + ze 
0 0 


do) 8W(u0,) = 0. 
da, %=0 


ili 


Đ Fu, V,) + 


Iz (44) nalazimo 


Lo JR 
(da, ka 2 + = oK (u, Vi) 
(a m 

0 42 
Neka je 
Pe 8w (u, V2) 
= 87(0 02)" 
Iz (48) —(50) dobivamo 
čPV(uv) +2180(u,v,) =0 
ili ' 


8(P+2P) (4% v0,) = 0, 


za svako dozvoljeno v,. Imamo, dakle, ovaj zaključak: 


01) dx — 0, 


o Mk 
&i=Q2=0 


| (45) 


(46) 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


Ako je funkcija u rješenje izoperimetričkog problema i ako ona ne zadovoljava 
varijacionu jednadžbu funkcionala D, onda postoji broj A takav da ona zadovoljava 


varijacionu jednadžbu funkcionala W +26. 


Iz varijacione jednadžbe (52) slijedi Eulerova jednadžba i i prirodni rubni uvjet 


funkcionala Y/ + 2 0. 


13.2. Primjer 


Teška nerastezljiva homogena nzz visi obješena o krajeve. Treba odrediti rav- 
notežni položaj niti polazeći od principa minimuma potencijalne energije. 
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Neka su krajevi niti točke (0, 0) i (x, y1)> X, > 0 (sl. 3); neka je duljina niti 
jednaka (, a njen položaj krivulja y = y (x). Potencijalna energija niti je 


Py) =o08[y()VIFy7(G) dx. (53) 
0 


0 Xi 


bi) xa. J1) 


Slika 3. 


Duljina krivulje » = y (x) je 


&() = [VTTy'7 (0) dz. (54) 
0 


Prema principu minimuma potencijalne energije ravnotežni položaj niti je funkcija 
y (x) koja zadovoljava rubne uvjete y (0) = 0, y (x4) = y, i minimizira funkcional 
(53) uz uvjet 


by) =1. (55) 


Znamo da postoji broj A takav da funkcija y zadovoljava Eulerovu jednadžbu fun- 
kcionala 


Py) +10(9) = [ (gy + DTP? dr. (56) 
0 


Budući da podintegralna funkcija ne ovisi eksplicitno o varijabli x, prema (14) 
vrijedi 
y>-b=aVTry7, (57) 


gdje je a konstanta i gdje smo stavili A = —p gb. Pretpostavljamo netrivijalan slučaj 
y # const. Zato je a # 0. Supstitucijom 


y—b=achz (58) 
dobivamo ' 


i aa. (59) 
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gdje je a konstanta. Iz (58) i (59) slijedi 


x—a 
y—b=ach S. 


(60) 


Krivulja (60) zove se lančanica (jer opisuje i ravnotežni položaj lanca obješena o 


dvije točke). Konstante a, b i a određuju se iz rubnih uvjeta 


9(0) =0,y() =9> . 
te iz uvjeta 


[Vi FP dx=, 
0 


koji zadaje duljinu niti. Imamo ove nelinearne jednadžbe: 


b= = gdi, 
a 


YI —b=ah"—=, 


Xi —a 
a 


sh 


Iz (63) i (64) dobivamo 


Yi = —2ashushv, 
gdje je 


Iz (65) slijedi 
I=2achushv. 
Iz (66) i (68) dobivamo 


Očigledno je 


I>VAT+yi >9y0 


nie sh O, x, <. 
a a 


(61) 


(62) 


(63) 


(64) 


(65) 


(66) 


(67) 


(68) 


(69) 


(70) 


pa imamo |y,/Z <1. Prema tome, jednadžba (69) ima jedinstveno rješenje u 


(sl. 4). Iz (66) i (68) slijedi 
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Vz—y =2ashy»; (71) 


Slika 4. 


dijeleći to sa x,/2a = v, dobivamo 


shv_ Viž=ga 
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; = (12) 
Iz Mac Laurinovog reda za funkciju sh » nalazimo 

ER E Li ARE (73 

sri sits . ) 


Iz toga zaključujemo da funkcija sh »/v strogo monotono raste od 1 do oo kad» 
raste od 0 do oo. Prema tome, funkcija sh »/v prima svaku pozitivnu vrijednost, 
i to samo jednom. Iz (70) slijedi 


a >1, (14) 
1 


pa jednadžba (72) ima jedinstveno rješenje » (sl. 5). Konstante a i a određujemo sa- 
da iz (67). 


Slika 5. 


13.3. Zadatak 
Izvedite jednadžbu (60) polazeći od jednadžbe ravnoteže niti. 
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Rješenje. Neka je mi kontaktna sila. Fednadžba ravnoteže je 


— 


dq + 
—— — = E 75 
a 08=9 os 


ds je element luka (koji se u ovom slučaju ne može zamijeniti sa dx). Srla g je tan- 
gencijalna na mt: 


16) =009769=4(09(% sa St ao (16) 


Rastavljanjem u komponente i upotrebom formule 


d 1 d 
ds ui 1 i (yi da (77) 
dobivamo 
d q 
=) (78) 
dx/1+(0) 
đ dy 
= ===>) = 08VI +09. (79) 
&lyreoy) 


To je sustav od dviju jednadžba s dvjema nepoznatim funkcijama y i q. Iz (78) 


dobivamo 
q = const. \1 --(y')?. (80) 


Uvrštavajući (80) u (79), za »y dobivamo jednadžbu drugog reda s općim rješenjem 
(60). 


$ 14. Metoda konačnih elemenata 


Promotrimo sada ravnotežni rubni problem s praktičnog gledišta. Do sada 
nismo (osim u slučaju b = 0, S 5) dali konstrukciju rješenja koja bi mogla zadovo- 
ljiti praktične potrebe. Iako je (10.7) eksplicitna formula za rješenje, primjena te 
formule pretpostavlja poznavanje Greenove funkcije; s druge strane, Greenovu 
funkciju u slučaju b :# 0 možemo eksplicitno izračunati samo u posebnim slučaje- 
vima. Preostaje da se rješenje rubnog problema odredi približno uz unaprijed za- 
danu točnost. U tu svrhu zgodno je poći od varijacione formulacije problema 
($ 11). 

Da bismo konstruirali približno rješenje, razdijelimo segment [0, 7] na n +1 
jednakih dijelova čvornim točkama (čvorovima) X1> X22 -..> Xn i stavimo xo=0, X,41 = 
= 1(sl. 1). Kažemo da je funkcija v po dijelovima linearna (uz razdiobu x,, £2,...Xn) 


"1 


ako je ona na svakom intervalu (%;, Xx+1) (& = 0, 1, ..., 1) linearna i ako je u svakom 
čvoru x, (& = 1, 2, ...,n) neprekidna (sl. 1). Skup svih takvih funkcija označimo 
sa V (0, 7) (kraće V/). Funkcija zv e V! jednoznačno je određena svojim vrijednos- 
tima 2, 01, -..> 0n+1 U ČVOrOViMA X9, X1 ...> Xna+1+ Proizvoljna superpozicija fun- 
kcija iz V, je funkcija iz V)*. Očekujemo da se svaka funkcija na segmentu [0, 1] 
može dobro aproksimirati nekom funkcijom iz V! ako se uzme dovoljno gusta razdioba, 
to jest veliko n (sl. 2). 


ELO Xi X2 i Xa i=XnH 


Slika 1. : Shka 2. 


Svakom čvoru x; (k = 0, 1, ...,n + 1) pridružit ćemo po jednu posebnu fun- 
Kkciju 9x € V;, definiranu formulom 


Oe. [Xk=13 Xka1] 


1 
P% (x) = I,“ = Kuli): x € [sx—a; X] (1) 


1 
Noa — x), x€ [X Zkanals 


gdje je h = ll(n + 1) (sl. 3). Funkcije 9, Py .--> Pn+1 imaju ovo svojstvo: 


Ozajsk 


2 
Ilzaj=k. & 


PQ) = 
Svaka funkcija w e V! može se prikazati kao superpozicija funkcija P;: 


n+1 j 
2 (x) — 2 9x (x) = 29 P0(X) +, Grid) e Peba Pr (€) + ja Rece LO): 
(3) 
gdje je 
u, = 0 (X). (4) 


Zaista, zbog svojstva (2) funkcije na lijevoj i desnoj strani u (3) podudaraju se u čvo- 
rovima pa se podudaraju svuda. Prikaz (3) je ga funkciju w e Vi; jedinstven, jer zbog 
svojstva (2) iz (3) slijedi jednakost (4). Prema (3) svaka funkcija iz V; je određena 
konačnim brojem (n -H 2) parametara pa je prikladna za numeričko računanje. 


. 


* To znači da je V,! linearni prostor. Očigledno je. Vl € #f! i VI = 302. 
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Navedimo još formulu za derivaciju funkcije w e V;. Iz (3) dobivamo 


n+1 
w'() = X 9(), (3) 
pri čemu je 
0, x € Is Kk4+1) 
1 


»X€ Boas Xu) (6) 


> xE [X Kk+1)- 


Funkcija ww je na svakom intervalu [X,_1, Xx) konstanta, a u čvorovima ima skokove 
(sl. 4). 


XK-1 *k XK 
Slika 3. Slika 4. 


I ovdje ćemo razmatrati problem s karakterističnim rubnim uvjetima : 


— (uY +bu=f (7) 
u(0) = 0 (8) 
yu (D+84(0)=0y>0. (9) 


Njegova varijaciona formulacija glasi ovako: odrediti dozvoljenu funkciju u koja 
zadovoljava varijacionu jednadžbu 


1 ' 1 
[rado + buv) dx +5a (2)u(l)v(1) = [redo za svako dozvoljeno v; (10) 
0 0 


pri tome su dozvoljene funkcije one koje zadovoljavaju rubni uvjet (8). Za svaku 
dozvoljenu funkciju w e V, vrijedi 


n+1 
0 () = 2 u 9x (0) (1) 


Budući da očekujemo da se svaka dozvoljena funkcija može dobro aproksimirati 
nekom funkcijom oblika (11), približno rješenje problema (10) tražit ćemo u tom 
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obliku. Reći ćemo da je dozvoljena funkcija # e Vi približno rješenje problema (10), 
ako vrijedi 
l 8 1 “ 
| (au'v' + buv) dx 1 (u ()v(l) = [od za svako dozvoljeno v € Vi, 
0 0 


(12) 


Uzimajući u obzir da se za veliko u skup V; »malo« razlikuje od skupa svih dozvolje- 
nih funkcija, možemo reći da se od približnog rješenja zahtijeva da »približno« 
zadovoljava varijacionu jednadžbu. 


Da bi gornja definicija imala smisla, potrebno je 


(i) dokazati da približno rješenje postoji, te naći efektivan postupak za njegovo 
određivanje; x 


(1i) dokazati da s rastućim n približno rješenje po volji dobro aproksimira pravo 
(točno) rješenje i po mogućnosti dati ocjenu greške. 


Približno rješenje u i proizvoljnu dozvoljenu funkciju v € V;, pišemo u obliku 


u = DA Pro x (13) 
k=1 i 
n+1 

v—žumo (14) 


pri čemu se niz %,, 42>...> Un+1 mora odrediti, a V,, 02, ..., Vn+1 SU proizvoljni. 
ZbOg 4x = u (xx) brojevi %, su približne vrijednosti rješenja u čvorovima x;. Uvr- 
štavajući (13) i (14) u (12), dobivamo 


1 n+1 n+1 l n+1_ n+1 
[aQ)% 140 (2) 5 VP;() dx + ji bi) S uno) Sv) + 
) kZi jši 0 kZi jzi 


fo) s I n+1 
+—a(?) Un+1 Vnai = IRAC) a V;Py (2) (15) 
7 0 k=1 
ili 
n+in+l e n+1 
Žž ž (4x +bju)ikvj= X1iV» (16) 
k=lj=1 i=i 
gdje je* 
, , O 
ax = fPa()P 0) P GRISE SE a Šk nan» (17) 
0 


* 9,, je Kroneckerov simbol: &,y = 0 zar sib, = 1. 


bu= [bm )y() do (18) 
0 
= [f69q()dx (19) 
0 


Budući da je izbor niza 2,, 92, ..., Vn+1 proizvoljan, možemo u (16) staviti redom 
(V,, V2, MARE) = (1, 0, 0, +...) 0) (20) 

(0, 1,0, ..., 0) 

(0, 0,0 ..., 1). 


“To daje 
n+1 ja 
Zla + bux) 4x =/f, 
n+1l s 
zla +04 = (21) 
n++1 o 
pao + Duran) Uk = Jn+1> 
ili 


n+1 
Zan +bg)te=fpi=12..n+1 (22) 


To je linearni sustav od n ++ 1 jednadžbi za n -+ 1 nepoznanica 21, 22, ...> pere 
Stavljajući 


Aj = Gjk + ba (23) 
sustav (22) pišemo u obliku 
mH : 
žAni=fni=b2on+1 (24) 
J= 


Primijetimo da je prema (17), (18) i (23) matrica A = (Ay) tog sustava simetrična: 
Aju = A, 43 (25) 


ona se obično naziva znatricom krutosti. 


Matrica A je regularna, pa sustav (24) ima jedno i samo jedno rješenje (4, 
Za» ...> Un+1)5 ono zadovoljava jednadžbu (16), to jest daje aproksimaciju rješenja 
problema (10) u čvorovima x,, 2, ...> Xn41- 
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Dokazat ćemo regularnost matrice A. Za bilo koji izbor 7,, V2, ...> Vn+1 Prema 
(17), (18) i (23) imamo 


n+1 n+1 l j : . fo) A 
S Auxuoj= X uv jan) PA) de+—a(2)V,41 + 
kj=1 kj=1 0 7 
n+1 l l n+i B n+i A 
+ Xao[(bn)pe)dr= [ac %9469)) ( o, 9) ()) dx + 
kj=1 0 0 kzi jzi 


b5aDoba + [oc iS %96)) (2970) da = 
1 "+1 : 2 Š 1 n+1 2 

= [209 (% mr4(9) dah =o()vja [o zo) dx. (26) 
0 kzi M 0 k=1 ' 


Pretpostavimo da je 


n+1 
S Axuk=031=12..n+ 1. (27) 
k=1 


Množenjem j-te jednakosti (27) sa v;i sumiranjem poj = 1,2, ...,1n -- 1 dobivamo 
n+1 
5 Ajax VV); = 0. (28) 


hk=1 


Iz (26) i (28) slijedi (zbog a (x) > 0, b (x) > 0) 


n+1 ; n+1 , 
> un=(Xun) =0 (29) 
kzi kzi 
ili 
n+1 
S, U Px = const. (30) 
k=1 
Zbog g, (0) = 0 imamo 
n+1 
Sun) =0 (31) 
k=1 


za sve x € [0, 7]. Iz toga dobivamo (stavljajući X = X,, X2, ...> Xn+1) 
VUUZE=V=.. =V41 =0. (32) 


Prema tome, iz (27) slijedi (32), a to znači da je matrica A regularna. 


Pozabavit ćemo se ukratko i efektivnim rješavanjem sustava (24). Primije- 
nit ćemo klasičnu Gaussovu metodu eliminacije bez permutiranja vedaka i 
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stupaca*. Situacija je računski posebno povoljna zbog toga što je matrica A trodija- 
gonalna, tj. vrijedi 


Ap =0za lj— ki > 1; (33) 


to slijedi iz činjenice da produkt dviju funkcija p; i p; (kao i njihovih derivacija) 
iščezava ako one nisu susjedi (Ik —;| > 1). Točnije, imamo 


(34) 


gdje je 


1 “k+1 1 o“ 
Ak = Au =3z 1) a (x)dx + 3 J b (x) (x — Xx. 1)? dx + 


*k-1 *k_1 


1 “+1 
zle Fo J b (2) (Xx+1—2)7 dx, 


%k 


k=1,2 51; (35) 
1! : fo) 1“ 2 
An41 E Anti n+i “miso ka a (X, — x)? dx, (36) 


lo*k+1 1 zk+1 
B, = Aa+i = — 3 j a(x)dx +52 1) b (2) (kar — 2) (X — Xx) de, 
ETA XR 


k=1,2,..,m. (37) 
Za sustav (24) moguće je Gaussove eliminacije izvesti na ovaj jednostavan način**: 


direktni hod: 
k=1,2,...,n: 
s=BdlA (38) 
Ani = Akti — SB 
Heri =Hra1 — Sf 


n+1 
* To je moguće zbog pogitivne definitnosti matrice A; naime, iz (26) slijedi 2 Ajxvxv;>0, 
jk=1 
a kao što smo vidjeli, (27) povlači (32). 


** U formulama koje slijede upotrijebili smo oznake uobičajene u opisivanju računskih pos- 
tupaka (algoritama), podešenih za računanje na stroju. Tako 4,41 = 441 — Sd, Znači: »izračunaj 
đk+1 — 5x i dobiveni rezultat upiši na mjesto gdje je stajao a, rt. 


TT 


obratni hod: 


k=nn—l1...1l: 


= —B (39) 
pz =/dA,. 


U formulama (19), (35)—(37) potrebno je izračunati zadane integrale, što za 
općeniti izbor podintegralnih funkcija zahtijeva upotrebu neke metode približne 
integracije*. Jednostavna mogućnost je zamijeniti funkcije a, b i f njihovim linear- 
nim interpolacijama na svakom segmentu [%;, Xx +1]. Tako umjesto f uzimamo 


N n+1 
HOI? (40) 
što umjesto f, daje 
s. “o ; hn fGax-D)+4VG0 +IlaHi) k=12..N 
he — 200 | U Peds — 6 \G() + 2f (a1) k=n +1, m 


pri čemu smo se koristili formulom 


4k=j=12..n 


PRP 72 k=j=0n+1 za 
0, (&—j] > 1. 


Na isti način postupamo i s funkcijama a ib. 
Izračunavanje približnog rješenja problema (10) koje smo ovdje opisali zove 
se metoda konačnih elemenata**. Funkcije px zovu se linearni elementi. 
Pogledajmo u kakvoj je vezi približno rješenje 2 s funkcionalom energije Q 
(8 12). Za proizvoljno v = vi; (i +V2P2... + Vn+1 Pn+1 izračunajmo vrijednost 
funkcionala D na funkciji % + v, uzimajući u obzir (16) te formule (17)—(19) i 
(23): 


BG+V)=0(8) + [ aa'v' + bižo fox +>a(DEDoD+ 
0 


n+1 s 
ra Age Uk Vi — 


j.k=1 


l < "2 2 0) 2 E EA 
poi: a de e () = 0 (4) + 


* Tu metodu biramo tako da pogreška koja pri tom nastaje bude istog »reda veličine« kao 
i onakoju činimo opisanom diskretizacijom (v. (55)). Točnija približna integracija samo bi povećala 
opseg računa, a da se ukupna pogreška ne bi bitno smanjila. 
** Riječ »konačan« znači da se umjesto »beskonačno malih« (diferencijalnih, infinitezimalnih) 
promatraju makroskopski komadi tijela (u našem slučaju segmenti [x,_ 1, *x]), kojima se unaprijed 
propisuje oblik pomaka (u našem slučaju linearni). 
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l 
Shah [ (002 +bod)dx + La(Deža = 0B) + 
0 y 
/ ' o) E 2 
pe o pov a (43) 


Odatle odmah slijedi 
PR) <a +, ' (44) 


pri čemu znak jednakosti vrijedi onda i samo onda ako je 
1 ' ' ' 
ji (av? + bo?) dx = 0, #441 =0, (45) 
0 


tj. v = 0. Prema tome, približno rješenje %# problema (10) minimizira funkciona 
energije D na skupu svih dozvoljenih funkcija iz Vi; i jedino je s tim svojstvom 
Uzimajući u obzir da točno rješenje minimizira funkcional energije na skupu svih 
dozvoljenih funkcija, možemo reći da gornji zaključak i na drugi način opravdava 
naziv »približno rješenje« za 7. 


Preostaje nam da preciznije istražimo vezu približnog rješenja x 
s točnim rješenjem u*. 


U skladu sa (40) linearnu interpolaciju proizvoljne funkcije v na 
segmentu [0, 7], uz razdiobu x,, x2, ..., X,» označavat ćemo sa 8. Za 
dozvoljenu funkciju v imamo 


TOB>OINO (46) 


14.1. Lema 


Za ve#? (0,1) € svako x € [0,1] vrijedi 


lo (x) — 2 (x)| < hmax |v" (&)|, (47) 
£e[0,1] 
h2 
lo (x) — 2 (I < g max (v" (£)|. (48) 
[0,1] 


Dokaz. Zbog v (x) — 2 (4x) = 0, V G41) — D (Xu+1) = 0 postoji 
točka &, € (Xx %4+ 1) takva da je 


v(&)—v(6) =0. (49) 


* Približno rješenje ovisi i o broju čvorova 2, tj. imamo posla s nizom približnih rješenja, 
.. . . m o . ve. .. .. . . 
pa bi bilo ispravnije stavljati # = u"; to ipak ne činimo radi jednostavnijeg pisanja. 
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Za x € (Xx, Xx41) imamo 
v' (a) — 860) = [(0"(&) —8" (8) df — [v"(d& > (50 
Ek Ek 


jer je 8" = 0. Iz toga slijedi 


lo" (2) — 8" (&)] < lx — &] maxje" (| < hmax |o" (&)|. (51) 
šel6k,x] £e[0,1] 


Desna strana te nejednakosti ne ovisi o & pa zato ona vrijedi za svako 
x € [0, /]; time dobivamo (47). Na segmentu [%, % +1] funkcija |v (x) — 
— 2 (x)| ima maksimum u nekoj točki €, sa svojstvom (49). Neka je, 


na primjer, 6, — %x < Si Prema Taylorovoj formuli imamo 
A A 1 Zara LSTNE: 
vx) —Vva)=0=Vv(6)—?2(8) + se Gi > S "ii 


_&)=v()—0E)+zo'(Š)G— E) = (52) 


gdje je x, < & < ć,. Iz toga dobivamo 


2 


TO E ZRELI BEUKOJ 


2 
(53) 
Za x € [X Xx41] imamo 
A h2 
Jo (2) — 2 (0) < Jo) — 8 (EVI < zomax lo" (GL (69 
x€eL0,, 


Desna strana te nejednakosti ne ovisi o k pa ona vrijedi za svako x € 
€ [0, 7); time dobivamo (48). 


14.2. Teorem 
Za svako x € [0,7] vrijedi 


lu (x) — u (x)| < Chmax |[u" (x)|, (55) 
x€[3,1] 


gdje je C > 0 konstanta neovisna 0 žh. 


Dokaz. Stavljajući u formuli (12.9) w = x, dobivamo 


I 
[Gal — 29% + PGu— BA + za DEO TOY - 
0 


— GR) — (u). (56) 
Budući da # minimizira D na skupu Vi! i budući da je # € Vi, vrijedi 
Du) — Du) <D() — 0 (u). (57) 


S druge strane, stavljajući u formuli (12.9) w» = 4 (i uzimajući u obzir 
da je 4 (1) = u (1)), dobivamo 


H 
PR)— o) = S [(a(0—4)2+b(u—2))dx. (58) 
0 
Uzimajući u obzir lemu 14.1 i (58), imamo 


h* 
DZ) — Du) < aa h? max |u"' (2)[2 + bo 2; max [u (x)l2), (59) 
2 xe[0,7] 64 er0,11 
gdje smo stavili 


4o =maxa(x), bo = maxb (2). (60) 
xe10,1] x€[0,1] 


Iz (56), (57) i (59) slijedi 
1 
[at —79+b(u— d+ Za DuO-ZO? < 


0 
bah? 
< (a Tb “are max |u" (x)|? < C, h2 max [u" (x)|*, (61) 
xe[0,1] xe[0,/] 


gdje je, na primjer, C, = 1 (ao -+ bol12/64). Iz toga dobivamo 


I 


il a(u' — 8")? dx < C,h? max (u" (x))2. (62) 
D xe[0,/] 
Stavljajući 
a, = Mina (x), (63) 
x€[0,1] Ž 
imamo 


1 1 
afu—q)?dx< fat —2)dx< Ch? max lu" (x)|, (64) 
0 0 x€[0,l] 


ili 
1 


C,h? 
[u —w)?d< —a. max lu" (|. (65) 
0 1 oxel0,l] 


6 Jednadžbe matematičke fizike 81 


1z jednakosti 
u(x) — u (a) = [ (w'(£) — 2 (Bač, (66) 
0 


koristeći se Cauchyjevom nejednakosti, dobivamo 


x x l 
(u(2) — 2 (0) < [a6 Put (£) — 2" (8) dE < 1 [lu (x) — 2" (29)? da. 
0 0 0 


(67) 
Iz (65) i (67) slijedi 


IC, 


uši 


(u(x) — (x)? < W2 max li '(2)l2. (68) 


Stavljajući MZ C,/a,= C, dobivamo (55). 


Iz teorema 14.2 slijedi da približno rješenje teži zočnom rješenju 
uniformno na [0,1] kad h >0 (tj. n > oo). Štoviše, nejednakost (55) 
omogućuje određivanje približnog rješenja s unaprijed zadanom točnosti. 


Metoda konačnih elemenata je poseban slučaj R:tz-Galerkinove 
metode, u kojoj se približno rješenje traži u obliku superpozicije »pro- 
izvoljnih« linearno nezavisnih funkcija Y1, X2>...> Xn. Odgovarajuća 
matrica A je tada općenito puna, što je numerički nepovoljno. Linear- 
ne elemente kao funkcije yx predložio je R. Courant 1943. Ta metoda 
je ponovno otkrivena 60-tih godina u okviru strukturalne mehanike* i 
dobila današnji naziv. Umjesto linearnih mogu se (a kod jednadžbi vi- 
šeg reda i moraju, v. $ 17) uzimati kvadratični, kubični itd. elementi. 
Time se povećava točnost približnog rješenja, ali matrica A postaje 
punija. Metoda konačnih elemenata za dvodimenzionalne i trodimen- 
zionalne probleme danas je osnovna numerička metoda u tehnici. 


14.3. Zadaci 


14.3.1. Opišite metodu konačnih elemenata za slučaj rubnih uvjeta u (0) = 


=u(l) =0. 
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Rješenje. Približno rješenje # dano je formulom 


u) = Sa, Py (2), (69) 
k=1 


* Strukturalna snehanika se bavi problemima ravnoteže i gibanja štapova, ploča i ljusaka. 


gdje je (24 22, ...> 4) rješenje linearnog sustava 
S Apxtuk =] je= 1/24 f (70) 
k=i 

i gdje su Aj: i f; dani formulama (17), (18), (23) i (19)*. 


14.3.2. Opišite metodu konačnih elemenata za slučaj rubnih uvjeta u' (0) = 
="(0)=0. 


Rješenje. Približno rješenje u dano je formulom 
s< n+1 pa 5 : 
u (X) = S 2:0 (2), (71) 
k=i 
gdje je (Zo #1> ...» 2n+1) rješenje linearnog sustava 
n+1 . 
S Axl = 19 1 WLaNPd (72) 
k=0 


i gdje je Ajoj, k > 1 dano formulama (17), (18)i(23)za 8 =0 i 


1 # 1% ; 
am kae jd a best dx, (73) 


h hk 
Aor = Ao = — ga [a)dx +33 [ 69) x(H —x)dz, (74) 
0 1 


dok je f; dano formulom (19). Ako je b #£ 0, matrica A = (Ay,) je pozitivno definitna 
(dokaz je isti kao u slučaju (8), (9)), pa dakle i regularna. Situacija je bitno drugačija 
u slučaju b = 0. Tada kao i u (26) vrijedi 


n+1 1 n+1 : ž 
ka Au VV, = f qa (x) ( kJ Vry Pr %) dx E 0, (75) 
0 k=0 


k,j=0 


pri čemu lijeva strana iščezava ako i samo ako je 
n+1 E 
> %xp4(2) =0 (76) 
k=0 


za sve x € [0, 7]; to je ispunjeno ako i samo ako je : 


V=V=... = V,+1. (77) 


* Matrica A = (Ay) je pozitivno defimitna i ocjena (55) vrijedi i u ovom slučaju. 
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Prema tome, matrica A je singularna. Zbog toga sustav (72) nije rješiv za bilo koju 
desnu stranu (fo> fu» -->fn+1). No, prema $ 5. i točno rješenje ne postoji zab = 0 
i za proizvoljno f, već f mora zadovoljavati uvjet (5. 12). Imamo 


S4=% [ronda l f(2) ba > piooja« (18) 
j=0 i=00 i= 
pa zbog 
n+1 
S ma) =1 (19) 
izo 
i (5.12) dobivamo 
Sj (80) 
jeo 


UZ taj uvjet sustav (72) je rješiv i rješenje mu je određeno do na aditivnu konstantu, 
istu u svakoj komponenti. Interpretacija ovakvog rješenja je ista kao ona u $ 5. 
Jedno rješenje sustava (72) može se naći tako da se uzme, na primjer, uo = 0, čime 
taj sustav prelazi u sustav dimenzije n -+ 1, s matricom koja se dobiva iz A brisanjem 
prvog retka i prvog stupca. Nova matrica je pozitivno definitna. 


14.3.3. Problem 
u =—1 (81) 


uO)=2:() =0 (82) 
riješite metodom konačnih elemenata za (1) n = 3, (ii) n = 6 i rezultate usporedite 
s točnim rješenjem. 

Rješenje. Točno rješenje je 


u(x) = st — 1) (83) 


(sl. 5). Približno rješenje dano je formulom (69). 


8 —4_ 0 IH 
A=|_4 8 —4b/=ZT1I1 (84) 
O.=4. # BI 


Rješenje sustava (70) je 


(2) Dobivamo 


(85) 
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0 1/2 1 


Slika 5. Slika 6. 


(sl. 6). Približno i točno rješenje se podudaraju u čvorovima. 


(22) Dobivamo 


M4. =. 0-0 BP 1 
= M4 7 BOO. 1 
je 0 —7 14 —7 0 0 pih 1 
O 04. =TFVU M4 ET 0 T|1 

O Be < 7 447 1 

O 100 0 =T 14 1 


- - - < 


Rješenje sustava (70) je 


0 VI 27.347 SA 6/7 1 


Slika 7. 


S 15. Ravnoteža štapa 


(86) 


(87) 


Štap je tijelo kojemu je jedna dimenzija — duljina— znatno veća od drugih. 
U nedeformiranom stanju štap opisujemo segmentom [0, 7] na osi x, koji se ovdje 
zove centralna limja, te poprečnim presjekom D (x) koji se dobiva tako da se štap pre- 
siječe ravninom okomitom na centralnu liniju na mjestu x e [0, 7]. Ako je poprečni 
presjek svuda isti, štap je cilindričan (sl. 1, 2). Pretpostavljat ćemo da u nedeformira- 
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x za Dix)=const. 


Slika 1. 


y yk 


x Z A Dix)=const. 


Slika 2. 


nom stanju štapa centralna linija prolazi kroz težište poprečnih presjeka. "Vreba napo- 
menuti da ta pretpostavka ne znači samo poseban izbor centralne linije, već i neko 
ograničenje na oblik štapa. Dva važna oblika za koja je takav izbor moguć jesu ci- 
lindrični štap i štap rotacionog oblika (sl. 3). 


\ 


x z (L\ D(x) # const. 


Slika 3. 


Ograničit ćemo se na promatranje male deformacije štapa, do koje prema iskus- 
tvu dolazi pri dovoljno slabome vanjskom djelovanju. 

Eksperimentalna opažanja kao i usporedbe s realističnijom teorijom elastič- 
nosti pokazuju da se pri proučavanju male deformacije štapa može poći od sljede- 
ćih pretpostavki, 


15.1. Pretpostavka 
Svaki poprečni presjek ostaje pri pomaku nedeformiran, tj. kruto se pomiče. 


15.2. Pretpostavka 


Svaki poprečni presjek poslije pomaka ostaje okomit na (pomaknutu) centralnu 
liniju. 
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15.3. Pretpostavka 


Pri deformaciji ne dolazi do torzije*. 


Sada ćemo dati strogu formulaciju gornjih pretpostavka. Označimo sa u (x) = 


=1u,(x) +ju, (x) + ku, (x) pomak točke x centralne linije. Pretpostavka da je 
deformacija mala izražava se ovim uvjetima**: 


i ()I<I (2) 


i l 
i"đ«> (2) 
za svako x € [0, /]. Iz (1) kao i kod žice (8 I) slijedi 


li) — KOI < i, (3) 
j 


Po pretpostavci 15.1 pomak točke (x, y, 2) € D (x) jednak je 
u(x) +0 (Xx (Uy + &2), (4) 
pri čemu o (x) zovemo vektor rotacije presjeka D (x). Pretpostavka 15.2 znači 


0; =4),0, = —u,, (5) 
a pretpostavka 15.3 daje 
VW, = 0. (6) 


Tako je položaj deformiranog štapa potpuno opisan pomakom centralne linije, 
to jest funkcijom 2. 


Slika 4. 


* Ova pretpostavka nije u proturječju s činjenicom da se štap može tordirati oko osi x, kako 
je to opisano u zadatku 1.2.2. Rješenje toga prividnog paradoksa jest činjenica da je za torziju štapa 
potrebno kudikamo jače vanjsko dijelovanje negoli za poprečnu ili longitudinalnu deformaciju, 
koje ovdje promatramo. 


5 
** |x| označuje duljinu vektora u. 
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Kao i kod žice, označimo sa P(x) pomaknuti položaj točke x e [0,7] (sl. 4). 
Uvedimo sljedeće oznake: 


q (x) — kontaktna sila na presjeku D (x), 
1(x) — gustoća* vanjske linijske sile, 

m (x) — kontaktni spin** na presjeku D (x), 
5 (x) — gustoća vanjskog limjskog spina. 


Kontaktni moment na presjeku D (2x) je 
m (2) + R(2) X q(2), (7) 
gdje je R (x) radijus-vektor točke P (x): 
RG)=x+u0). (8) 
Analogno, gustoća vanjskog linijskog momenta je 


s(2) +R (2) x f(0). (9) 


Princip ravnoteže glasi: Ako je tijelo u ravnoteži, onda su ukupna sila i ukup- 
ni moment koji djeluju na bilo koji komad tijela jednaki nuli. Prema tome, za štap 
u ravnoteži vrijedi 


20) —a(0+ [fd —0, (10) 
0 


mo)+Ra) xd) —mO+[GO+RĐOxfE)dE=0 (M) 
0 


pri čemu smo zbog (1) diferencijal duljine luka na deformiranoj centralnoj liniji 
zamijenili diferencijalom luka na osi x. Deriviranjem po x € (0, 7) izlazi 


T+f=0, | (12) 
> —> > > —> —> > — 
m+RxXq+Rxq+s+Rxf=0. (13) 

Zbog (1) stavljamo 
RK xq=ixq (14) 


* Gustoće su uzete na jedinicu duljine nedeformiranog štapa. 
** Engleski to spin znači vrtjeti. 
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tako da iz (13) (uzimajući u obzir (12)) dobivamo 
-> > > > 
m +ixq+1s=0. (15) 


Jednadžbe (12) i (15) glase u komponentama* ovako: 


ax +f.=0 (16) 
q+f,=0 (17) 
a+f.=0 (18) 
m—a+s=0 (19) 
mp +q,+s, = 0. (20) 


To su zakoni ravnoteže za štap. Sad ćemo opisati zakone ponašanja štapa. Primije- 
timo odmah da će ponašanje za g, i q, biti zadano ponašanjem za 21, odnosno 1, 
kako slijedi iz (19) i (20). 

Za q, pretpostavljamo da se ponaša po Hookeovom zakonu (v. zadatak 1.2.1). 
To ovdje znači da je površinska gustoća uzdužne kontaktne sile (sila po jedinici 
površine) na pomaknutom presjeku D (x) proporcionalna derivaciji (po x) uzdužne 
komponente pomaka (4), to jest jednaka 


Eu) +02 —0o(Wy) =EU()—u()Z—w(x)9) (2) 


gdje je E > 0 Youngov modul, koji je za homogen štap konstanta, a inče je funkcija 
od x. Kontaktna sila g, na D (x) dobiva se integracijom izraza (21) po presjeku (zbog 
male deformacije svi integrali uzimaju se po nedeformiranom štapu): 


a0) = [[ Eu) —u(z—w(y)dyda= 
Di) 


—EG)SAUL)—E(u"(x) il | zdyde —E(x)u" (x) / il dy dz, (22) 
DCx) D(x) 


gdje je S (x) površina presjeka D (x). Pretpostavljamo da je S (x) > 0 za svako 
x e [0, 7]. Pošto smo pretpostavili da se težišta poprečnih presjeka nalaze na cen- 
tralnoj liniji, integrali na desnoj strani u (22) iščezavaju pa imamo 
dq, = ESu,. (23) 
=> 
Spin m (x) je ukupni moment sile (21), što daje 


m(x) = [[ (9 + zh) x Bi(uz (0) — u2(2) 2 — 15 (1) 9) dy de — 
D(x) 


= E(—u2(1,(0) —15(0) De GI + E2091), (0) +18 (0) 7, (YE, (24) 


* x-komponentu jednadžbe (15) ne pišemo jer se odnosi na torziju, koja se po pretpostavci 
15.3 zanemaruje. 
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gdje je 


1,() = [[ 97 dy dz, (25) 
D(x) 
1. (2) = [[ 2 dy de, (26) 
D(x) 
I, () = [[ ya dy dz. (27) 
D(x) 


U daljnjem ćemo pretpostavljati da se koordinatne osi y i z mogu odabrati tako da je 
I, (&) =0 (28) 
za sve x e [0,7]. Taj uvjet je ispunjen za cilindrični i rotacioni štap (v. zadatak 
15.1.1). Uz (28) iz (24) izlazi 

m, = El, u, (29) 
m,= —ELu. (30) 

Uvrštavajući (29) i (30) u (20) odnosno (19), dobivamo 
Q=—(ELu)—s, 1) 
&=—(ELu)y+s, (32) 


Jednadžbe (23), (29)—(32) su zakoni ponašanja za štap. Iz (16)—(18) i zakona (23), 
(31) i (32) izlazi 


(ESuy' +, =0, (33) 
—_ELuY=-atj,=9 (34) 
—(ELu)"+s+f.=0. (35) 


To su jednadžbe ravnoteže štapa. Jednadžba (33) je neovisna o drugim dvjema i 

opisuje uzdužnu deformaciju koju smo razmatrali ranije (zadatak 1.2.1). U daljnjem 

promatramo samo jednadžbe (34) i (35), koje opisuju poprečnu deformaciju (progib 
> => 

štapa je 4,j -+ u,k). One su također neovisne (to je posljedica pretpostavke (28)) i 

istog su tipa. Zato možemo bez gubitka općenitosti pretpostavljati da se štap progi- 

ba, na primjer, u ravnini xy. Uvodeći kraće oznake, 


u=u,q=qam=mf=f,s=sA=El,, (36) 
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odgovarajuće jednadžbe (31), (29) i (34) pišemo u obliku 


q=—(AuY—s (37) 
= Ai (38) 
(Ajtjges (39) 


Ako se štap nalazi u elastičnom sredstvu s koeficijentom elastičnosti b (x) > 0, 
umjesto (39) imamo 


(Au")"+tbu=f—s'. (40) 


"To je obična linearna diferencijalna jednadžba četvrtog reda. Funkcija koja zadovo- 
ljava jednadžbu (40) zove se ravnotešni pregib. 


U daljnjem pretpostavljamo da je s = 0, jer je to rEdoi slučaj (za izuzetak 
v. zadatak 20.3.2). 


15.4. Zadaci 


15.4.1. Dokažite da je u slučaju cilindričnog i rotacionog štapa moguće oda- 
brati osi y i z tako da vrijedi (28). i 


Rješenje. Ako je štap rotacionog oblika, onda je D (x) krug pa je 7,,(x) = 0 
za svaki izbor koordinatnih osi y, g i za svako x € [0,7]. Ako je štap cilindričan, 
onda je D (x) isto područje za svako x pa je 1,, = const. Prijelaz na bilo koji drugi 
sustav (y', 2) dan je transformacijom 


y =yc0sp — Zsin p, (41) 
zZ =ysinp + zcosp. (42) 
U novom sustavu imamo 
1): = NI Y z dy ds, (43) 
E 


gdje je D" transformirano područje D. Uzimajući u obzir da je Jacobijan transfor- 
macije (41), (42) jednak 1, dobivamo 


ip i | ( yz cos 29 + 2 in 2) dy de. (44) 
D 
Izborom 
21), za om 


izlazi 1,» = 0; taj izbor je očigledno neovisan o x. 
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15.4.2. Izvedite jednadžbu poprečne ravnoteže štapa uz pretpostavku da na 
njega djeluje velka longitudinalna sila. 


Rješenje. Ako na štap djeluje velika vanjska longitudinalna sila, onda je 
qx = a veliko (a dobiva se rješavanjem jednadžbe (33) uz zadane longitudinalne 


=> 
rubne uvjete). U tom slučaju aproksimacija (14) nije opravdana, jer q ima veliku 
uzdužnu komponentu. Pretpostavivši da je 


—> —> > —> > > —> > > 7 
R'xq=(1+ uy) X (2 +3 agh) dyR — Qe — au,k, (46) 
iz (13) dobivamo (uz oznake (36)) 


a=— m E mi= (Au'")' + au. (47) 
Umjesto (39) imamo 
(Au")" — (uy =f. (48) 


$ 16. Rubni problemi za štap 


Kod štapa imamo dvije geometrijske veličine (progtb i rotaciju) i dvije dinamičke 
veličine (stu i moment), pa će izbor rubnih uvjeta biti bogatiji nego kod žice. Na 
svakom kraju možemo zadati ili obje geometrijske veličine, ili obje dinamičke ve- 
ličine, ili jednu geometrijsku veličinu (progib odnosno rotaciju) i jednu dinamičku 
veličinu (moment odnosno silu). Uvjet u kojem se zadaje progib odnosno rotacija 
zove se geometrijski, kinematički ili Dirichletov, a uvjet u kojem se zadaje moment 
odnosno sila zove se prirodni, dinamički ili Neumannov. Ako je zadana vrijednost 
geometrijske odnosno dinamičke veličine jednaka nuli, uvjet je komogen, inače je 
nehomogen. U daljnjem promatrano samo homogene uvjete, jer i ovdje vrijedi 
princip superpozicije pa se rubni uvjeti mogu hkomogenizirati. 

Ako je u (0) = 0 odnosno o (0) = 0, kažemo da je kraj x = 0 učvršćen na pro- 
gib odnosno učvršćen na rotaciju; ako je q (0) = 0 odnosno m (0) = 0, kažemo da je 
kraj x = 0 slobodan na progib odnosno slobodan na rotaciju. 

U primjerima koji slijede opisani su najinteresantniji slučajevi rubnih uvjeta ; 
uzimamo kraj x = 0. 


\' 


16.1. Primjer 


Kraj je ukliješten (to jest učvršćen na progib i na ro- 
taciju, sl. 1): 


u (0) = 0, o (0) = ' (0) = 0, (1) 


Slika 1. 
16.2. Primjer 


Kraj je slobodan (tj. slobodan i na progib i na rotaciju): +: 


d(0) = — (Au) (0) =0, m(0) = (4u") (0) =0. (2) oke 
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16.3. Primjer 
Kraj je učvršćen šarkom (tj. učvršćen na progib i slobodan na rotaciju sl. 2): 
u (0) = 0, m (0) = (Au'") (0) = 0. (3) 


Često se efekt šarke postiže (kod jednostrano opterećenih štapova) slobodnim oslon- 
cem (sl. 3, v. zadatak 16.10.5). 


Shka 3. 


16.4. Primjer 


Kraj je užlijebljen (to jest učvršćen na rotaciju i slo- 
bodan na progib sl. 4): 


o(0) ="(0)=0,3(0)=—(AW"Y (0) =0 (4) 


Slika 4. 
16.5. Primjer 
Kraj je elastično vezan na progib i slobodan na rotaciju: 
—q(0) = —x4(0), m(0) =0 (x, >0) (5) 
ili ekvivalentno 
(Au (0) + 2, u(0) =0, 2" (0) = 0. (6) 


Slika 5. 


Takav uvjet može se ostvariti vezanjem kraja na elastično pero (sl. 5). 


16.6. Primjer 


Kraj je učvršćen na progib, 
u(0) =0, (7) 
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1 elastično vezan na rotaciju: 
—m (0) = —x2 o (0) (8) 
ili, ekvivalentno 


(Au"') (0) — x, u' (0) = 0. (9) 


"Takav uvjet može se ostvariti šarkom i kružnim perom (sl. 6). 


7 Skka 6. 


Rubni problem za štap sastoji se u određivanju funkcije u koja zadovoljava jed- 
nadžbu ravnoteže (15.40) i rubne uvjete na krajevima. 


Budući da je jednadžba (15.40) četvrtog reda, očekujemo veći 
stupanj glatkoće rješenja nego kod žice. Pretpostavljat ćemo 


Ae: (0, 1),b € 3 (0, 1,f €" (0, 0). (10) 


Pod rješenjem jednadžbe (15.40) na segmentu [0, 7] razumijevat ćemo 
funkciju u € %* (0, 7), koja u svakoj točki x € [0, 7] u kojoj ima 4. de- 
rivaciju zadovoljava jednakost (15.40). Svako rješenje jednadžbe (15.40) 
na segmentu [0, 1] ima 4. derivaciju u svakoj točki u kojoj su funkcije 
A", b i f neprekidne (v. Dodatak, $ 2). Precizna postavka rubnog pro- 
blema glasi ovako: Odrediti funkciju u e * (0, 1) koja je rješenje jed- 
nadžbe (15.40) na segmentu [0, 1] € koja zadovoljava rubne uvjete. 


Razmotrimo pitanje jedinstvenosti rješenja rubnog problema za štap. Zbog 
principa superpozicije dovoljno je analizirati Homogeni problem, tj. homogenu 
jednadžbu 


(Au')"+bu=0 (11) 
uz homogene rubne uvjete. Ako taj problem ima samo trivijalno rješenje u = 0, 


nehomogeni problem ima najviše jedno rješenje; u protivnom za nehomogeni pro- 
blem jedinstvenost ne vrijedi. 


Pomnožimo jednadžbu (11) sa » i integrirajamo od 0 do /: 


l l 
i (Au")" udx + | bu? dx = 0. (12) 
0 o 
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Prvi član na lijevoj strani transformirajmo parcijalnom integracijom : 
1 H 


| (Au")"udx = | Au"? dx + (Au) (lu (2) — (Au")' (0) u (0) — 
0 0 


— (Au) (1) u (1) + (Au) (0) u (0). 
Iz (12) i (13) izlazi 
1 


[ (Au'"2 + bu2) dx + (Au") (1) u(D) — (Au) (0) u (0) — 


0 


, — (Au) (Du (1) + (Au) (0) u (0) = 0. 


Sad treba uvažiti rubne uvjete. Razmotrit ćemo tri karakteristična slučaja. 


16.7. Primjer 
Kraj x = 0 ukliješten, a kraj x = I slobodan: 
u(0) = u' (0) u" (2) = (Au"Y (1) =0. 
Iz (14) i (15) slijedi 
j (Au"? + bu?) dx = 0. 


0 


Zbog A > 0ib > 0 imamo 


u = Cix + C,, 


gdje su C, i C, konstante; iz (15) dobivamo C, = C, =0,u = 0. 


16.8. Primjer 
Na oba kraja šarke: 
u(0) = u" (0) Z u()=u"(2D) =0. 


Iz (14) i (19) dobivamo opet u = 0. 


16.9. Primjer 


Oba kraja slobodna na progib i elastično vezana na rotaciju: 


(13) 


(14) 


(15) 


(16) 


(17) 
(18) 


(19) 
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(Au") (0) = (Au (D) =0, (20) 
(Au) (0) — zv (0) = 0, (Au) (1) +22u' (1) =0, x, >0. (21) 
Iz (14, (20) i (21) dobivamo 


I 
| (Au"? e bu2) dx +242 (0) + 2242 (1) =0. (22) 
0 


Iz toga zaključujemo da vrijedi (18) i 
buž =0, (23) 
(0) =u(1) =0. (24) 


Ako je b #0, onda zbog (23) u iščezava na nekom intervalu pa je C, = C, =0, 
tj. u = 0. Ako je b = 0, onda iz (24) slijedi C, = 0, to jest u = C2; u tom slučaju 
rješenje je, dakle, određeno do na kruti pomak. 


Lako se provjerava ovo pravilo: 


U slučaju b s O rubni problem ima najviše jedno rješenje. U slučaju b = 0 rubni 
problem ima najviše jedno rješenje ako 1 samo ako su barem dva rubna uvjeta geometrijs- 
ka, od kojih barem jedan sadrži progib. 


Kao i kod žice, u slučaju nejedinstvenosti rješenja postoje nužni uvjeti za pos- 
zojanje rješenja. Integrirajući jednadžbu 
(Au")" =/f (25) 
d 0 do (, dobivamo 


l 


(Au) (1) — (Au) (0) = [1 ds. (26) 


0 
Ako je, na primjer, q(0) = —(Au"Y (0) =0, q(D) = —(Au"Y (D =0, iz (26) slijedi 


4 


| fdx=0. (27) 


0 


Pomnožimo li (25) sa x i integriramo od 0 do /, dobivamo 


I 
1(Au"Y () — (Au) (D) + (Au) (0) = lEj dx. (28) 


0 
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Ako je na primjer, m(0) = (A4) (O) =0 mD=(AW")(D) =0, q(D = 
= — (Au"Y (D = 0, iz (28) slijedi 
4 


i xf dx = 0. (29) 


0 


16.10. Zadaci 


16.10.1. Pretpostavljajući A = 1 provedite homogenizaciju rubnih uvjeta 
u' (0) =m,—u"O=quD="(D)=0. (30) 


Rješenje. Neka je h(x) polinom 3. stupnja koji zadovoljava uvjet (30). Zbog 
uvjeta u x = / imamo 


h()=a(x—1)?2+B(x—1)?, (31) 

a zbog uvjeta u x = 0 vrijedi 
2a == 6 IB = ma, (32) 
—6 8 = q. (33) 


Iz toga se jednoznačno određuju u i 8. Rješenje problema bit će 


u=v+h, (34) 

gdje je v rješenje problema 
vi kbv=f—bh, (35) 
vO=v"O)=vO=v(D=0. (36) 


16.10.2. Riješite rubni problem za štap u slučaju b = 0 uz ove rubne uvjete: 


(2) oba kraja učvršćena šarkama ; 

(1) oba kraja ukliještena; 

(122) lijevi kraj ukliješten, desni slobodan ; 

(20) lijevi kraj ukliješten, desni učvršćen šarkom. 
Rješenje. Iz (15.40) dobivamo 


(uy e)=[f&dE+C, (37) 
0 


(u) = [dn[|FG)AE+Ca+C,= [(2— DfE)AE+Cix+C,, (38) 
0 0 0 2 


bio \3A0 i ih n.dn 
DECE oroe+a [140 fra. (9) 
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x 


Piko Pata lo- oroa+o |, Dan + 


x 


+C, Io + Cax + Ca, 
gdje su C,, ..., Ca konstante. 
(2) Iz ž 

u' (0) =4(0) = 0, 

u"D=uQ) =0 
dobivamo 

C, = Ci == 0, 
OPE Enes 
l 
1 
C, = + Jrztala — FE) dE + 
0 
gdje je 
1 1 
P= [7(2)dx, L = | x 7(2) dx. 
0 0 
(2) Iz 

(O) =u(0) =0, 

vO=u(D=0 
dobivamo 


C; = C4 =0, 
dok se C, i C, dobivaju kao rješenje linearnog sustava 
d1iCi +a12C2 =b, 


dnCi +42,C2 =», 


l 
: Zb 
lair A(x) 
0 


gdje je 
1 


eis ] 1 


0 


98 


1 


i zako 


I— 
dđ 
A mira 1 dn, 


(40) 


(41) 


(42) 
(43) 


(44) 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 


Za determinantu tog sustava imamo 


11 
tu —es=7| |K de dy > (max 4(9)-13 
0 0 


pa su C, i C, jednoznačno određeni. 


(iti) Iz 
"(0) =u(0) =0, 
(Au"YO=""0=0 
dobivamo 
C; = Cu, =0, 
Ci, = —rF, 
C, = L. 


(zv) Iz (47) dobivamo (49), a iz (42) 


l y t/ 
= 4 s 
C, u] Zagoro |roax+ [ira IF), 


C=—C+L— 
gdje je 


>0, 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 


(59) 
(60) 


(61) 


(62) 
(63) 
(64) 
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16.10.3. Pokažite da su uvjeti (27) # (29) dovoljni za postojanje rješenja rubnog 
problema u slučaju kad su svi rubni uvjeti dinamički i b == 0. 


Rješenje. Iz (Au"Y (0) = 0 i (37) dobivamo C, = 0. Iz (Au")(0) =0 i 
(38) dobivamo C, = 0. Formula (40) daje rješenje do na progib C;x + C., gdje 
su C, i Cx proizvoljne konstante. 


16.10.4. Odredite maksimalnu (apsolutnu) vrijednost kontaktne sile i kontak- 
tnog momenta teškoga homogenog štapa čiji su krajevi (1) učvršćeni šarkama, (77) 
ukliješteni. 


Rješenje. (2) Prema (37), (38), (43) i (44) dobivamo 
16) = —(4tY 6) =—og(+—), (65) 


m (x) = (Au) (x) = i x — ZZ x. (66) 


q postiže maksimum zax =0ix =, 


1 
alne = 55; . (67) 
me X l 
Im| postiže maksimum za x = 5 
o gl? 
mlaz = E (68) 
(#2) Prema (37), (38) i (50)—(55) dobivamo 
ti\! l 
= —deyw=-e (3) (69) 
_ og“ ogl ,ogi 
me= sr DR. (70) 


lql postiže maksimum zax = 0ix = 1/2, 


1 
ale =: (11) 
|m| postiže maksimum za x = 1/2, 
12 
lm Ike == SIE . (72) 


16.10.5. Postavite i uz pretpostavku b = 0 riješite problem ravnoteže štapa 
koji je na lijevom kraju ukliješten, a na desnom slobodno oslonjen (sl. 7). 
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UA 
Slika 7. 


Rješenje. Uvjeti na lijevom kraju su 
u(0) = u (0) =0. 613) 


Za oslonjeni kraj postoje u ravnoteži dvije mogućnosti: ili je on u kontaktu 

s osloncem ili nije. U: prvom slučaju je.u (2) = 0, q (1) > 0; u drugom slučaju je 

u (1) > 0, a (1) = 0. U obama slučajevima je m (1) = 0. Prema tome u svakom. slu- 
čaju je 

u'(D=0, (74) 


u (1) >0, (Au) (D < 0, u (D (Au) (D) =0. (75) 


Uvjeti (75) su jednostrani i oni su nelinearni. Moramo razmotriti dvije mogućnosti. 


(0) (Au"Y (D) =0. Kraj je tada slobodan pa je rješenje dano formulama (40), 
(59) —(61). Nužno je u (1) > 0, iz čega dobivamo 


1 l 
1—y 
Jage le — (2) ni m 


9 
(7) u(l) = 0. Kraj zadovoljava uvjete šarke, pa je rješenje dano formulama 
(40), (62)—(64). Nužno je (Au")' (1) < 0, iz čega dobivamo 
1 1 ' 


i 
]aoje] 6-96 dx <0. (77) 


Vidimo da o tome koja će od dviju mogućnosti biti na kraju x = / realizirana 
odlučuje izraz na lijevoj strani u (76) odnosno (77). Dovoljan uvjet da vrijedi (76) 
odnosno (77) je, na primjer, nenegativnost odnosno nepozitivnost ove funkcije 
(varijable y): 


l HA l 
[ (a—9)f()dx= [ af (2) da —y [f (1) dx. (78) 


3 


$ 17. Greenova funkcija i rješenje rubnog problema za 
štap 


Razmatrat ćemo -rubni problem za štap pretpostavljajući da odgovarajući ho- 
mogeni problem ima samo trivijalno rješenje. Prema $ 16 to je ekvivalentno pretpostav- 
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ci da je ili b E O ili da su bar dva rubna uvjeta geometrijska, od kojih bar jedan sadr- 
ži progib. 


Označimo sa F (x) vanjsku linijsku silu koja djeluje na komad (0, x) štapa; 
tada na komad (x, x2) djeluje linijska sila 


F(x2) — F(x,). (1) 
Umjesto (15.10) dobivamo općenitiji zakon ravnoteže 
x2 : 
(Au"Y' (x2) — (Au"'Y (x) + [ bud =F(x2) —F(x,). (2) 
Xi 


Kao i u slučaju žice (8 8), jedinična sila koncentrirana u točki x' € (0, 7) je fun- 
kcija 


x: I 0,x< 
o=kio o 
Uz takvu silu iz (2) kao i prije dobivamo 

(A (x)u" ()" +b(x)u(x) =0 za x £ x'. (4) 


Pretpostavimo da su u točki x' neprekidni progib, rotacija i moment, to jest fun- 
kcija ui njezina 1. i 2. derivacija: 


u(x' +0) —u(x' —0) = 0, (5) 
u(x +0) —vw(x' —0=0 (6) 
u'(x' +0)—u"'(x' —0=0. \ (7) 


Iz (2) i (7) slijedi da 3. derivacija funkcije u (kontaktna sila) ima u točki x' skok: 


u"! (x sA 0) —u" (x 2 0) = Pa (8) 


Funkciju u (x) koja za dano x' €(0,/) zadovoljava uvjete (4)—(8) i zadane rubne 
uvjete označimo sa G (x, x"); ako ovdje i x' smatramo varijablom, dobivamo Gree- 
novu funkciju za Štap. 


17.1. Primjer 


Odredimo Greenovu funkciju u slučaju A = 1 i b = 0, ako je lijevi kraj štapa 
ukliješten, a desni slobodan. Imamo ove uvjete: 
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G (0, x") = 


Iz (9) slijedi 


G (xx) = 


*G (x, x") 
x+ 


G(x'+0x)—G—0x)=0 


=0xEx 


2G(x +0x)_2G(—0,2) _4 


Ox 


02G (x' +0, x) še 0?G (x' — 0, x) _ 


0x 


0x? 


PG POx) Vo —Da) 


0x 


0x 


0 x2 


0.x2 


gdjesu C;i D; (€ =0, ..., 3) funkcije od x“. Neka je 


CG, — D, = B, i=0,..,3. 


Iz (10) — (13) slijedi 


Iz toga dobivamo 


By = 


B4 + Bux -- B.x'? -- B, 3 = 0 


B,+2B,x' +3B3x'2=0 


B,+3B;x' =0 


Iz (14) i (21) slijedi 


6B; = —1 
bj HE 
radi are B, = 

Do=-<#% Di=3x" 


Iz (16) dobivamo 


čR=dt=ti a = 


2 


%, C, = 


0 


1. 


2G(0,x') _ 82G(1x) _ BG(1,2") _ 
Eda o Tagan 
Ca + Cix + Co? + C42%, x < x 
Do+Dx+Dx" +D32%,x>x, 


“- 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 
(17) 
(18) 
(19). 
(20) 


(21) 


(22) 


(23) 
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Prema tome je 


2 2x —x),x<g 
G(x,x) = , (24) 

1 12 ! I 

be Sx—x), x>x'. 

Funkcija (24) je simetrična, to jest vrijedi 
G(xx)=G(x,x za xx e€(0,); (25) 
ona je također pozitivna: 

G (xx) > 0 za x, x' € (0,7). (26) 


Svojstva (25) i (26) ima Greenova funkcija i u općem slučaju. 


Uz pretpostavku jedinstvenosti 1 uz homogene rubne uvjete rješenje rubnog pro- 
blema za štap dano je formulom 


l 
u (x) = | G (2,2) f(2') dx“. (27) 
i | 


Precizna definicija Greenove funkcije za štap potpuno je analogna 
onoj u slučaju žice (89). Posebno, pretpostavlja se da za svako x' € 
(0, 7) vrijedi 

G( ho? (0,0) NEO K) e D (28) 


Prelazimo na konstrukciju Greenove funkcije u općem slučaju. Pri- 
mijetimo najprije da rubne uvjete za štap možemo zapisati u obliku 


3 
U,(uj= rav) he? =0:=0123 (29) 
i=o 


gdje su a;;i B;; zadani brojevi. Prema $ 16 za uvjete na lijevom od- 
nosno desnom kraju je f;; = 0 odnosno a;y = 0. Međutim, za kon- 
strukciju koja slijedi takve pretpostavke nisu nužne, pa naše daljnje 
zaključivanje vrijedi za rubne uvjete mnogo općenitije nego one koje 
smo naveli u $ 16. Pretpostavit ćemo da homogeni problem 


(A"') +bu=0 (30) 
U (W)=0:1=0,1,23 (31) 


ima samo trivijalno rješenje u = 0. Neka su 44, R =0, 1, 2, 3, linearno 
nezavisna rješenja jednadžbe (30) na intervalu (0, 7); tada je matrica 
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(U; (249): =0,..3 (32) 


regularna (inače bi postojala netrivijalna kombinacija rješenja u, koja 
bi zadovoljavala uvjete (31), to jest postojalo bi netrivijalno rješenje 
problema (30), (31))*. Greenova funkcija je nužno oblika 


3 
X CGuG)x<x 
G(x,x) = ša (33) 
by Di u; (6), X > X 
i=0 
gdje su C; i D, neke funkcije od x" € (0, 7). Stavljajući 
C, = D, == B, (34) 
i uzimajući u obzir uvjete (5)—(8) dobivamo 
3 
Bu (x) =0 (35) 
i=0 
u(x)=0 (36) 
3 
xBuj (x) =0 (37) 
i=0 
3 B mis (38 
i s x Zave I\* : 
KP 69 — To) 2 
Uvjeti (29) daju 
3 
p> U; (u) D, = — L P> mi ui & (0)) B1i=0,...,3. (39) 
=0 


Uvjeti (35)—(39) predstavljaju sustav od osam linearnih jednadžba 
s osam nepoznanica B, D, £=0,...,3. Pokazat ćemo da taj sustav 
ima jedinstveno rješenje. Diremtinsnia sustava (35)—(38) je Wronski- 
jan (linearno nezavisnih) rješenja %, ..., ua jednadžbe (30), pa se na 
segmentu [0, /] nigdje ne poništava (v. Dodatak, $ 2). Odatle slijedi 
da sustav (35)—(38) ima jedinstveno rješenje B,, #=0,...,3. Time 
(39) postaje sustav s nepoznanicama D;, # = 0, ..., 3; njegova matrica 
(32) je regularna, pa on ima jedinstveno rješenje B;, D;, 1 =0, ...,3 


* Vrijedi i obrat: ako je za bilo koja četiri linearno nezavisna rješenja 1, k = 0,..., 3 jednadž- 
be (30) matrica (32) regularna, problem (30), (31) ima samo trivijalno rješenje u = 0; ako je rang 
te matrice jednak r < 4, problem ima 4-r linearno nezavisnih rješenja. 


105 


106 


Iz (34) dobivamo C;, € = 0, ..., 3. Tako određena funkcija (33) zadovo- 
ljava sve uvjete koji se traže od Greenove funkcije. Rješenje sustava 
(35)—(38) lako nalazimo pomoću Cramerovog pravila; koristeći se još 
Liouvilleovom formulom (v. Dodatak, 2), dobivamo 


_ O WEE NAGE. 
B-- POVO 7%? (40) 


gdje je W;(x'") algebarski. komplement elementa u; (x) u determi- 
nanti W (x) = W (uo, ..., 435%); bez smanjenja općenitosti možemo 
pretpostaviti da je 


A2(0W(O) = —1, (41) 
pa dobivamo 
Bi = A (29) W,(2). (42) 
Iz toga slijedi da je 
By, Ce Di € 32 (0,1), 120;:45% 3. (43) 


Dokazat ćemo još formulu (27). Prema (33) imamo 
3 ž l 
u(s)= Su (9 (| Da)fG)dx+ [ GG) f(a)ds'). (44) 
i=0 G : 


Neka je x € [0, /] točka neprekidnosti funkcija A", bi f. Uzimajući u 
obzir uvjete (35)—(38), iz (44) dobivamo 


(du) (9 => A" G9([ DG) F6) dx + 
Ša 0 


l 
+ [ G(B)] (0) dx") +(0) (45) 
Zbog jednakosti 
(Au ()= —b(uE)i=0,..3 (46) 
iz (45) slijedi 
(Au")" (g) +b (Yu) =7(2) (47) 
Dalje imamo 
l 
U, (= | U(GC,2) fx) dx —0, (48) 
0 


jer je U(G(.,2)) = 0 za svako x' € (0,7). 


Rezultati ovog paragrafa prenose se s nebitnim promjenama na 
rubni problem za linearnu diferencijalnu jednadžbu n-tog reda (n > 2) 


ao(x)uUP (x) +a(H)u"-P (x) +. traku) =f(9 (49) 


na segmentu [0, /]; ovdje su a4, 41, ...> 4, i f zadane, a u nepoznata 
funkcija. Pretpostavlja se 


a9(x) £0zaxel[0,/]. . (50) 


Rubni uvjeti za jednadžbu (49) glase 
n—1 
U; (W) =ža;40 (0) + Bu Q=01=0,1.,n—1. (51 
j=0 


17.2. Zadatak 


Dokažite korektnost rubnog problema za štap. 
Rješenje. Iz (27) slijedi 


l 


lu (2) < max |f(2)lmax | |G(x,1)lds', (52) 
xel[0,/1 xe[0,7] < 
lu" (2)| malj (| max [lee s: s) | dx, (53) 
xe[0,7] 


lu" (|< ni x 1769) max [|2260 God dx. (54 


xe[0,f] 2 0x2 


8 18. Varijaciona jednadžba i funkcional energije za štap 


Varijaciona i energetska formulacija rubnih uvjeta za štap dobiva se na isti 


način kao u slučaju žice. Zato se ograničujemo na izvođenje osnovnih formula, a 
dokaze prepuštamo čitatelju. Također se ograničujemo na rubne uvjete 


uQ)=w"(0)=0 (1) 


u" (1) =(Au"Y(D=0; (2) 


postupak je isti i u drugim slučajevima (u kojima je rješenje jedinstveno). Za fun 
kciju v na segmentu (0, 7) kažemo da je dozvoljena ako zadovoljava zadane (homoge 
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ne) geometrijske uvjete rubnog problema*; u našem slučaju to znači v (0) = v' (0) = 
= 0. Pomnožimo li jednadžbu (15.40) s takvom funkcijom, integriramo po seg- 
mentu [0, 2] i iskoristimo formulu (16.21) te uzimamo u obzir (1) i (2), zaključujemo 
da vrijedi 


1 
| (Au" 0" + buv'— fv) dx = 0 za svako dozvoljeno v. (3) 
0 


Prema tome, akoje funkcija u rješenje problema (15.40), (1), (2), ona zadovoljava va- 
vijacionu jednadžbu (3). Pokazuje se da i ovdje vrijedi obrat: Ako dozvoljena fun- 
kcija u zadovoljava varijacionu jednadžbu (3), ona je rješenje problema (15.40), (1), 
(2). 


Funkcional energije ima oblik 


4 H 
D(w) = = | (A ("")2 + bzv2) dx — | fwdx, (4) 
0 0 


a promatra se na dozvoljenim funkcijama w. Pokazuje se da rješenje u rubnog problema 
(15.40), (1), (2) minimizira funkcional energije (4) (na skupu dozvoljenih funkcija) i 
obratno, da funkcija u koja minimizira funkcional energije (4) (na skupu dozvoljenih 
funkcija) predstavlja rješenje rubnog problema (15.40), (1), (2). 


Gornja razmatranja prenose se bez bitnih promjena na rubni 
problem za linearnu diferencijalnu jednadžbu reda 2m(m > 2): 


že Lj (aphjneb G) n = s 


U ovom slučaju geometrijski (Dirichletovi) rubni uvjeti su oni u kojima 
se pojavljuju derivacije funkcije u reda najviše m — 1; ostali rubni 
uvjeti su prirodni (Neumannovi). : 


18.1. Zadatak 
Napišite varijacionu jednadžbu i funkcional energije za štap za ove slučajeve 


rubnih uvjeta: 


(2) šarka na lijevom, a žlijeb na desnom kraju, 

(72) šarka na oba kraja, 

(22) lijevi kraj slobodan na progib i elastično vezan na rotaciju, a desni kraj 
slobodan na rotaciju i elastično vezan na progib. 

Rješenje. (2) Dozvoljene funkcije: v (0) = v' (2) = 0; varijaciona jednadžba 
je (3), a funkcional energije (4). (27) Dozvoljene funkcije: v (0) = v (7) = 0; varija- 
ciona jednadžba je (3), a funkcional energije (4). 


* Ovdje se pretpostavlja u € #2 (0,1). 
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(21) Dozvoljene su sve funkcije: varijaciona jednadžba: 


I u 


[ (Au"v" + buv) dx = ) fod— u (0)V(0) +24 (2) 0 (1); (6) 
0 < 0 


funkcional energije: 
1 1 
b(w) = 3 | (Aw"2 + bož) dx — | fudx + Zo (0))2— ZW). 
0 0 ' 


S 19. Metoda konačnih elemenata za štap 


Pozabavit ćemo se približnim rješavanjem problema ravnoteže štapa. Izabrat 
ćemo tipične rubne uvjete 


u(0) = u (0) =0 .(D 
u D=(Au""YQD+s(D=0, (2) 


što odgovara varijacionom problemu: odrediti dozvoljenu funkciju u (to jest onu 
koja zadovoljava uvjete (1)) takvu da je 


1 1 
je (Au'"v" + buv) dx = [ (fv + sv") dx za svako dozvoljeno v. (3) 
0 0 


Razdijelit ćemo segment [0, /] na x» + 1 jednakih dijelova (duljine h = 1/(:-- 1)) 
ČVOVOVIMA X, X25 .-.> Xn (Stavljamo xg = 0, X,.g1 = 1). Na žalost, ovdje nije moguće 
primijeniti skup V ; iz $ 14, jer se u (3) traži da u i v imaju konačne druge derivacije; 
to za funkcije iz V,, općenito nije slučaj (funkcija 9, ima u ČVorovima X;_,, X) *k41 
beskonačnu drugu derivaciju). Zato ćemo koristiti po dijelovima kubične funkcije. 
Funkcija w je po dijelovima kubična za razdiobu Xi, X2> ...> Xn ako je na svakom od 
diobenih intervala (2%_ ,, Xx) polinom stupnja najviše 3 i ako je u svakom čvoru xy, 
(k =1,2,...,1) neprekidna zajedno s derivacijom (sl. 1). Skup svih-takvih funkcija 
označujemo sa V? (0, 7) (kraće sa V2). Proizvoljna superpozicija funkcija iz V: jest 
funkcija iz V*. 


| - x 
O=xo Xi X2 *n l=Xn41 Slika 1. 


* To znači da je V? linearni prostor. Očigledno je VŽ e Sf? iV?2 cl: 4+. 
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Polinom stupnja 3 ima oblik 
a0%5 + a2 + ax + 03. (4) 


Za određivanje koeficijenata a4, ..., 4; potrebna su četiri podatka. Posebno, ti koe- 
ficijenti su jednoznačno određeni vrijednostima polinoma (4) i njegovih derivacija 
u dvjema različitim točkama (v. zadatak 19.1.1). Iz toga slijedi da je funkcija zv e 
e Vi; potpuno određena vrijednostima 20 (xx) i 0 (4), k =0,1,..,n +1. 

Svakom čvoru xx» &k=0,1,..,n +1 pridružit ćemo par funkcija W.y 
V+1 € V2, koje imaju ova svojstva : 


0,j#k 


Pa (X) zi 1 j = 


> Zi (;) =0 za svako j, (5) 


0,j/zk 


Vax+1 (X) = 0 za svako j, Paua (29) = Pr (6) 


(sl. 2). "Te funkcije su dane formulama 


"x-x(25%) (D 
Van) =hx (ž ; x), (8) 
0 Paket Slika 2. 
gdje je 

(lx] — 1)*(21x1 +1), xe[—L1] 
2 e e vr < 

x (lx] — 1)? > xe[—l,1] 
X (2) -( 0 oSAELI (10) 


(sl. 3). Funkcije Y/,, # = 0, 1, ..., 2n ++ 3 zovu se hermitski kubični elementi. Svaku 
funkciju zv e V2 možemo prikazati kao superpoziciju funkcija P,: 


Slika 3. 
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2n+3 n+1 
vw (x) = že "09 mao Pax (2) + azkga Poznan (2), (11) 
gdje je 
ax =), dara = U (4). (12) 
Zaista, zbog (5) i (6) funkcije na lijevoj i desnoj strani u (11) podudaraju se zajedno 
s derivacijama u svim čvorovima pa su jednake. Prikaz (10) je za funkciju we V; 


jedinstven, jer zbog (5) i (6) iz (11) slijedi (12). Dozvoljena funkcija zv e V? ima pri- 
kaz 


2n+3 n+1 
wx) = > PH(2) — Dva VP, (£) + Gra Pong (0). (13) 


Kažemo da je dozvoljena funkcija # e V: približno rješenje problema (3) ako 
vrijedi 


1 1 
[ (Av + biv) dx = | f vdx za svako dozvoljeno v e Vi. (14) 
0 0 
Funkcije # i v pišemo u obliku 
2n+3 
ž()=Xa (2), (15) 
r=2 
2n+3 
v (x) == 2 # # (2), (16) 


pri čemu koeficijente a, treba odrediti, a koeficijenti f, su proizvoljni. Posve ana- 
logno kao u $ 14 iz (14)—(16) dobivamo 


2n+3 
>O. =h$=23..,2 48% (17) 
gdje je 
C = A “E bsp (18) 
l i 
da = [A PP" dx, o, = [09v W, dx, (19) 
0 0 sa 
1 m 
f=[10,d6 (20) 
0 


C =(C,,) je matrica krutosti. Sustav (17) sadrži 2n -- 2 linearnih jednadžba s isto 
toliko nepoznanica a, ..., Q2n+3- Matrica C je simetrična. Ona je i regularna*, 
pa sustav (17) ima jedno i samo jedno rješenje (a2, &3, ...> 42,+3)- Pri tome je prema 


* Matrica C je pozitivno definitna. Dokaz je posve analogan onom u $ 14. 


» 


lll 


(12) ax = 8(X4) dogri = #' (X,), pa za problem (3) dobivamo aproksimaciju rje- 
šenja i njegovih derivacija u čvorovima x,, x2, ..-> Xnr1*. 

Iz (19) lako zaključujemo da je C sedamdijagonalna (to jest da ima samo sedam 
dijagonala različitih od nule) i da ima ovu strukturu: 


X X X X 
X X x X 


(21) 


(XX XXX 


x 
pas 


XX xXx XxXxX X 


Đ XXX xX 
i DOKA 


19.1. Zadatak 


Dokažite da je polinom 3. stupnja određen svojim vrijednostima i vrijednostima 
svojih derivacija u dvjema različitim točkama. 


Rješenje. Neka je P(x) = ao2? + 0127 + a,x + a3, O (x)= Boi B,2x2 + 
+P2x+Bsinekaje P(x) =0(x), P'(x) = 0 (x) 1= 1,2, gdje je x, # x». 
Tada za polinom T= P— Q vrijedi T(x) = T'(x) =0,1: = 1,2. Prema tome, 
x, i x, su dvostruke nultočke polinoma T koji je najviše 3. stupnja. Iz toga slijedi 
T=0. 


8 20. Stabilnost štapa 


Ako je uzdužno djelovanje na štap slabo, ono ne utječe na poprečni progib 
(a = q, ne ulazi u (15.37) —(15.39)). Ako je to djelovanje nenegativno, ono također 
ne izaziva poprečni progib. Negativna i po apsolutnoj vrijednosti dovoljno velika 
longitudinalna sila izaziva poprečni progib (izvijanje štapa). Proračun progiba u 
tom slučaju stvar je nelinearne teorije. Međutim, u sklopu linearne teorije, koristeći 
se modelom koji smo izveli u zadatku 15.1.2, moguće je proračunati pri kojoj lon- 
gitudinalnoj sili će doći do izvijanja. Razmatranja su slična onima u $ 7. Polazimo 
od jednadžbe (15.48). Neka je a= —pu, u = const. > 0 (sl. 1). Homogena jednadžba 
(15.48) glasi 


(Au) + uu" = 0. (1) 


Pretpostavljat ćemo da su uz tu jednadžbu dani homogeni rubni uvjeti i daza u =0 
odgovarajući rubni problem ima samo trivijalno rješenje. Očekujemo da za dovoljno 
malo u > 0 problem također ima samo trivijalno rješenje. Najmanja vrijednost pa- 
rametra u > 0 za koju problem ima netrivijalno rješenje zove se kritična sila; obilje- 


Žavat ćemo je sa ui. 
x 
u a ananaičinak -p ' 


Slika 1. 


* Slično kao u & 14 dokazuje se i ovdje da približno rješenje minimizira funkcional energije 
(18.4) na prostoru V?. 
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Ako jeu < u,, ravnotežni položaj je u = 0. Za u = u, ravnotežni položaj u = 
= 0 je nestabilan, jer postoji i ravnotežni položaj u, s 0 (izvijeni štap). U daljnjem 


pretpostavljamo da je štap cilindričan i homogen (E = const.). 


20.1. Primjer 


Odredit ćemo kritičnu silu za štap kojem je lijevi kraj ukliješten, a desni učvrš- 


ćen šarkom (sl. 2). Neka je uJA = A. Imamo problem 
(ud) + 2" =0 
u(0) ="(0) =0 


ud)=u"(D=0. 


Slika 2. 


Opće rješenje jednadžbe (2) je 
u(&) = C,x + C, + Gz cos /2x + Ci sin /2 x, 


gdjesu C,, ..., Ca proizvoljne konstante. Iz (3) slijedi 


ča tita 2 


Va 


&. 
Iz (4), (5) i (6) dobivamo 


(1 e ini 1Va) C,+(1—cos2/4) C, =0 


Va 
VasiniV/A +C, +AcoslVA +C, =0. 


Taj sistem ima netrivijalno rješenje (C,, C») ako i samo ako je 


1— 1 sini VT 1—coslj/4 
M _o 


VA sini VA AcosiV2 
č= tg €, 


Iz toga dobivamo 


8 Jednadžbe matematičke fizike 


(2) 
(3) 
(62) 


(5) 


(6) 


(D) 


(8) 


(9) 


(10) 
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gdje je 
E=IVA (11) 
Najmanji pozitivni korijen jednadžbe (10) je €, = 4,49 (sl. 3). Dobivamo 


&A 
m =" (12) 


Odgovarajuće rješenje je 


(13) 


Ho: : 
u) = c( — g, 8m-r — 2/sin Zi 


dark 55) 


gdje je C proizvoljna konstanta (sl. 4). 


Slika 4. 


20.2. Primjer 


Odredit ćemo kritičnu silu za štap kojem je lijevi kraj ukliješten, a desni slo- 
bodan (sl. 5). Neka je u|A = A. Imamo problem 


ud) + žu =0 (14) 
u) ="(0) =0 (15) 
uD=0 u" D+" D=0; (16) 


ll4 


Slika 5. 


primijetimo da se u ovom slučaju parametar A pojavljuje i u rubnom uvjetu (v. 
(15.47)). Za 4 dobivamo uvjet 


cos 1 /A = 0. (17) 
Iz toga slijedi 
IVA = (2n— DZ n=12.. (18) 
ili 
jo =) 21 (19) 
Najmanja vrijednost dobiva se za n = 1, pa imamo 
n= a (20) 
Odgovarajuće rješenje je 
u(x)=C (1 _ cos 3) (21) 


gdje je C proizvoljna konstanta (sl. 6). 


0 ' => x 


Slika 6. 


- 


Zadatak o kritičnoj sili je još jedan primjer problema svojstvenih vrijednosti. 
Svojstvena vrijednost jednadžbe (1) (uz dane homogene rubne uvjete) je ona vrijednost 
parametra u za koju ta jednadžba ima netrivijalno rješenje; odgovarajuće rješenje je 
svojstvena funkcija. U primjerima 20.1 i 20.2 svojstvenih vrijednosti ima prebro- 
jivo mnogo; to su (do na faktor) kvadrati korijena jednadžbe (10) (v. sl. 3) odnosno 
vrijednosti (19). Kao rješenje homogenog problema, svojstvena funkcija je određena 
do na konstantan faktor. Problem svojstvenih vrijednosti za jednadžbu (1) pojavit 
će se i u veži s titranjima štapa. 


20.3. Zadaci 


20.3.1. Odredite kritičnu silu za štap kojem su oba kraja učvršćena šarkama 
(sl. 7). 
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Slika 7. 


Rješenje. Imamo problem 


ut) k2u=0 ' (22) 
u) =2"(0) =0 (23) 
w="0=0, (24) 
gdje je A = ul A. Dobivamo 
n="f (25) 
Odgovarajuće rješenje je 
u (x) = Csin S, ' (26) 


gdje je C proizvoljna konstanta (sl. 8). 


Skika 8.. 


20.3.2. Štap kojem su oba kraja učvršćena šarkama rotira jednoliko oko cen- 
tralne linije kutnom brzinom 2. Odredite kritičnu vrijednost te brzine (to jest naj- 
manju vrijednost pri kojoj je ravnotežni položaj u = 0 nestabilan). 


Rješenje. Na štap djeluje centrifugalna sila s gustoćom (sl. 9) 


f[&)=72 [[ Ges) +uG)dyde= 20746) (27) 
D(x) 


Slika 9. 
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i spin s gustoćom 


s(x) = £2 [| 9sino (0) eoso() +u(s))dydsm 


D(x) 


2 so [|yowo+uG)drd= 


D(x) 


205 [] yu +uG)dyda= LOŠI u (s 
D(x 


(28) 


ovdje je e linijska gustoća mase (zbog homogenosti je g = const.), a S površina 


poprečnog presjeka. Jednadžba (15.39) glasi 


0: 


ui u -5« 


gdje je 


Uz oznake 


m fra 1+ V1+£5), r-Vžl+ Vi mi 


opće rješenje jednadžbe (29) je 


u()=Circhax+C,shax + C; cos Bx + Cusin B x, 


gdje su C,, ... Ca proizvoljne konstante. Uvjeti u (0) = u" (0) = Odaju C, = 


= 0. Uvjeti u (1) = u" (2) = 0 daju 


\ 


Cosh al -+-Cysin BI = 
C,a?shal — C,B7sin fl =0. 
Taj sistem ima netrivijalno rješenje ako i samo ako je 


shal -sinB1 -+(a2 +87) =0. 
Iz toga slijedi 
snpi=0 
ili 
Pil=nmn=1,2,.. 


(29) 


(30) 


(31) 


(32) 


C; = 


(33) 


(34) 


(35) 
(36) 


(37) 
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Imamo jednadžbu 


Vk slt+ +55) na (38) 


Najmanji korijen te jednadžbe dobiva se za n = 1 i jednak je 


av deva“ 
Kritična kutna brzina je 


VE = (VV. M, 
2, Ki (5) VG) (40) 


20.3.3. Štap kružnog presjeka je ukliješten na obama krajevima. Lijevi kraj 
je učvršćen na torziju, a desni zakrenut (oko osi x) za kut r (sl. 10). Odredite Rrizičnu 
vrijednost tog kuta (to jest najmanju vrijednost pri kojoj nedeformirano ravnotežno 
stanje u = 0 postaje nestabilno). 


i 


Slika 10. 


Rješenje. Ako štap nije izvijen, za longitudinalni kontaktni spin dobivamo 
m,= Cr, C=u -*R*a|/21 (v. zadatak 1.2.2). Za izvijen štap prirodno je usvojiti 
taj zakon za komponentu spina u smjeru tangente na izvijenu centralnu liniju. Tada 
umjesto (15.29) i (15.30) imamo 


m =EIu6+Cru (41) 
m =—EIui+Cru (42) 


gdje je /=1I, = I,. Jednadžbe ravnoteže glase 


u +lu"=0, (43) 
us Eas 1 uy = 0, (44) 
gdje je 
CT 
ir (45) 


Sustav (43), (44) riješit ćemo uvođenjem nove funkcije 
U=u,+iu,. (46) 
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Ona zadovoljava jednadžbu 
Geo —:1U "= (47) 


Opće rješenje te jednadžbe je 
U = C,x2 + Ca + C3 + Cd, (48) 


gdje su C,,..., C, proizvoljne konstante. Uvjeti U (0) = U' (0) = 0 daju Cz = 
= —C4 Ca= —i2C,. Uvjeti U (1) = U (1) = 0 daju 


PC, +(—i4l—1+e2)C,=0, (49) 
NC, +(—iA +12 04) C, =0. (50) 


Taj sustav ima netrivijalno rješenje (C,, C4) ako i samo ako je 


2+:241 
il _ 
Seanu GI 
Iz toga slijedi 
X 11 
gi=bi=7 (52) 


Najmanji pozitivni korijen te jednadžbe je ć, 24,49. Kritična vrijednost kuta je 


EI, _EI2 
LAT 


Ti = 


Ši (53) 


S 21. Nelinearni problemi 


Pozornijim promatranjem izvoda naših rubnih problema u $1 i 6 15 uvjerit 
ćemo se da je linearnost posljedica pretpostavke 


(a) da su deformacije male, 
(b) da svojstva materijala (na primjer, modul elastičnosti) ne ovise o progibu. 


Linearnost omogućuje — bar načelno — elegantno i pregledno rješavanje rub- 
nog problema, o kojem smo govorili u prethodnim paragrafima. 

S druge strane, ako nisu zadovoljene pretpostavke tipa (a) odnosno (b), pot- 
rebno je posegnuti za nelinearnim modelom, za koji se pretpostavlja da bolje odgova- 
ra stvarnosti. Isto tako, matematičare će zanimati pitanje koliko je izvod linearnog 
modela bio korektan. To znači da treba strogo ocijeniti grešku linearnog modela pre- 
ma točnijem nelinearnom. U svakom slučaju nelinearni model najprije treba posta- 
viti i riješiti. Teorija je ovdje, dakako, mnogo zamršenija, a isto tako i metode prak- 
tičnog rješavanja. I jedno i drugo je još uvijek predmet intenzivnog proučavanja kako 
matematičara tako i inženjera, iako su neke klase problema već dobro proučene. 
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U $ 13 vidjeli smo da problem minimizacije nekvadratičnih funkcionala oblika 
(13.1) vodi na nelinearne rubne probleme. Sustavno izlaganje metoda rješavanja ne- 
linearnih jednadžba daleko bi prelazilo okvire ove knjige. Mi ćemo se stoga ograni- 
čiti na jednostavnije slučajeve, posebno one koji se ne udaljuju mnogo od linearnosti. 


21.1. Primjer 


Promotrimo jednadžbu provođenja topline 


—_auwY=f (1) 
uz rubne uvjete 
u(Q) =0,2(]) =d (2) 


i uz pretpostavku da je koeficijent provođenja a funkcija kako mjesta x tako i same 
temperature u, tako da na mjestu x koeficijent provođenja ima vrijednost a (u (x), x). 
Pretpostavljamo 


a(u,x) > 0 (3) 
za sve x € [0, /] i za sve u. Integracijom iz (1) dobivamo 


ii — FG G 
i JEŠJE TE (4) 


gdje je 


Pol== | fd, (5) 
0 


a konstantu C treba odrediti iz rubnih uvjeta. Jednadžba (4) je obična diferencijalna 
jednadžba prvog reda, koja se općenito ne da svesti na kvadrature*, tj. rješenje jo 
se ne da prikazati pomoću neodređenih integrala zadanih funkcija. Osim toga nije' 
jasno postoji li njezino rješenje i je li jedinstveno. 

Ako je štap komogen, onda a ne ovisi o x, pa iz jednadžbe (4) slijedi 


u 


[aQ«)du=([ F&AE+Cx+C,. (6) 
0 


0 
Vidimo da smo sad jednadžbu (1) »sveli na kvadraturet, pri čemu jednakost (6) 


određuje u zmplicitno kao funkciju od x. Konstante C i C, određuju se iz rubnih 
uvjeta. Uzmemo li, na primjer, 


a(u) =a, + u, (7) 


gdje su a, i € pozitivne konstante, prelazi (6) u 


* U starijoj matematičkoj terminologiji »kvadratura« je sinonim za »integracija«. 
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dju aa a = s =/[r F (£) dE + Cx + Ci. (8) 
0 


Uvrštavanje rubnih uvjeta (2) daje 


1 
c= rlad+2f— | Fe) a),c—0. (9) 
l 3 

0 


Zbog pozitivnosti konstanata a, i s moguće je za svako x naći točno jedno u tako da 
bude ispunjena jednadžba (8), pa je time definirano rješenje u. 


21.2. Zadaci 


21.2.1. Dokažite da za svaku neprekidnu funkciju a = a (1) takvu da je 


= mfa(n) > 0. (10) 
diferencijalna jednadžba 
—(a(u(2)u)Y =1() (11) 


uz rubne uvjete (2) ima rješenje i da je to rješenje jedinstveno, 


Rješenje. Stavimo li 
u 


p(W=[a()dq (12) 


0 
vidimo da uvjet (10) povlači da je iznverzna funkcija p-1 defnirana na (—oo, 00). 


Tako iz (6) slijedi 
u() = 9-1(1(), (3) 
gdje je 
x l 
ne = | Feat+F00— | rea) a4 
0 


0 


Jedinstvenost dobivenog rješenja u slijedi iz same konstrukcije. 


21.2.2. Odredite rješenje u rubnog problema (1), (2) uz pretpostavku (7), 
akojeaa =e=d=1=1,f(x) = 1. Predočite to rješenje grafički. 


Rješenje. Imamo 


b2 


|roae- ča =5 (15) 
0 


0 
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i prema (9) 


4 1 5 
vise kae kir s (19) 
Iz (8) i (9) izlazi jednadžba 
u 2 5x 


koja zmplicitno definira funkciju u. Vidimo da je prema očekivanju # (0) = 0, u (1) = 


=laosimtogaiu (5) = 1/2. Koristeći se formulom (4) za derivaciju 


2 
x+ zA 
deri ue 
: ; 5 .,(1 16 _, 1 i : DRE 
dobivamo u' (0) = rak 5) = 1“ (1) = 7» 14 osnovi čega se može načiniti 
gruba skica krivulje u = u (x) (sl. 1). Točnija slika dobiva se iz prikaza inverzne 
funkcije 
4 5 E: uš 


Slika 1. 


Ovo je bio jedan od rijetkih slučajeva egzaktno rješivog nelinearnog rubnog 
problema. U sljedećem primjeru upoznat ćemo jednu vrstu rubnih problema, koji 
se mogu aproksimirati linearnim. 


21.3. Primjer (Nelinearna linijska sila.) 


Neka je f (u, £) funkcija varijabla # € R, € e [0,7]. Promotrimo rubni problem 
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—(a(2)u' (9) =f(u(2)%) (20) 


u(0) =0y»vw(7) +81(7) =0, (21) 


pri čemu za y i & pretpostavljamo uvjete (2.9). Kažemo da je ovaj problem slabo 
nelinearan ako je 


IT) =fQ0) —bO)u+/2(6 2%), (22) 


gdje je b 20, a funkcija f» mala po apsolutnoj vrijednosti. U tom slučaju imamo 
jednadžbu : 


—(auY +bu=f1(%) +f2(6 2). (23) 


Problem (23), (21) ekvivalentan je integralnoj jednadžbi 


I 
u) =g(0)+[GxD)fu(E)ĐdE, (24) 
0 


gdje je G odgovarajuća Greenova funkcija i 
1 
g()=[G(8)/, (dt. (25) 
0 


Funkcija (25) je rješenje »nesmetane« linearne jednadžbe 
—(au') +bu=f,(2) (26) 
uz rubne uvjete (21). Smatramo li f, »malom smetnjom«, uzet ćemo funkciju 
u =g (27) 


za prvu aproksimaciju rješenja jednadžbe (23). Za drugu aproksimaciju stavljamo 


l 
u2() =8()+ [ G(a,Đ)7, (u (8), &)dE, (28) 
0 
i slično 


l 
"ar =gB+ [GG D)f, (E) k=2,3,.. (29) 
0 


Očekujemo da će uy davati to bolju aproksimaciju točnog rješenja u što je k veće. 
Ovaj se postupak zove metoda sukcesivnih aproksimacija. Za praktične svrhe 


123 


ona je preporučljiva samo ako je fx toliko maleno da možemo očekivati da će već 
druga aproksimacija biti zadovoljavajuća te ako se integrali u jednakostima (25) i 
(28) daju eksplicitno izračunati. Ovaj posljednji zahtjev uključuje i eksplicitno poz- 
navanje Greenove funkcije. 
21.4. Zadaci 
21.4.1. Izračunajte prve dvije aproksimacije u primjeru 21.3, ako je 
P=08#0,/1Q) =f1 0=04=4, (30) 
ln(ux) = —eu, (31) 


gdjesuao >0,€ >0,i fo konstante. 


Rješenje. Greenovu funkciju nesmetanog problema znamo od ranije (v. 
zadatak 9.2.1): 


ži 1 (1-$)e«<: 


G(x&)=—. (32) 
“A (1 -3) ča>E 
U tom slučaju je nesmetano rješenje (prva aprok$šimacija) dano formulom 
m&)=eg=đ21-a)s, (33) 


a za drugu aproksimaciju imamo 


fo 


PJE 0-+ [obe BaE=f,(1—s)— 


ali (- 7 Pago kme doe a) 
ŠTA x S +3 Suki TE +5“ : (34) 


21.4.2. Ispitajte konvergenciju miza u,, 42, ... suRkcesivnih aproksi- 
macija definiranog formulama (27)—(29). 


Rješenje. Postupak je posve analogan dokazu konvergencije 
sukcesivnih aproksimacija u Dodatku. Iz nejednakosti 


l 
lee) — GL < | GG 8) [2 ((2), 8) — (35) 
0 
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ti 
—f2 (1 (8), ĐLAE < K | G(x, 6) dć max lx() — ug (20)| 
0 i 
uz pretpostavku da je 
K = sup 


se 
—oo<u<oo 


(36) 


= e 5) eĐl< o, 


zaključujemo da niz u,, 42, ... konvergira jednoliko na [0,7], ako je 


i 
8=K[G(x5)d6<1. (37) 
0 


Jednako kao i u Dodatku dokazuje se da niZ 41, 42, ... Ronvergira prema 
rješenju u jednadžbe (24) i da je to rješenje jedinstveno. Ocjena greške 
aproksimacije uy bit će 


max lu(9) — (91 < 38 max lux (9 — uu (9. (9 


Primijetimo da za funkciju f» oblika (31) supremum (36) nije 
konačan, pa prema tome naše razmatranje nije primjenjivo na taj 
slučaj ma kako mali bio broj e. 


21.4.3. Primijenite metodu sukcesivnih aproksimacija na rubni problem u 
primjeru 21.1 


Rješenje. Pretpostavljamo da je ovisnost funkcije a o varijabli u slaba. "To 
možemo izreći tako da stavimo 


1 


FICO NE a, 5 + a2 (u, X), (39) 


pri čemu je a, pozitivno, a a» dosta maleno. Integrirajući jednadžbu (4) i uzima- 
jući u obzir rubne uvjete, dobivamo 


(F&+C 

u(x) = FIČTINAN u&.9 dić, (40) 
0 
d- [-FOa_ 

ono. (41) 


o 4(u(£),6) 
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Mogli bismo uvrstiti (41) u (40) i tako dobiti jednu jednadžbu za nepoznatu fun- 
kciju u. Budući da je ta jednadžba dosta zamršena, rješavamo (40), (41) kao sustav 
jednadžba za nepoznato C i u. 

Budući da u formuli (39) i funkciju a, smatramo »malenom«, prvu aproksima- 
ciju rješenja definirat ćemo tako da u (40) odnosno (41) umjesto a uvrstimo a,. 
Tako dobivamo 


FF(&)dE 
d— [5 
G pm zm i &) (42) 
[= 
0 
u) = ngća dč. (43) 
0 


Sljedeće sukcesivne aproksimacije definiramo ovako: 


i F(&)di 
š- 1 a5 _ Ba(u(£),&) &) 
Kazan dč 
0 a(u(£)) 8) 


NLO 
Ui1 (9) = Ul EFICNGNJR de, (45) 


Ga = (44) 


zak=1,2,.... Primijetimo da tako definirane aproksimacije zadovoljavaju oba 
rubna uvjeta, to jest vrijedi x (0) = 0 i u, (1) = d za sve k. 


21.4.4. Ispitajte postojanje i jedinstvenost rješenja rubnog problema 
iz primjera 21.1 koristeći se sukcesivnim aproksimacijama definiranim 
u zadatku 21.4.3 za slučaj da je 


a (u, x) =a (x) + € u (46) 
gdjejeg >0i 


=infa(x) >0, 6 =supa(x) < o. (47) 
Rješenje. Prema (39) i (46) imamo 


guž 


026%) — — ZT Fen (9) 


Najprije ćemo dokazati da su nizovi C;, #4, ograničeni. Zaista, iz (42) 
i (44) slijedi 
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t 
lo < o 2 (4+ | 1re)la), b=lbi« (49) 
0 


a iz (43) i (45) slijedi 


l 
lu, (x)|< gebe -| IFE)ld&k=12.. (50) 
44 (X) 
0 


Pišući kraće a, (€) = a(u, (6), £), imamo 


1 
duct (z- l ) aš + 
4% 0% -1 


i 1 -2)a (51) 


i odatle 


[Cca — O] £ ŽEK(0) | IPIdE max (5) — 1 (91 + 
š 0Zx<i 


0Zxql 


1 
a (a ia) IF| de) (2) 1K (g) max [u (2) — ura (X)| = 


1 
2 KO (40 d + (1 +20) il IF| d& max |u, (2) — wx-1 (| = 
49 ž 0zx<l 
= K, (€) max | (Xx) — 4-1 (2)|. (52) 


Pri tome smo se koristili jednakošću 


1 za 1 | 2 
arta a8)+7(0 7 
< K odre lv (m — 22 (ml, (53) 


koja vrijedi za |v,| S U, |v2| S U, gdje je 


1 


2 2U 
dn a, (£) + en? 


< (59 


2en 
K (£) = max I 
2 nl<U (a9 + £n)“ 
O<E< 


< max 
hni<U 


1217 
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Mo=K 02(a ia (1 +2) | IF at). (55) 


Odatle i iz (45) izlazi 
: l 
lesa 6) — GL < K (2) [(1F] + Č) dč max lu (E) — ua (1+ 
0 == 


+2 Gea — GO] £ Ka (0) max ha () — ua (Bh (59 


gdje je 
K,(€)=K(£) (č1+(14+2 pa 5+ (2) JI IF| ae). (57) 
Rješenje postoji i jedinstveno je za 


K,()<1. (58) 


Pitanje korektnosti može se odnositi na dvije stvari: (£) ponašanje 
rješenja za e -> 0, to jest kad nelinearnost postaje slabija, (27) ponašanje 
rješenja za mala vanjska djelovanja, to jest za maleno F i d. U obama 
slučajevima moguće je postići da bude ispunjena nejednakost (58). 

(2) Primijetimo da je baš naša prva aproksimacija u, rješenje li- 
nearnog problema za € = 0, to jest rubnog problema 


—(auY=f,u0=0,u(D)=d. (59) 
Iz (40) i (43) dobivamo slično kao gore 


E3 


F(&)+C F(&)+C 


MOE a 19% 


luQ) — u (< 


0 


< | drI+ jat Koma lu (+2 10-01 (60) 
M sxs 0 


IC — C,| < K, (€) max [u (2)|. (61) 
0Zx<l 
Odatle izlazi konačno 


max |u(x) — u, (&)| SK, e nn Idi (|. (62) 


0gxql 


Vidimo da za e -> 0 rješenje u (koje također ovisi o e!) jednoliko konver- 
gira prema »nesmetanom« rješenju , linearnog problema. Tako je 
dokazano prvo svojstvo korektnosti. 


(27) Ovdje smatramo € čvrstim i puštamo da veličine 


Jo = max |/(&)|, d (63) 


o<tel 
teže k nuli. Tada je, kao što se lako vidi, 


max |u (x) — u, (x)| < const (d + fo)", (64) 
0<x<I 


pri čemu konstanta na desnoj strani ne ovisi o di fo. Tako smo provje- 
rili korektnost i u drugom slučaju. Činjenica da na desnoj strani ne- 
jednakosti stoji kvadrat »male veličine« d + fy znači da je aproksima- 
cija u, točnog rješenja osobito dobra. "Takva kvadratična greška je u 
vezi s postupkom linearizacije, koji ćemo opisati na kraju paragrafa. 


Iz razmatranih nelinearnih rubnih problema razabiremo da način rješavanja 
jako ovisi o konkretnom obliku diferencijalne jednadžbe. Posebno je važno ispitati 
može li se jednadžba bar »djelomično« integrirati (kao u primjeru 21.1) prije negoli 
se definiraju sukcesivne aproksimacije. Osim toga, važno je znati odrediti »bliski« 
problem, kojega rješenje uzimamo kao prvu aproksimaciju. Svakako, najvažniji 
način zamjenjivanja nelinearnog problema linearnim jest linearizacija. Promotrimo 
nelinearnu diferencijalnu jednadžbu 


D(u",u,ux)=0, (65) 
uz rubne uvjete : 
u(O) =0,u(]) =d. (66) 


Upotrebom Taylorove formule funkciju D za fiksiranu vrijednost varijable 4 = X 
aproksimiramo linearnom funkcijom 


3 Vo 
D (na, 022 )3> Ma) 2 D (91 729 1132 14) + EI (074 972 1152 994) (Mu — 1) 
i= 


oko točke (2%, !/> 119). Koristeći se tom formulom, aproksimiramo jednadžbu (65) 
jednadžbom 


DB (u8 (2), 6 (0 o (99,3) + (18 (6 (0) o (0,9) "+ 


+ 


06 00 (9) o (9) 9) 0 + (69) 


+ a (uo (2)> 40 (2) 40 (2) X) V =0. 
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Pri tome će 4 + v biti približno rješenje. Rubne uvjete za v postavljamo ovako: 
v(O) = —w(0)vO) =d—u (0) (69) 


(tako se osigurava da približno rješenje zadovoljava isteuvjete kao i točno). Tako smo 
dobili linearni rubni problem (68), (69). Kažemo da smo problem (65), (66) linea- 
rizivali oko funkcije ug. Taj se postupak može ponoviti tako da za dobiveno v iz 
(68), (69) provedemo linearizaciju oko “o + v itd. Tako se dobiva niz sukcesivnih 
aproksimacija za točno rješenje problema (65), (66). Ova metoda približnog rješa- 
vanja nelinearnog problema zove se Newtonova metoda zbog analogije s pozna- 
tom istoimenom metodom za rješavanje jednadžbe f(x) = 0 


21.5. Primjer 


U primjeru 21.1. imamo 


D (1),> 7922 1032 14) = N91 4 (M32 173) bo o 94) 12 +3 


se Na) Nn2 +J (2). (70) 


Upotreba formule (68) uz = 0 vodi na rubni problem 
—(avY=fvO)=0,vQ)=d, (71) 
uz a, (x) = a (0, x). Rješenje v ovoga linearnog problema identično je s rješenjem 


uy iz formule (43), ako stavimo a (u, x) = a, (x) + € u“. Provjeru te činjenice pre- 
puštamo čitatelju. 


21.6. Primjer 


Prva linearna aproksimacija problema (20), (21) dobivena linearizacijom u 
nuli daje u formuli (22) 


fi) =f(0,2) b() = —Z10, 29). (72) 


(Vidi posebni slučaj zadan formulom (31)). 


Primijetimo da se sukcesivne aproksimacije u primjeru 3 i zadatku 4.3 raz- 
likuju od sukcesivnih aproksimacija Newtonove metode, iako se prve aproksimacije 
mogu poklapati. 


Za Newtonovu metodu je bitno da se na svakom koraku rješava novi linearni 
rubni problem. 
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DODATAK 


1. Cauchyjev problem za obične diferencijalne jednadžbe 
Normalni sustav običnih diferencijalnih jednadžba prvog reda ima oblik 


Y =f1(90922 ++>9m X), i=1,2,..,N; (1) 
pretpostavljat ćemo da su funkcije fix 1 =1,2,...,n definirane na R" x [a,b] (a, 


beR, a <b). Pod rješenjem sustava (1) na segmentu [a, b] razumijevamo n- 
-torku funkcija y1, y2> <-> 9n €" (a,b), za koje vrijedi 


YI (x) =f, (M (x), Y2 (x), +++>9n (2), x)i1=12..n, (2) 

u svakoj točki x € [a, 6] u kojoj su one diferencijabilne*. Za proizvoljno x, € [a, b), 

No€R,2=1,2,..,n postavljamo za sustav (1) Cauchyjev problem: odrediti 
rješenje sustava (1) (na segmentu [a, b)) koje zadovoljava početne uvjete 

NK) =9o i=1, pe (E (3) 


Normalna obična diferencijalna jednadžba n-tog reda ima oblik 


PRE F9: 2019999, 093 (4) 


pretpostavljat ćemo da je funkcija f(Y,>92> .-.>9m X) definirana na R" x [a,b] 
(a,beR,a <). Pod rješenjem jednadžbe (4) na segmentu [a, b] razumijevamo 
funkciju y € 7" (a, b), za koju vrijedi 

99 (9) =f(09 (09 (2) ..097-0 (), 2) (5) 


u svakoj točki x e [a, b] u kojoj ona ima x-tu derivaciju. Za proizvoljno x, € [a, bl, 
No€R, £=1,2,...,n postavljamo za jednadžbu (5) Cauchyjev problem: 
odrediti rješenje jednadžbe (5) (na segmentu [a, b]) koje zadovoljava početne uvjete 


yt-D (x) =MNot=12.oN. (6) 
Problem (4), (6) može se svesti na problem (1), (3). Ako je y rješenje problema 

(4), (6) i ako je 
H=9"0i=12..m (7) 


* Diferencijabilnost u točki a odnosno 2 označuje postojanje desne odnosno lijeve derivacije. 


131 


onda vrijedi 
YA = (2), 92 Goss 9 (09), 2) (8) 


u svakoj točki x € [a, b] u kojoj y ima n-tu derivaciju ; iz toga slijedi da je x-torka 
(Y1>922 -.>Yn) rješenje Cauchyjevog problema 


NE=Nini=1,2...n— 1 9 = (90:92 nala) (9) 
Ji (%0) =yo1i=12.n (10) 


Obratno, ako je (Y1,92> ...»Yn) rješenje problema (9), (10), onda je y = y, rješenje 
problema (4), (6). 


Neka je funkcija f (91, 92» ...»9m X) definirana na R" x [a bl(aabeR,a< 
< b); kažemo da je ona po dijelovima neprekidna po varijabli x ako je neprekidna svu- 


da osim na R" X S, gdje je skup skokova S c (a, b) konačan, i ako za svako (y,, 
Y2> -.>Yn) € R" i svako x € S postoje limesi 


I(Y1>92> ..>9m X + 0) = lim FACIR N2> .+++>"m €), (1 1) 
UT sh <. >. 8 
ć > 
(9192 ++>I9m* — 0) = lim In» 1122 +++>Nm Š). (12) 
m9 i=1,2,..,1 
š<x 6x 


1.1. Teorem 


Neka su funkcije f(V12922 -->9m2)x 1=1,2,...,n definirane na R" x 
[a,b] (a, beR,a <b): po dijelovima neprekidne po varijabli x, sa skupom sko- 


kova S; neka su parcijalne dervacije i .43=12...,n neprekidne na R" x [a,b] 
j 
\ S) z neka je 


K = sup 


91:922:.s9n€R 
xe [a,b1/ 


(52) (PisY2> 99 2) | < 09% (13) 


£,j=1,2,...,n 


Tada problem (1), (3) ima jedno i samo jedno rješenje (y1>92> ...>Yn), Pri čemu su 
funkcijeyu it =1,2,...,n diferencijabilne na [a, DI \ S. 


Dokaz. Problem (1), (3) ekvivalentan je sustavu integralnih jednadžbi 
(2) = o + [f(91(8)292 (8) 091 (8), A6, x €[a,b), (14) 
0 


i=1,2,..)n 
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u ovom smislu: ako je (Y,, 22 ...>Yn) rješenje problema (1), (3), onda vrijedi (14); 
ako su funkcije y; neprekidne na [a, b] i ako vrijedi (14), onda je (Y1>92> <-+> In) 
rješenje problema (1), (3). Zaista, iz pretpostavke y€2f1(a,b) i1=1,2,..>n 
i jednakosti (2) i (3) neposredno slijedi (14). Obratno, ako su funkcije y,» = 1, 
2,...,n neprekidne na [a, 8] i ako zadovoljavaju (14), onda one zadovoljavaju (3), 
diferencijabilne su na [a, b] \ S i na tom skupu zadovoljavaju (2); iz toga za x € S 
dobivamo 


9 (X +0) =f(91 (992 (2), -->9n (2) * +0), (15) 
pajeme#F(a,b),i=1,2,...,n. Prema tome, dovoljno je dokazati da sustav 


(14) ima jedinstveno (neprekidno) rješenje. Definirajmo niz sukcesivnih aproksima- 
cija: 


yi = Jo (16) 
HU) =Y0 + [n (91 (£), 92 (8), 294. (6), £) dE, (17) 
k=0,1,2,.... 


Funkcije (16), (17) su neprekidne na [a, b]. Dokazat ćemo da ng y,%k =0,1,... 
jednoliko konvergira na [a,b]. Zaxela,b]\Sy,;NER1=1,2,...,n iz (13) 
slijedi (prema Teoremu srednje vrijednosti) 
n 
Ifa (F1 F2 +.>9m x) ZI (Y1>92> ++>9m x| < KX 19% —y| A 


< Kn Ka li — ul. (18) 


Neka je x € [a, bl, x & xo;5 iz (17) i (18) za & > 1 dobivamo 
ly" (0) — (< | Lfe 2 (2), 9 (E), &) — 
— 1 (sto 2 (6): . sMO (6), €)| di < (19) 


<K e zik <.) BV (€) — 171 (6) | e-%(€-x0) ga(€—x0) de, 


mr e)=peleKku. max (ly) = 
zh la, s 
yit-21 (£)| e-eE-x0)) = (20) 


|yP+t (x) — yi (x) | e-e(x-x0) < 


< max |yP(&) — yi (Bl e-e6-a0, 1) 
i=1,2,..., 


šela, Dj 
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Nejednakost (21) dobiva seizaxe[a,b]l,x < x,, pa za x € [a, b] imamo 
Lo" (8) — g (0) e-ate—aol < 


< max [of (B) — ofo" ()]e-elt=a0). (22) 
i=1,2,. 


čela, ši 


Iz toga indukcijom dobivamo 
lo" (8) — yi (0)| eele—sol < 


k_1 
< (2) max (Dy"() — 919 (8) | e-ule-sal), (23) 
Eat 


Neka je a > Kn. Iz (23) slijedi 


PŽ dot) — gi Bl erat < max (li"() — (24) 
re EEK 
e (Kn\'-" 
= Bljećem) 2 . — max (|P! (D— gl (£)| e—eit—*ol) ' 
k=i \q i=1,2,...>1 1 Kn 
žela,bl Bisu 
pa red 
X e-ale=a0l (yi (a) — yf" (x) (25) 
konvergira jednoliko na [a, b]. Parcijalne sume tog reda čine niz 
e-%lx-xoly(x), k =1,2,... (26) 


koji konvergira jednoliko na [a, bl, za svako z= 1, 2,...,n. Lako se vidi da isti 
zaključak nE izanizyr CG) k=.4 2,04 ZbOg neprekidnosti funkcija y!“, 
k=1|,2,..., neprekidne su na [a,bli funkcije 


yn (x) = lim a (Q)i=12..,n (27) 


Dokazat ćemo da te funkcije zadovoljavaju (14). Koristeći se još jednom nejednakošću 
(18) dobivamo 


[G, (YN (£), *++>Yn (6) £) JI (9 (8), .. M pa (E), Ehdč |= 


<G—a)Kn. max ba (0) — yi (&)l: (28) 
Sa [3 a id 
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desna strana u toj nejednakosti teži nuli za & —> oo pa zaključujemo 
| FS (8), 98! (8), E) AE > [f1(1(Ž), 99 (E), AE > 00. (29) 
ET) 


Uzimajući sad u (17) limes k —> oo, dobivamo jednakost (14). Pretpostavimo da 
osim rješenja (27) sustav (14) ima još neko rješenje 21, 22, ...» Zn» to jest da vrijedi 


ZI (x) = Jio + [4 (ži (£), ...> Šn (6), ć) dč, x€e [a, b), 1 = 1, Žarko n, (30) 


*o 


pri čemu su funkcije 24, € =1,2,...,n neprekidne na [a, 6]. Iz (14) i (30) slijedi 


l91() — 26091 = [ 191 (8) 99 1(8)8) —Ji(2a (8), > 2 (8), LE. B1) 


*o 


Iz toga slično kao u dokazu nejednakosti (22) dobivamo 


194 (3) — 2 (0)| e—sle—sol < SE mak (lpu(E) — m (B)le-at-e0, GB) 
eb &) 


ME JONG) EO gm nje (33) 


i= 
nE bI 


2 max Glan (E) — a, (&)| e-elf-x0l); 


i=1,2,.. 
čela, 6] 


zbog Kn/a < 1 zaključujemo da je y; (x) — & (x) = 0 za svako x € [4,0]. 
Za problem (1.4), (1.6) slijedi iz teorema 1.1. ovaj rezultat. 


1.2. Teorem 


Neka je funkcija f (Y1>92> --->Yn» X) definirana na R" x [2b] (abeR,a< 
< b)1 po dijelovima neprekidna po varijabli x, sa skupom skokova S; neka su parcijal- 


ne derivacije soi =12,...,n neprekidne na R" x [(a,b]'\ S) rnekaje 
j 


K= sup 


91:922...>InER 
xela,bINS 


of 
(52) (Y1> D92> +++>9m x) < o, (34) 


j=1,2,...,n 


Tada problem (4), (6) ima jedno 1 samo jedno rješenje, pri čemu funkcija y ima n-tu 
derivaciju na [a, b|\ S. 
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2. Linearne jednadžbe. Liouvilleova formula 
Sustav (1.1) je linearan ako je 


fi (Yi YV2> +...) Ym x) o ay (X)9 + b, (x), (1) 


gdje su a;; i b; funkcije na [a, b). Prema tome, normalni linearni sustav prvog reda 
ima oblik 


VE 2, đu (2) 9; +, (x), = 1/2,..3N. (2) 
= 
Iz teorema 1.1 dobivamo za sustav (2) ovaj rezultat. 


2.41. Teorem 


Ako je ai» bje#(a,b), tj=12,...,n, Cauchyjev problem za sustav a 
ima jedno i samo jedno rješenje (y1, Y2> . pa. pri čemu su funkcije Yi=1,2,. 
diferencijabilne u svakoj točki x e [a, b) u kojoj su funkcije ap x b3=1,2,. 
neprekidne. 

Sustav (2) je komogenako jeb; =0,1=1,2,...,n, inače je nehomogen. Ako rje- 
šenje (Y1>92> .+>Yn) Romogenog sustava 


= >4()y» i=12..,n (3) 
ij= 


zadovoljava homogene početne uvjete y; (0) = 0, i = 1, 2,...,n, ondajey; =0 
(to jesty,;(x) =0 za svako x € [4,b),f =1,2,...,n). Ako su 


(91,923. LA k=1,2 1 (4) 


rješenja sustava (3), onda je i njihova /inearna kombinacija s proizvoljnim koeficijen- 
tima eR, k =1,2,..., 7, to jest funkcija 


r 


2% 01 yi SM. >y), (5) 


također rješenje sustava (3). Rješenja (4) su žinearno nezavisna ako linearna kombi- 
nacija (5) iščezava samo u slučaju & =0,k = |,2,...,n; u protivnom rješenja (4) 
su linearno zavisna. Ako su za neko x, € [a, b] vektori 


I (0), 9" (20) «> . ade 1 (0)) k = 1, 2, +>7 (6) 


linearno zavisni, onda su i rješenja (4) linearno zavisna; drugim riječima, ako su 
rješenja (4) linearno nezavisna, onda su vektori (6) linearno nezavisni za svako X, € 
€ [a, b). Zaista, ako su vektori (6) linearno zavisni, onda postoje brojevi cx,k = 1, 
2, ...,7 od kojih je bar jedan različit od nule, takvi da je 
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r 
ža E" (05 92" (Ko) +9. (X0)) = 0. (7) 
Neka je 


=>4 m. (8) 
k=1 


Tada je (Y1>92> ...29n) dame sustava (3), a zbog (7) ono zadovoljava lane: 
početne uvjete y; (x0) =0,;=1,2,...,n;iztogaslijedidajey;, =0,2/ =1,2,. 


Svaki skup od » linearno nezađlenik rješenja 


Kik 


(yi »92>»- «ra == ...,H (9) 


sustava (3) je fundamentalan skup (sustav) rješenja. Dokažimo da fundamen- 
talan sustav postoji. Neka je 


(a9, a, ..,a9)k=1,2..n (10) 


proizvoljan skup od nn linearno nezavisnih vektora. Neka je (y%, 9%, ...,»%) rje- 


šenje sustava (3) uz početne uvjete 


Ne) =a) = et (11) 
Tada su vektori 
(st (£9) 92" A as (X0x £=1,2,...)1 (12) 


linearno nezavisni, pa su i rješenja (y1', 92), ..> 99) k=1,2..>n linearno ne- 
zavisna. 


Ako je (9) fundamentalni sustav rješenja, onda je opće rješenje sustava (3) 
dano formulom 


(912922 ...>9n) — 2% (95 Ny, . sate (13) 


gdje su c, proizvoljne konstante; drugim riječima, za svako rješenje (Y,; Y2> ...>Yn) 
postoje brojevi cx, k = 1, 2, ..., m takvi da vrijedi (13). Zaista, ako je (y1> 922 .-->9n) 
rješenje sustava (3) i x € [a, b], onda se vektor (y, (X0):92 (X0)> -.. 9n (X0)) može 
razviti po bazičnim vektorima (12); neka su cok=1,2,..,n koeficijenti tog raz- 
voja: 


(Yi (X9)> Y2 (X0), *<<>Yn (x0)) > Ck (yi, EVA (X9)> «<. ge (X0)). (14) 
Tada rješenje (91 922 .-.,9n) i rješenje 
n 
ž (992 +399) (15) 
k=1 
zadovoljavaju isti početni uvjet pa se podudaraju, to jest vrijedi (13). 
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Neka su (9) rješenja sustava (3); promotrimo matricu kojoj su ta rješenja stupci: 


Kr BEBE 
DOJE + + g 
LPE, A (16) 
LJA- 9 osi VR 


Determinantu matrice (16) označujemo s W i zovemo determinanta Wronskog 
ili Wronskijan rješenja (9). Iz prethodnih razmatranja slijedi ova činjenica: Ako 
je (9) fundamentalan sustav, onda je W (x) # 0 za svako x € [a,b]; ako su rješenja 
(9) linearno zavisna, onda je W = 0. 


2.2. Teorem 


Neka je W Wronskijan nekog fundamentalnog skupa rješenja sustava (3) i neka 
je x, € [a,b]. Tada vrijedi Liouvilleova formula 


dsoa 
WG)=W(ex)e (17) 
gdje je 
S (x) = an (18) 


Dokaz. Neka je (9) fundamentalni sustav i W njegov Wronskijan. Funkcija W 
je diferencijabilna u svakoj točki x € [a, b] u kojoj su diferencijabilne funkcije (9); 
izračunat ćemo W' (x). Neka je W, determinanta matrice koja se dobiva iz matrice 
(16) tako da se njezin 7-ti redak zamijeni retkom 


(0) OP), +9 (GY). (19) 


Tada prema pravilu za deriviranje determinante imamo 


W'G)= ZM). (20) 
Iz jednakosti 
(yi) = žut i=1,2,..,n (21) 
slijedi 
(98) OP) o (S) = ža (9,92, o 9), je 


pa zaključujemo da je redak (19) linearna kombinacija redaka matrice (16). Budući 
da se determinanta ne mijenja ako se nekom retku pribroji linearna kombinacija 
ostalih redaka, to je 
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W, (x) = au (2) W (29). (23) 
Iz (20) i (23) dobivamo 


W' (x) = S(x) W(2). (24) 
Iz (24) slijedi (17). 


Normalna linearna jednadžba n-tog reda ima oblik 
yo=%0()y%-9 +4(29, (25) 
j=1 


djesu c; i d funkcije na [a, b]. Obično se umjesto (25) promatra jednadžba 
d0(2)y9 +4,(9"-9 +... +a,(My=1(%) (26) 


gdje su a9, 4,> ...» G, 1f funkcije na [q, b] i ao (x) & 0 za svako x € [a,b]. Iz teorema 
1.2 dobivamo za jednadžbu (26) ovaj rezultat. 


2.3. Teorem 


Akojea,e€#(a,b),i=0,1,..,nbeEX%(a,b)iakoje 
min [ao (x)| > 0, (27) 


xe[a,b] 


Cauchyjev problem za jednadžbu (26) ima jedno i samo jedno rješenje 1 ono ima n-tu 
derivaciju u svakoj točki u kojoj su funkcije api =1,2,...,nif neprekidne. 


U daljnjem pretpostavljamo da su zadovoljeni uvjeti teorema 2,3. 
Uvodeći prema $ 1 funkcije 
=a =1,2.;8 (28) 
jednadžbu (26) svodimo na sustav 


N=Yinni=1l2..n—1, (29) 


si = — zoglen 6991+ aga (992 e +01 (0) ga — 1609). (30) 


Jednadžba (26) je komogena ako je f = 0, inače je nehomogena. Ako rješenje 
homogene jednadžbe 


20(K)yP +a()y"-0 +... +40(Qy=0 (31) 


zadovoljava homogene početne uvjete y(x0) =y (x0) =... = y9%-2 (x) =0, 
onda je y = 0. Ako su 
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sa k=N%or (32) 


rješenja jednadžbe (31), onda je i njihova linearna kombinacija s proizvoljnim. koe- 
ficijentima 6; € R, 4 =1,2,...,r, to jest funkcija 


Nag (33) 
k=1 


također rješenje jednadžbe (31). Rješenja (32) su linearno nezavisna, ako linearna 
kombinacija (33) iščezava samo u slučaju 6, =0,k =1,2,...,75 U protivnom rje- 
šenja (32) su Zinearno zavisna. Rješenja (32) su linearno nezavisna onda i samo onda 
kad su linearno nezavisna odgovarajuća rješenja homogenog sustava (29), (30). 
Svaki skup od # linearno nezavisnih rješenja. 


ELE =1 2.01 (34) 


jednadžbe (31) je fundamentalan sustav rješenja. Iz korespondencije jednadžbe 
(25) i sustava (29), (30) slijedi egzistencija fundamentalnog sustava. Ako je sustav 
(34) fundamentalan, onda je opće rješenje jednadžbe (31) dano formulom 


v=žaym, (35) 


gdje su cy proizvoljne konstante; drugim riječima, za svako rješenje y jednadžbe 
(31) postoje brojevi c;, R = 1, 2, ..., ntakvi da vrijedi (35). 
Neka su (34) rješenja jednadžbe (31); promotrimo matricu 


| yit y121 .. tal 
[2y [21)! 2. (VLRIY 
ga i a) a (36) 


(y121)"- 1) (9122) md)... (gl) 1) 


Determinantu matrice (36) označujemo sa W i zovemo determinanta Wronskog 

ili Wronskijan sustava (34). Wronskijan sustava (34) jednak je Wronskijanu skupa 

odgovarajućih rješenja sustava (29), (30). Iz toga slijedi ova činjenica: Ako je (34) 

fundamentalan sustav, onda je W (x) #£ 0 za svako x € [a, bl; ako su rješenja (34) 
linearno zavisna, onda je W = 0. Iz teorema 2.2 slijedi za jednadžbu (31) ovaj 

rezultat. 


2.3. Teorem 


Neka je W Wronskijan nekoga fundamentalnog sustava rješenja jednadžbe (31) 
i xy €[a, b]. Tada vrijedi Liouvilleova formula 


W (x) = W (x) e a (37) 
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